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Vol, 03, PART 7 1 July 1950 No. 367A 


Theory of Longitudinal Photons in Quantum Electrodynamics 


By SURAJ.N. GUPTA * 
Dublin Institute for Advanced Studies, Dublin, Eire + 


Communicated by W. Heitler; MS. recetved 20th Fune 1949 and in final form 
13th December 1949 


ABSTRACT. ‘he radiation field is quantized by introducing four types of photons— 
two transverse, one longitudinal, and one scalar. ‘The scalar photons are treated by using 
an indefinite metric, and it is found necessary to modify the usual supplementary condition 
slightly. ‘The present theory offers a justification for the symmetrical treatment of the four 
components of the electromagnetic potential, recently applied by a number of authors, and 
proves to be very convenient in applications. The results of physical interest, however, 
are the same as obtained from the ordinary formulation. 


Si) INTRODUCTION 


N quantum electrodynamics usually the longitudinal part of the electromagnetic 

] field is separated from the transverse part and replaced by the Coulomb 
interaction, which not only destroys the relativistic appearance of the inter- 
action but is also inconvenient in application. It has been shown by Majumdar 
and the present writer (1947) that the self-energy of an electron may be easily 
obtained by treating the longitudinal field in the same way as the transverse field, 
and more recently the same device has been used by Feynman (1948). In some 
applications, Schwinger (1949) and Tomonaga and co-workers (Tati and 
Tomonaga 1948) have also adopted a symmetrical treatment for the four com- 
ponents of the electromagnetic potential. But so far no satisfactory justification 
for this procedure has been given. We, therefore, give a theory of quantizing 
the radiation field by introducing four types of photons corresponding to the 
four components of the electromagnetic potential. ‘T'wo of these photons are 
transverse, one longitudinal and one scalar. Owing to the change in sign in the 
commutation relations for scalar photons, it has often been proposed to reverse 
the rdles of emission and absorption operators for scalar photons (see Dirac 1947). 
But this destroys the symmetry in the treatment of the four types of photons, 
and Jeads to difficulties (Ma 1949). Therefore we treat the scalar photons by 
using an indefinite metric and modify the usual supplementary condition. ‘The 
use of an indefinite metric for the scalar photons does not raise any difficulty of 
physical interpretation, For, as a consequence of the supplementary condition, 


* State Scholar of the Government of India. 
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only transverse photons are observable, while the longitudinal and the scalar 
photons exist only in virtual states. For simplicity, the familiar momentum 
representation is used, and to show the advantage of the present treatment some 
applications are also given. 

In the present paper the following notations will be used: Latin indices take 
the values 1, 2, 3, while Greek indices take the values 1, 2,3,4. A repeated Greek 
index is to be summed over all the four values. We shall often (explicitly) use 
the index 0 instead of 4, the relation between the two being A,=7A,. For 
convenience, a Greek index will be used also for the quantity «,, where «,, &5, & 
are the Dirac matrices, and «,=2 or a»=1. An asterisk denotes a complex 
conjugate or a Hermitian conjugate, while a dagger denotes an adjoint (for an 
indefinite metric the two concepts are different). We use natural units, i.e. we 
bake Ga —= Ie 


§2, QUANTIZATLION—OF THE RADIATION FIED 


For the reason already mentioned, the present treatment will differ from the 
usual one in the following two respects: 

(1) We take the electromagnetic 4-potential A, to be a complex 4-vector with 
three Hermitian space-like components and one antiHermitian time-like 
component, i.e. we take A, to be of the form 


A, =NOA®+NPA®+NOA@MLNOAM (1) 


where A®, A®, A®, A are scalars, AM, A®@, A® being Hermitian and A anti- 
Hermitian, and N/?, N, N©, N© are orthogonal unit 4-vectors, Ni), NN: NS 
being space-like and N{ time-like. It should be noted that we can choose a 
privileged frame of reference such that 4;=A®, 4,=A, and N= (OO Ors 


In general the components of N and N satisfy the relations 
(1) AJ) (2) \7 (2) 3) NT (3 0 0 
NON®+NONO+N@N®—NONM=3, (2) 


(2) We use an indefinite metric (Dirac 1942, 1943, Pauli 1943), i.e. we take 
the expectation value of an operator i(A,,) for a state © of the radiation field to 
be given by 

A) = OVA )O 0 = Oa) eee (3) 


where 7 is a Hermitian operator to be chosen later on in such a way that 4 is 
self-adjoint in the sense: 


A, Ay = Alt ee ee re (4) 


It follows from (4) that 7A, is Hermitian, and therefore the expectation values of 
A,, are real. ‘ 
The commutation relations for A,, are given by 


[AA*), 48) =@/2)8,.A(w—x'), (5) 


where A(x—.’) is the four-dimensional delta-function, and x denotes a space- 


time 4-vector. From. (1) and (5) we obtain forthe commutation relations of 
the A’s 


AM), AM] = (ACK), AMY] = (AM), Ax) = /2)A@—H)) 
[A(x), A(a’)] = ~~ (/2)A(a@ —»’), 
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We assume in the usual way that the radiation field is enclosed in a box of unit 
volume with periodic boundary conditions, and expand the A’s as follows: 


AO=E si (Z) (C,(k) exp [i(k . x — At)] +Cl(k) exp [ —i(k. r — ke)]} | 


A= > | (=) {Cy(k) exp [i(k.x — kt)] + Ci(k) exp[ — i(k. —A2)]}, 


rae (7) 
where GAN C; (ke iC ik y= = Ci k)e 2) in Aen (8) 
We then obtain from (6) 
[Ci(k), Ci(k)] =[C,(k), C2(k)] =[C3(k), C3(K)] =1 l eae (9) 
[Co(k), Co(k)] = —1, J 
all other commutators being zero. Thus, we have 
C,(k) = ae a Ck) == Tee) Sie wie (10) 


The Lagrangian density for the radiation field may be taken to be 


fp (eAnren 1 (OA, 0A 
e=- <1 (5) +3(F e =) \. iat (11) 


This gives in the usual way for the equation of motion and the Hamiltonian 
density respectively 


G0 rr en eee (12 a) 
and 
2 CA,\" 
= a { (3) +(curl A)? +(div A)? — (2) —(grad Ay ey a (125) 


Substituting the Fourier expansions for A,, from (1) and (7), using the commuta- 
tion relations (9), and omitting the zero-point energy, we get for the Hamiltonian 


{ He dv =Dk§ Ci(ke)C (ke) + Ch(k) Cot) + CH(k)Ca(k) — Ch(k)Co(k)}. wee (13) 


Corresponding to each value of k, the Hamiltonian describes four ‘radiation 
oscillators’, which may be treated in the usual way. 

We shall now confine our attention to a particular value of k, and drop the 
suffix k. We first consider the component NjVA“ of the radiation field. The 
quantity CjC,, which gives the number of photons associated with this component, 


has the eigenvalues 0, 1, 2,.... For a state containing m, such photons the 
eigenfunction $,(7,) and its adjoint #/(m,) are given by 

TAC YET Ai oes GHEE C30 el (221) igs Os a (14) 
with b()bi()= orig (15) 


The operators C, and Cj have the simple meaning of the absorption and the 
emission operators respectively, for 


Cyb1(m) = h(t — 1), Cidi(m) =(%41 + 1)¥9b,(m, +1). ...--- (16) 
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The matrix elements of C, and Ci, as defined below, are given by 
One a b}(m,)C\b3() a yO m,—1 
Ken nym, = b}(ny)Cibs(my) ma GER b, Ny—-1 j 
For the components NA and N\A), we have only to replace the suffix 1 by 2 
and 3 respectively in the above results. 

The treatment for the component NA is somewhat different due to the 
difference of sign in (9), which makes the use of an indefinite metric necessary. 
In this case the quantity C}C, has the eigenvalues 0, —1, —2,...., and therefore 
the number of photons is given by —C}C). In spite of this difference, according 
to (13) the eigenvalues of the energy of these photons are also positive. For a 
state containing m) such photons the eigenfunction 49(m) and its adjoint 4)(m9) 
are again given by 


Po(Mo) = (Mo! V?CN™; — Pi(Mo) = (Mo! CG, sev ees (18) 


with a different normalization rule 


Pi(Mo)Po(%)=(—1)™Singnge tte (19) 


Again, Cy and Cj are the absorption and the emission operators respectively, for 
Cobo(tto) = — Mol ?bo(tto'— 1), Chbo(tty) = (5 PL) eg 1). eee (20) 
The matrix elements of Cy and Cj, as defined below, are given by 
[Colne = Po) Cofo(%o) =( —1)""79"?8ing, ng—1 
[Col nsms = PUM) Ch po(™Mo) =( —1)""6"8m., ng—1 J 


Since, according to (3), we have 


D7 = £1 (My) P2(M2)h3(123)ho(Mo)h4 (12, )ho(I)h3(I3)ho(IMp) 
= Pj (My) P3(Me) P4(123)h)()h1 (my )bo( Mz) b3(mMs)ho(7Mo); 


it follows from (15) and (19) that the operator 7 is given by 


{NyNzN3zNo | 1) | MMMM, ) a ( am LOO pon Oe ace. ereyeys (22) 


This shows that states corresponding to odd values of mp occur with a negative 
normalization ©"7@. But this does not raise any difficulty of physical inter- 
pretation, as we shall see below. 


§3. THE SUPPLEMENTARY CONDITION 


We can choose the space-like unit 4-vectors NV), N°), N°) in an infinite number 
of ways such that they are orthogonal to each other as well as to a given INO’ 
In particular, we can choose them in such a way that for the k -component 

10 
of the field 
Ur, Ano 3 
k,Ni =k, Ni? =0, k, NO = —=k,,N io, 


pats case, the photons corresponding to the k,-components of N\A and 
NA® may be called the transverse photons, and those corresponding to the 
k,-components of NPA® and NA may be called the longitudinal and the scalar 
photons respectively. Thus, corresponding to each value of k,, we have two 
transverse one longitudinal, and one scalar photon. We note that the 
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k,-components of NYAY4+NEA®, N@A®, and NA are not transverse, 
longitudinal and scalar in the ordinary sense. ‘hey may be called generalized 
transverse, longitudinal, and scalar components respectively, and they reduce to 
the ordinary meaning in the special frame of reference in which N= (0,050.7); 

In order to account for the experimental fact that only transverse photons are 
observed, we have to impose some supplementary condition which will make 
longitudinal and scalar photons unobservable. We choose the supplementary 


condition to be 
0A, + 0A, 
| Oa) ale SRM, Sed en See oes (24) 


CA? OAS] : se ‘ 
where Ax and 3x | are respectively the positive and the negative 
i u 


frequency parts of 04,/dx,. This condition is weaker than the usual supple- 
mentary condition (04,,/0x,,)® =0, which is too restrictive to allow any state of 
the radiation field to satisfy it. 

Using (1), (7) and (23), we obtain from (24) for a pure radiation field 


(Ca CD Ee a ee (25) 
Now, a general state of the radiation field will be of the form 
Oe epee, (45 2a May Nose ees (26) 


where the A’s are numerical coefficients, and O(n, 15, 3,9) is a normalized 
eigenstate containing 7, and n, transverse, m3 longitudinal, and my scalar photons. 
Applying the condition (25) to (26), we get 

where 73, )=0, 1, 2, ...., and the coefficient A with a negative sufhx is to be 


taken as zero. It follows that the states of the pure radiation field allowed by the 
supplementary condition are 


D, =O(n,, No, 0, 0) 
DO, =O(n,, No, 1,0) —D(n,, 25, 0, 1) 
r= O71 1 NAOH ae 1s O71, He, Us 2) tC. Oban (27) 


1 +N rs Ae No, N3— 1p = 9; 


4, M2, N35 Mo— 


These states can be added with arbitrary coefficients, and the resultant state 
normalized. But, for all of these states, except ©), the normalization ®'® 
vanishes on account of (15) and (19). Therefore, a state of the radiation field 
satisfying (24) will be of the form 


D—O, + 4,0, +a,Do+.05, > tee ane (28) 
where the a’s are arbitrary coefficients, and 
O'd =O/®,. 


The terms with arbitrary coefficients in (28) do not contribute to the expectation 
value of the Hamiltonian of the radiation field. We call them the redundant terms, 
because they do not give rise to any observable effect. Thus, in a pure radiation 
field only transverse photons are present, and we may take 


Tees Bear R(20) 
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It should be noted that our supplementary condition (24) may still be regarded 
as the quantum analogue of the classical Lorentz condition 0A,,/0x, =—0, because 
according to (24) the expectation value of 0A,,/0x, vanishes, 1.e. 


bio 0), | ake (30a) 


OX, 


We may also obtain the quantum analogues of the Maxwell’s equations according 
to the present theory. The electric and the magnetic fields E and H are given by 


E=— = = gfad Aired = curl A. tee ener (30 d) 
From (12a), (30a) and (306), we obtain 
div H =0, | 
a +curlE =0, 


ot 
> 
i (div E)® =0, 


, (OE 
®* (— —curlH }) ®=0. 
ot ‘ 

§4. INTERACTION OF ELECTRONS WITH RADIATION 
So far we have considered a pure radiation field. The interaction of electrons 
with radiation may be treated in the usual way. However, to understand the 


meaning of our supplementary condition, the following remarks should first be 
noted: 


(1) The supplementary condition (24), being weaker than the usual one, is 
evidently compatible with the equation of motion of the system. ‘Therefore, 
if (24) holds for the initial state, it will also hold for all the subsequent states. 
A formal proof of its consistency may be given in the usual way. 

(2) As shown in the last section, in the absence of interaction our supple- 
mentary condition implies the absence of longitudinal or scalar photons. 

(3) In the presence of interaction the supplementary condition still holds, 
but (25) does not follow from it because now the field equation takes the form 
[_?A,,= —j,,and we cannot use the expansions (7) with their time factors. There- 
fore, in the presence of interaction the supplementary condition does not imply 
the absence of longitudinal or scalar photons. 

In most cases of practical interest, the initial and the final states may be taken 
to be interaction-free, while the interaction may be treated as a perturbation, 
which is switched on adiabatically in the initial ‘state and switched off in the 
final state. In such cases, it follows from the above remarks that longitudinal 
or scalar photons do not exist in the initial state, and though they may occur in 
intermediate states, they will again disappear in the final state. We have also 
verified this situation by direct calculations in a number of cases. In fact, it is 
found that the effect of the emission of real longitudinal or scalar photons is merely 
to change the coefficients of the redundant terms in the final state. But since the 
redundant terms have zero normalization, the transition probability for the real 
emission of longitudinal or scalar photons vanishes. In other words, the 
probability of the emission of a real longitudinal photon is cancelled by the 
‘negative probability’ of the emission of a corresponding scalar photon. 
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The Hamiltonian H for a system of electrons and the radiation field is given by 
Veur Te Rea EL mo Merce ee eee (31) 


where H, is the sum of the Hamiltonians of the electron and the radiation field, 
and H is their mutual interaction given by 


H=—| A,jyde, j= eho Aga (32) 
The expectation value of His 
GRU s ine Beep hp ts ene (33) 


Using the ordinary perturbation theory or the method of canonical transformation, 
it may be easily seen that 


Hy =Hygy Hy = bees Me oye Vinee haan a (34) 


where H, and H, are the expectation values of the first and the second order 
interaction energies respectively for the transition from a state @ to a state b, and 


Hy,=OnH®,, Nag=On®, Eyg=O,H iP. «+--+: (35) 

Now, we have 

ee 33 Aya =o» where TiO Ome hee Patek, (36) 
Since, according to (22) and (29), 7,, =1 for a final state, and n=, =! for any state, 
(34) reduces to 
ighiva le 
E,—E,. 
Thus, in spite of the use of the indefinite metric, 7 does not occur in the pertur- 
bation formulae. 


Using Fourier expansions similar to (7), but without their time factors, we 
obtain from (32) in the usual way 


Go Oe ute 


bar 


ClO." eles Meu ae.: (37) 


vay: 


vel =e e (=z) LICE) ip, +k, pp Ir [Ci(kK)]ia5p,,, ppt} Nn (use utea) 
+ {[Co(k)] ba Dy +k, pp 7s [Ci(K)] ae Pp NO (ua Mead) 


where p,, and p, are the momenta of the electron, and u,, and u, are its Dirac 
amplitudes, while [C,(k)]j,, and [Ci(k)];,, are given by 


[Ci(k) ia =P,C(k)Oq =a"? ong, nav 
[C}(q) oq =P C](K) Og = M9" Sng, mgt 
[Co(k)]in =P5Co(K) Pu = — Ma Ong, na-D 
[Cin =P, CHK)O a HM ny, mgt 


The occurrence of a minus sign in the matrix elements of [Co(k)]i0 is very important, 
and is due to the fact that A is antiHermitian. 

-. Ag already mentioned, according to (2), there exists a special frame of reference 

in which A,=A® and A,=A®:. Since there is no physical reason to choose this 

special frame of reference in any particular way, we may choose it arbitrarily 
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according to convenience. ‘The most convenient choice, of course, will be to 
take this special frame of reference to be the same as that in which the calculations 
are to be carried out. In that case (38) reduces to 


Qn TERE t 
Hha= = 6 | (FE) LLCO py rp FLU Pp nye ete) 
—{[Co(k)]indpy +4, pp tLCO) ind py, pps (Me%o%a)]- +++ 2+ + (40) 
Using these matrix elements, we shall now consider some applications. 


§5. INTERACTION OF TWO ELECTRONS 

We consider the interaction between two electrons whose initial momenta are 
p, and p, and final momenta p{ and p;. Let the Dirac amplitudes of the electrons 
be wu, and w, in the initial state, and uw and uj in the final one. Neglecting the 
exchange effect, the following two transitions are possible: 

(1) The electron p, emits a photon k=p, —p; going to the state pj; the 
electron p, absorbs the photon k going to the state pj. 

(2) The electron p, emits a photon —k=p, —p going to the state p3; the 
electron p, absorbs the photon —k going to the state pj. 

The total interaction energy is given by 


Hi, Hi, 
a Fe alert aus ee EG Gd 1 
s “EE, G ) 


where i, n and f refer to the initial, intermediate and final states respectively. 
We add the contributions of the four types of photons, bearing in mind, according 
to (39), the minus sign in the matrix elements of the absorption operator for a 
scalar photon. We thus obtain 


Te 27e? 1 1 y x y® i: 1* 

= FE (gaeca + orca) {3 , uted a) — (uae) ate) 

Since the energy in the initial and the final states must be the same, we have 
E, —E, = —(E, —E3) =e, say. 


ATRL II cas ve 
Hence ie 2 (tip cc 1j the 206, tks) ye Ot we eae (42) 


which is Moller’s well-known result (cf. Heitler 1944). 
We further see that the contribution of the longitudinal and the scalar photons 
alone is 


/ 4 os y* * ee * k 
P= 3 l(ua. nu,)(us’o. nu.) (uu, )(us"m)], =F, 


Acre” 
are (1; "t4y)(wg"te) 


477e? 


ay (u(k, ry —1)}. (Wy) (Ys"po). 


Therefore, the above exchange of the longitudinal and the scalar photons is 
equivalent to an interaction 


47re* 
V.= a exp {a(k, ry —1,)}. 
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Summing over all values of k, we get 
e 
Sly atta es ie aera 
eo | Dae se (a(k, tr, —1,)} dk 
— e 
[Poet 
which is the Coulomb interaction. 
It should be noted that if we use the more general matrix elements (38), we 
obtain for the Moller interaction 
Are” 
_—_ 
e?—k 


which, on account of (2), is the same as (42). 


r(2) 72 3) NB 0) O)\(4)* * 
(NON 4+NYNY 4 NON® —NON)) (uy otes)(Ue 0, Ul) 


§6. SELF-ENERGY OF THE ELECTRON 


We shall now discuss the self-energy of the electron, which is defined as the 
interaction energy when the electron is present minus the interaction energy 
in vacuum. To calculate the self-energy of an electron with momentum p, 
consider the following transitions: 

(1) The electron p emits a photon —k going to the state (p +k); the electron 
(p +k) re-absorbs the photon —k going to the original state p. 

- (2) An electron p and a positron —(p +k) are created with the emission of a 
photon k; the pair is annihilated with the re-absorption of the photon k. 

We notice that the first transition can occur only when the electron is present, 
while the second transition can occur only in its absence. Therefore, the self- 
energy is equal to the difference of the contributions of the above two transitions, 
and is given by 

W=> 2ne* | s* (Uj %,t)(U" ot) SC (wp%ytt)(U" atl) | ’ 
ee ky—-E-k —fy—-E=R || 


where u, is the Dirac amplitude of the positive energy electron p, and u and uw 
are the Dirac amplitudes of the electron (p +k) in positive and negative energy 
states respectively; S+ and S~ denote summations over the spin directions of u 
and u; k& is the energy of the virtual photon, E={p? +(p +k)?}4?, and 
E,={u2+p?}#2. After the usual simplifications we get 


e& dk[p2?+2u2—-EE,+(k.p) ee 


Caen, Ek E+k-E, E+k+£, 
Rare (44) 
For the special case p=0, (44) gives the well-known results of Weisskopf 
(959): : 
Seas |= wen) (45 
w=" | aka + aE >) 


For the general case, in order to obtain the result in a covariant form, we 
evaluate (44) by a method similar to that used by Pauli and Rose (1930 )am We 
‘first carry out the transformation from the variables k to the variables ¢, v and w, 
where ¢ is the azimuthal angle around the axis parallel to p, and 


Lape EI ee ae (46) 
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Computing the functional determinant, we obtain for the elements of volume the 
relation 


ak eh 
— == fh IR i ithe FNP ect 47 
Ek — 2p dv dw dd (47) 
The limits of integration of the new variables are 

0<¢<2z, ] 

UO DE ae ee! teats te), eee (48) 
w —p? —p 
(ue +p ye <w<o. 


The above limits of integration of v are obtained from the inequality 
; 1 
(cw —p?— p2)? = 7 (k. p)*(o —w)*<p'(v —w)? 


or, (w? —p2)v2 —2ep2v + (ut +p4 + 2p%u2—p'w2)<0. sean. (49) 


Transforming (44) to the new coordinates, we obtain the simple expression - 
2p?w 
ds 
=e samp | aa b(o+ a) 


~Hew(?, (3 1 50) 
TE a ie ae(2 3) Sioarhets (50) 


where e =| — pe ee eee Gi 


It should be noted that the calculation of the self-energy involves a summation 
over states containing a virtual pair of electron and photon with momenta (k +p) 
and —k and energies E and k respectively. Since (k+p,/) and (—k,k) are 
4-vectors, their sum (p, £ +A) is also a 4-vector, and hence « is a Lorentz-invariant 
buantity. The self-energy (50) may, therefore, be expressed as 


Wares toe = bu | y Bi do ee (52) 


where dy, though divergent, is Lorentz-invariant. This shows that the effect of 
self-energy is equivalent to a change in the rest-mass of the electron. ‘The same 
conclusion has also been arrived at by Schwinger and Tomonaga. 


§7. CONCLUSION 


For the purpose of applications, the main results of the present theory may be 
stated as follows: 


(1) The electromagnetic 4-potential A,, may be quantized as 
4-5] (FZ) {C:) exp [ie »)] + C10 exp [ i(k. x), | 
k 
2 :; +o 
Bis Hl (Z) {Co(k) exp [2(k . r)] + Ci(k) exp[ —i(k.r)]}, 
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where C; and C) are the absorption operators, and C! and C} the emission operators 
for four types of photons of positive energy. 


(2) Ordinary perturbation formulae may be applied, taking the above absorp- 
tion and emission operators to be the matrices 


[CG] an wi [Co] age at! 70) 
[Cl] wm = [Ci] of a nv? 


m,n-1) 


Mm, Le 3 


(3) All the four types of photons may be present in virtual states, but longi- 
tudinal photons (corresponding to the component of C,(k) along k) and scalar 
photons do not exist in a pure radiation field. 


Using the above results, we can easily treat the interaction of the radiation 
field with particles of any spin. 


Note added in proof. In the present paper we have used the ordinary 
representation for the treatment of interaction. The extension of the present 
theory to ‘interaction representation’, which requires the generalization of our 
supplementary condition, has been given by K. Bleuler in a forthcoming note 
in the Physical Review. 
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ABSTRACT. When the energy of «-particles emitted by radioactive nuclei is compared 
with the energy calculated from the usual semi-empirical formula for nuclear binding 
energies, systematic variations ranging from 4°5 Mev. to —1-:0 Mey. are found, which 
exhibit in a very clear way the special importance of configurations containing 82 protons 
and/or 126 neutrons. The general features are readily explained in terms of the shell 
model. The analysis can be used to predict unknown a-energies. Calculated a-energies 
are tabulated for the heavy nuclei (Z >78) and for nuclei in the rare-earth region. 


$1. INTRODUCTION 


T is well known that the general trend of the mass-defects of nuclei can be 
] satisfactorily represented by a formula of the type 


M(A, Z)=ZM,+NM,,—«A+B(N—Z)2/A +yA28 4 8224-48, 


where Z is the number of protons, N of neutrons, in the nucleus, d=Z+N is 
the mass-number, M(A, Z) is the mass of the nucleus, M,, M,, are the masses of 
the proton and neutron, and «, B, y, 5 are constants. Such a formula represents 
an average variation of nuclear masses, and does not take into account the 
individual characteristics which arise from special groupings in the nuclear 
structure. In particular, it does not, as it stands, represent the fact that there are 
systematic variations with the parity of A, Z and N, nor the special importance of 
‘magic numbers’. On the other hand, it is sufficiently accurate to give a very 
good account of the beta-stability of nuclei. The difference between the actual 
mass of a nucleus and the mass calculated from the formula is in fact rather small, 
and shows systematic variations reflecting the individual characteristics of the nuclei, 
which presumably depend on combinatorial properties of A, Z and N, such as 
the completion of closed shells. Unfortunately only the masses of the lighter 
nuclei are known with sufficient accuracy to make a significant comparison 
directly, for the deviations are of the order of a few m.M.U., and the precision with 
which the masses are known decreases below this limit with increasing mass- 
number. 

The energies of the alpha-particles emitted by radioactive nuclei, on the other 
hand, are known with very high precision and furnish a sensitive comparison with 
the calculated values of the difference M(A, Z)— M(A-—4,Z—2). The energy 


of disintegration (i.e. energy of the «-particle plus energy of recoil) in the reaction 
(A, Z)>(A—4, Z—2) +a-particle is given by 


E,(A, Z)=[M(4, 2) MA 4,7) i (2) 
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The difference between the observed value of £,, and the value calculated from (1) 
we shall call the excess and denote by A(A, Z): 


A(A, Z) - Le cots ies Teac: Cot oc (3) 
It reflects the characteristics of the two nuclei involved in the disintegration. 


For the purpose of the present calculation we shall take (1) in the numerical 
form 


[M(A, Z)-44M(4, 2)]c? = —7-0A4 +0-39(N — Z) +19-6(N — Z)?/A 
14-047? ¢0:6177A-t8 Mev 1.5 (4) 
These coefficients differ slightly from those given in the literature. They were 
chosen because tabulated values of M(A, Z)—+4A(4, 2) were available for them, 
having been computed in 1945 for use in the Montreal Laboratory of the Canadian 
National Research Council, they then being the best values available; and it did 
not seem worth while recalculating them with newer values. ‘The conclusions 
from the present analysis are essentially unaffected by the particular choice of 
constants. 
S28) EE AWAYaN UCI El 

The results of the comparison for the heavy radioactive nuclei (Z>78) are 
shown in Table 1. The nuclei are tabulated according to Z and N. For each 
nucleus there is listed: 

(a) Its chemical symbol (or classical radioactive symbol where appropriate) 

and mass-number: e.g. U1238. 

(b) The disintegration energy F,,,,, where known (in roman type): e.g. 4-251. 

(c) The calculated value, E,,.9;, (in italic type): e.g. 5-249. 

(d) The excess A=E,,,,— Exec (in Clarendon type): e.g. — 0-998. 

Energies are in Mev. Calculated values are believed accurate to two or three 
key., and better in the middle part of the table. ‘They are quoted to this accuracy 
so that they may be used for extrapolation, should occasion demand it. Experi- 
mental values are from Seaborg and Perlman (1948), Feather.(1949) and Perlman, 
Ghiorso and Seaborg (1950), and are corrected for recoil; except for the «emitting 
isotopes of gold and (?) mercury, which come from Thompson, Ghiorso, 
Rasmussen and Seaborg (1949). The energies for RaE 210, Pu241 and Am 242 
have been inferred from closed decay cycles, and that of At 212 from the life-time, 
following Broda and Feather (1947) and Perlman, Ghiorso and Seaborg 
(1948, 1950). The gold and mercury isotopes have been assigned, arbitrarily, 
the mass-number 186. 

The table shows some remarkable regularities: 

1. For the seventy-two nuclei with Z>84, N>128 (those in region A of 
Figure 1), A decreases reasonably smoothly with increasing mass-number, starting 
at 4-525 mev. for ThC’ (Z=84, N=128), decreasing at first by 0-33 Mev. per unit 
of mass-number, and then levelling off at approximately —1-0mev. for the 
heaviest «emitters. The excess appears to be closely correlated with the mass- 
number, as may be seen from Figure 2, in which it is plotted against A for the 
nuclei in the range from A=202 to A=244. The points for region A all lie 
reasonably close to one curve. There is no such close correlation with Z, N or 
N—Z. Nor are the departures from the curve correlated in any obvious way 


with Z, N or N—Z. 
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2. There is a sharp break between Z=84 and Z=83, and a similar break 
between N=128 and N=127. In the ‘tram-line’ regions B, and Bs, corre- 
sponding to Z=83 and N= 127 respectively, six experimental points are available. 
The excess starts at about 2-8 Mev. and decreases with increasing mass-number 
with a slope somewhat smaller than in region A. 

3. There is another break between N=127 and N=126, and, although no 
experimental points are available in region D, to prove it directly, there is probably 
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Figure 1. The classification of nuclei according to regions. 
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Figure 2. Plot of the excess of observed over calculated alpha-energy against mass-number. 


a similar break between Z=83 and Z=82. ‘There are sixteen points, some of 
which are still doubtful, in region D, (Z>84, N<126), which furnish the evidence 
for the former break. 

The various regions have a definite significance in relation to a closed shell of 
82 protons and a closed shell of 126 neutrons. In region A, both parent and 


£ 
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daughter nucleus have at least 82 protons and at least 126 neutrons. In region B, 
the parent has more, and the daughter less, than 82 protons, while both have at 
least 126 neutrons. In B, both have at least 82 protons, but the parent has 
more, and the daughter less, than 126 neutrons. And so on in a way which is 
clearly seen in Figure 1, if one remembers that the daughter is the nucleus with 
two protons and two neutrons less than the parent nucleus, which is the one 
indicated. 

It will be seen from Figure 2 that the points belonging to region A clearly 
fall into one group, those in regions B, and B, into a second group lying close to 
a single curve, and those of regions C and D, fall into a third group, also lying 
close to a single curve. ‘The data for the neutron-deficient bismuth isotopes, 
which lie in region E,, have not been plotted, their mass-number assignment being 
uncertain, but it is probable that they belong to a fourth group, particularly as 
Bi 209 (also in E,) is apparently stable, from which it can be inferred that its excess 
is less than —0-3 Mev. 

It is proper at this point to say something of what can be inferred concerning 
the apparently stable nuclei. On the basis of the relation between energy and 
lifetime it has been assumed in the table that if the «-particles had an energy greater 
than 4-0 Mey. they would have been detected. This estimate is open to question, 
for it depends on a specific assumption concerning the radius of the daughter 
nucleus, which is rather uncertain for the nuclei in question, probably being 
considerably smaller than for the classical radio-elements. And for rare isotopes 
of an element, «-emission would only be observed if the specific activity were 
higher than is necessary for detection in an abundant isotope. It has not been 
thought worth while to make a detailed estimate for each isotope, however, on 
account of the vagueness of the conclusions to be drawn in this way. An upper 
limit of 4-0 Mev. has therefore been assumed throughout for the elements from 
platinum to bismuth. This gives a few upper limits to A in the region F. 
Positive information comes from the newly discovered «-radioactive isotopes of 
gold and (?) mercury, whose mass-numbers are not yet certain. The gold 
isotope is quoted (Thompson et al. 1949) as lying in the mass-range 185-188, and 
for the purpose of tabulation it has been listed under 186. 


§3. THEORETICAL DISCUSSION 


The general features described above find a ready explanation in terms of the 
nuclear shell model (Mayer 1949, Nordheim 1949, Feenberg and Hammack 1949). 
If we assume that 82 protons form a closed proton shell, and 126 neutrons a closed 
neutron shell, and that the binding energy of a nucleon inside a closed shell is 
greater than of a nucleon outside a shell, by about 2 Mev., we can understand the 
breaks between the regions A, B, C, D etc. The decrease in A with increasing 
mass-number in regions A and B is also easily explained (cf. Glueckauf 1948) if we 
suppose that, as more and more nucleons are added outside a closed shell, their 
mutual binding energy increases, until the irregularity arising from the proximity 
of the closed shell is smoothed out, apart from minor fluctuations. 

We may try to derive quantitative results from a quite naive application of the 
shell idea. Let us suppose that outside a closed shell the energy of a nucleus 
differs from the smoothed value (1) bya term proportional to the number of ‘loose’ 
nucleons, representing the lower binding outside the shell, together with a term 
which is quadratic in the number of loose nucleons, representing their mutual 
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binding; and that inside the closed shell it differs by a term proportional to the 
number of vacancies, representing the tighter binding inside the shell. The 
quadratic term will, of course, only be a good representation near a complete 
shell. We then get the following: 


(M— My.)c2 = «(A — 208) — $a(A — 208)(A — 209), Z>83, N>127; 
=<(N~ 126) e'(82—Z) =4a(N — 126 N 127) oe 
= ¢(Z—82)=<'(126— N)—4a(Z— 82\(Z 83), 7 8a ie 
Se 0s Za 82.2 N10" 


A constant term representing the discrepancy for the nucleus Z=82, N=126 
(Pb 208) has been omitted, as it contributes nothing to the energy differences. 
The terms $a(A — 208)(A — 209) etc. are proportional to the number of pairs of 
loose nucleons. It is assumed that in the interaction between loose nucleons it 
is immaterial whether they are protons or neutrons. 

This gives the following results for the excess: 


Region Excess 
A 4 —4a(A —210-5) 
Bie Be 3¢ +’ —3a(A —210) 
C (RaE 210) 2e +2e’—a 
D, 2e + 2e’ —Za(N — 127-5) 
D, 2e + 2<’ —2a(Z — 83-5) 
E,, E, e+3e’ 
F 4e’, 


If we take the slope of the smooth curve drawn through the points of region A 
in Figure 2 to determine a, we find the value a=0-083 Mev. The line through 
the points of the regions B has been drawn with slope 3a, according to (6), and it 
will be seen that the fit is reasonable. The values of e and «’ from these data are 
« =1:26Mev., and e’ = — 0-68 Mey. 

This predicts the value A=1-08 Mev. for RaE and the polonium isotopes of 
region D,, and 0-92 Mey. for the astatine isotopes. The agreement at this point 
is poor, but this is not surprising in view of the crudity of the assumptions. As will 
be seen later the discrepancy is readily accounted for when one makes more 
realistic assumptions. 

We might refine our asumptions in several ways. In the first place, we could 
introduce different coefficients ¢,, <«, and €,', €, for protons and neutrons. 
The good correlation between A and A in region A indicates, however, that «, and 
€, cannot in fact be very different. And the similar-behaviour of regions B, and 
B, indicates further that ¢,’ and ¢,’ are nearly equal. We could also introduce 
different coefficients @,,, dy, Ay, for the mutual interaction of loose nucleons. 
Again, the correlation of A and 4 indicates that a,,, and a,,, are not very different; 
and the relative slopes of the curves through the points of regions A and B indicate 
that a,,, is approximately equal to a,, and a,,,, though this last datum is not very 
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significant. We should also introduce a term depending on the parity of Z and 
N; this is particularly important in connection with regions B and E and the 
discrepancy just noted. Finally we might introduce quadratic terms inside the 
closed shell to represent the variation of the true smoothed mass from the calcu- 
lated smooth mass, and bilinear terms for nuclei having loose protons and an 
incomplete neutron shell, and vice versa. With so many parameters, remembering 
that there must also be minor fluctuations of a combinatorial character, it is 
not worth while to pursue the matter very much further until more reliable data 
are available for regions D, E and F. 

There are a few points, however, which may profitably be noted. It will be 
seen that region D, is represented in Figure 2 by points which fall clearly into two 
groups—those nuclei with 126 neutrons have a value of A around 0-7 Mey., and 
those with 125 neutrons or less around 0( + 0-2 Mev.). This is readily explained 
if one assumes that the last neutron completing the shell is less tightly bound 
(by 0-7—0-8 Mev) than the previous ones, though still much more tightly than 
neutrons outside the shell. This is in line with the recent observation that the 
y-n threshold of Bi209 is abnormally low (Collie and Parsons 1950).* 

It is commonplace that nuclei with unpaired nucleons systematically have 
higher energies than nuclei with all nucleons paired. If we assume that this effect 
is more marked with nucleons outside a closed shell than inside we can understand 
why it is that the points of region B (extrapolated to 4 =210) lie higher than 
predicted on the simple theory, namely, halfway between A and D. Further, 
the assumption that this is particularly marked when there is only one loose 
nucleon, as in the regions E, would explain the apparent fact that the excess for the 
neutron-deficient bismuth isotopes of region E, is not so low as would have been 
expected on the simple theory. 

We have said that there is no other obvious correlation of A than with mass- 
number. ‘This is not altogether true, for it can be seen that the actinium points 
tend to lie below the curve, and the plutonium points to lie above it. And if one 
looks carefully, one can detect other indications of regularities, but it is doubtful 
whether they are very significant. It may be remarked that the low value for 
U 233 may be due toan error in the experimental value, for there is a y-ray associated 
with the disintegration, whose energy ought perhaps to be added to the «-energy. 
The standard deviation of the points from the smooth curve is 0-15 Mev., but there 
is a tendency for greater scatter in the region around A=234, ‘This may be 
connected with the incidence of isomerism in this region. 

That the general behaviour of the excess is very little sensitive to the choice 
of numerical constants in (1) may be seen from the following argument. If we 
consider the change in E,,,,;, caused by a change da etc. in the constants, we can 
estimate it sufficiently accurately by writing 4c?(9M/0A)y_, for E,. This gives 

2 

ee (=) 9B pA My 4474-38 (1 = a) ee (7) 
The constants given by Rosenfeld (1948) differ the most of any published from 
(4), namely «=14-66mev., B=20-54Mev., y=15-4Mev., 5=0-602mev. The 
constants in (4) are, taking 28:2 Mev. as the binding energy of the «-particle, 
a = 14-05 mev., 8 =19-6 Mev., y=14-0 Mev., 5 =0-61 Mev. This means that, had 

* Note added in proof. The y-n thresholds of Pb 206, Pb 207, Pb 208 recently measured by 
Hanson and Palevsky (private communication from G. T. Seaborg) confirm this further. 
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we adopted Rosenfeld’s values, we should have obtained an «-energy lower, 
and so an excess higher, by 


Nay : : Z —1/3 
2:44 +3 76(~>") _ E 13 =” 0322 (1-3) | MEV. .-eees (8) 


The values of this, for three widely spaced nuclei, are, Sm 144: 2-11 Mev., Pb 208: 
2:38 mev., Cm242: 2-43mev. Although the discrepancies are considerable, 
they are nearly the same in all three cases. ‘The major discrepancy comes from 
the value of «, which gives a constant difference of 2-44mev. We may therefore 
regard the variation of A(A, Z) which we have evaluated to be a correct repre- 
sentation of the facts, but place no reliance on the zero with respect to which it is 
referred, as this is too sensitive to the choice of «, the mean binding energy per 
nucleon, which is not known sufficiently accurately for the purpose. 


64. RARE EARTH ISOTOPES 
The existence of a natural «-emitting isotope of samarium, and the recent 
discovery (Thompson et al. 1949) of artificially produced «emitting rare-earth 
isotopes, makes it desirable to extend the calculation to this region. ‘The 
incidence of u-radioactivity in this region is presumably closely connected with 
the closing of a shell of 82 neutrons. This is done in Table 2, which lists £,,.).- 


Table 2. E,,.;. in the Rare-earth Region 


Element Nd Pm Sor Eu Gdeee ibe DD yaeto Bee. Ahan, — Wo 
Z 60 61 62 63 64 65 66 67 68 69 70 


N 

98 200" - 2:43. 12:95 
97 PINES DA) jo 
96 2:30 a2iO seo 
95 206% "2754.23 Oils es"46 
94 1762-242 ZS nO, 
93 WAS = 194)" 92-422 SONS :3 eseo 
2 1145 1-64. 92-13" 226053005 Slee 5 
91 O83 S41 S83 2-8 ee ia ee Oem cal) 
90 02514 1-03 e153) © 2-012 4 2-0 54 OO Ome 
89 0-19 .0°71-°1:22 1°71 2:19 2°66 ~ 3-12 53°56 ~ 3-99) 4-44 
88 = 0-14 -0°39.- 0-91 2°41— 91°90) 2237 2584-9 328 = 68/2, at eee 
87 0°07> 0°60 =A TIS 16022-05225 563-00 G4 eS Onn > Omer) 
86 0:28) 0:802% 1:30) 1-79) 222653 2:72) 3-17 3-61 at-03 54 eee 
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83 O°S9 1-3 8 1S 865233952718 22 ee OOM One 4 Ono 
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Beta-stable isotopes are indicated in heavy type. Sm 1521s naturally radioactive, 
with a disintegration energy of 2-4 Mev., corresponding to an excess of 1-9 Mev. 
The study of the artificial radioactivities will no doubt give more data in the near 
future. 
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§5. CONCLUSION 


The advantage of subtracting the calculated from the observed energy of 
a-particles is that one then obtains a quantity which varies fairly smoothly with 
mass-number, independently of the element. ‘The sharp breaks between proton- 
numbers 84 and 83, and between neutron-numbers 128 and 127, and between 
127 and 126, are some ten times larger than the deviations from smooth behaviour, 
and are clearly real. ‘here can in fact be no doubt concerning the special impor- 
tance of a configuration of 82 protons, or of 126 neutrons. ‘he reasonable success 
of the analysis of §3, based on a very simple application of the shell model cannot, 
however, be taken as evidence for any very simple theory of the nature of the closed 
shells, though it may give one some confidence in the predictions concerning the 
regions where little or no experimental data are available. 

The fact that the deviations from the smooth curve for the nuclei with Z>84, 
N>=128 are of order 0-15 Mev. indicates that one is able to predict the energies 
of unknown «-decays in this region to about this accuracy. This, of course, 
is no improvement on the plots of Glueckauf (1948) and Perlman et al. (1948). 

The isotopes of lead (7=82) present a case of some interest. So far, no 
w-active lead isotopes are known. ‘The position of RaD (Pb 210), occupying in 
region D, the position corresponding to Po 210 in Dg, suggests that its excess is 
probably about the same (compare the similarity of regions B, and B,), namely 
about 0-7 Mev., and hence that it emits «-particles of energy about 4-1 Mev. 
(disintegration energy 4-2mev.). The partial lifetime for «-decay is obviously 
very long, but it may nevertheless be possible to detect «-decay competing with 
the 22-year B-decay if a careful search is made. The neutron-deficient lead 
isotopes lie in region F, and on the theory of §3 would be expected to have an 
excess of about — 3 Mev., though the data for the gold isotope suggests that further 
away from the closed shells it may be considerably higher. If the former were 
the case, Pb 180 would emit «-particles of about 5 mev., and Pb170 of 8Mev., 
and in view of the gold data, probably higher. It is therefore possible that 
w-emitters will be found among the very neutron-deficient lead isotopes. 

The present method of comparing observed with calculated disintegration 
energies can obviously also be applied to beta-decay. The evidence there is not 
so clear-cut; for the B-energies are not known with such certainty as «-energies, 
particularly where there are y-rays also ; where electron capture is taking place there 
is so far no means of determining the disintegration energy directly; and for even 
nuclei there is the alternation of even—even and odd-odd nuclei. A preliminary 
analysis along these lines has been made by Dr T. R. Kaiser, who found certain 
sharp discontinuities for the heavier elements. In particular there is a sharp 
increase in the B-excess as one passes through the neutron-number 82. (‘This is 
clearly shown by the decay sequence for A = 137, the B-energy of Csi37 (N= 32) 
being 1:-4mev., and of Xe137 (N=83) being much higher, namely 4 Mev.) 
What one expects on the shell model is a sharp increase in going through a neutron 
shell, and a similar sharp decrease in going through a proton shell. ‘This compari- 
son does not give any useful evidence at present for the heavy nuclei we have been 
discussing, for the line of stable nuclei goes simultaneously through a proton and 
a neutron shell at Pb 208, and the effects cancel out. ‘The only reported f-emitter 
which should show an anomaly is T1209, but it does not. This may well be due 
to the presence of a so far unreported y-ray. 

47-2 


700 W. E. Moffitt 


ACKNOWLEDGMENTS 
I would like to acknowledge my debt to Dr. A. N. May, who in 1945 first called 
my attention to the interesting results of the present type of analysis, and to the 
members of the Computing Division of the Canadian National Research Council’s 
Montreal Laboratory, without whose tabulation of formula (4) the present work 
would not have been undertaken. 


REFERENCES 


Bropa, E., and Featuer, N., 1947, Proc. Roy. Soc. A, 190, 20. 

Co.uig, C. H., and Parsons, R. W., 1950, Proc. Phys. Soc. A, 63 (in the press). 

Featuer, N., 1949, Nucleonics, 5, 22. 

FrenberG, E., and Hammack, K. C., 1949, Phys. Rev., 74, 1877. 

Guueckaur, E., 1948, Proc. Phys. Soc., 61, 25. 

Mayer, M. G., 1949, Phys, Rev.. 74, 1969. 

NorpuEm, L. W., 1949, Phys. Rev., 74, 1894. 

PERLMAN, I., GHrorsO, A., and SzaBorc, G.T., 1948, Phys. Rev., 75, 17303 1950, Lbid., 77, 26, 

RosENFELD, L., 1948, Nuclear Forces (Amsterdam : North Holland Publ. Co.). 

Srasorc, G. T., and PeRLMAN,! ., 1948, Rev. Mod. Phys., 20, 585. 

Tuompson, S. G., Guiorso, A., RASMUSSEN, J. O., and SEABORG, G. T., 1949, Phys. Rev., 
76, 1406. 


Excited Electronic Levels in Conjugated Molecules— 
IV: Symmetrical Cyanine Dyes* 
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ABSTRACT. ‘Topics in the electronic theory of the colour of symmetrical cyanine dyes 
have been discussed by means of resonance formalism. Our conclusions for the formamid- 
inium ion have been confirmed by use of the molecular orbital method. As a result of this 
analysis, certain modifications in the qualitative symbolism of colour theory have been put 
forward in the hope that the factors governing absorption may be more clearly understood. 


§1. INTRODUCTION 


HIS note is concerned with the theory of the colour of symmetrical cationic 

| dyes of the cyanine type. Maccoll (1947) has recently reviewed the present 
status of this theory and has distinguished, very pertinently, between two 
different approaches to this subject. The dye chemist, on the one hand, is 
primarily concerned with a qualitative interpretation of colour on the basis of 
resonance theory, which may enable him to understand the relationship between 
colour and chemical constitution. ‘The molecules with which he is dealing are 
large and complicated, so that he requires a symbolism of wide generality. On 
the other hand, the quantum theorist is still perfecting his methods for dealing 
with the simplest of molecules, and may find himself embarrassed, if only 


temporarily, by the increasing demands of the dye chemist for such an all- 
«Pp ber IIT j i ot No : : 
en Rea. in this series is concerned with naphthalene, by Miss J. Jacobs, Proc. Phys. 
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encompassing theory. Nevertheless a synthesis of these different approaches 
has been attempted (Brooker 1942, 1945, Herzfeld and Sklar 1942) and proved 
to be quite successful. 

It seems, however, that the ready acceptance with which Brooker’s views have 
been greeted is liable to propagate an undesirable misunderstanding. For, 
as we shall show, the qualitative resonance symbolism which has been evolved 
is not entirely sound and, if the dye chemist is to rely on the results of qualitative 
considerations, it is surely desirable that the concepts from which these arise 
should be fundamentally well established. We shall endeavour to throw a 
clearer light on some of the factors governing the colour of symmetrical cationic 
dyes. It is not our intention to make any new calculations, but rather to develop 
a more reliable picture of the absorption process. 


§2. THE FORMAMIDINIUM ION 


Let us first consider the current theory of the colour of the simple forma- 
midinium ion: 


eM J + 
R,N=CH—NR,, a; R,N—CH=NR,, 2’. 


This molecule has two ‘extreme’ structures, a and a’, which do not ‘interact’ 
directly in the quantum-mechanical sense—that is, the matrix element 


Hg Wert eds 


vanishes in the usual approximation. Both structures have the same energy, 
E,,=E,,, so that the system is said to be degenerate. Another structure of higher 
energy E,, namely 


+ 
RACH Neen 


may be drawn, and this structure interacts with both a and a’. Brooker (1945) 
now envisages a resonance splitting of the degenerate energy levels of the extreme 
structures a, a’, in which the intervention of this ‘ excited’ structure is treated 
as a perturbation. Graphically such splitting is represented by Figure 1. ‘The 


Energy 


Figure 1. Erroneous splitting of levels of the formamidinium ion. 


resultant states e, g—of which g has the lower energy and is therefore identified 
with the ground state—are together considered as responsible for the first 
absorption maximum of the formamidinium ion. (Since (E,-E.)<(E, —E,), 
it is also deduced that the structure r is more important 1n the excited state than 
it is in the ground state.) The analogy between such symmetrical ‘splitting’ 
of the extreme energy levels, brought about by the contributions of r to the 
resonance hydrids, and the splitting of doubly degenerate 1H states of diatomic 
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molecules in the presence of magnetic fields, seems inescapable—nevertheless 
it is quite incorrect. 

Structures such as a, a’ do not, of course, represent actual states any more than 
the excited structure r. Spectroscopic states, represented in atomic orbital 
approximation by linear combinations of structures, may contain more of a, a’ 
thanofr. And, in perturbation theory, it is legitimate to use simplified functions, 
such as the structures a, a’, which may be considered as ‘splitting’ under the 
influence of the perturbing Hamiltonian. But sucha picture is apt to be misleading, 
especially when this ‘splitting’ is attributed to the intervention of the excited 
structures, such as r, and not to the Hamiltonian. 

In the first place neither of the structures a, a’ reflect the two-fold symmetry 
(Cy,) of the molecule, and are therefore unsatisfactory zeroth-order functions, 
mathematically. We use o to denote the plane of the molecule; the plane of 
symmetry which contains the central carbon atom and is perpendicular to « we 
shall denote by v. It is better to construct linear combinations 


1 1 
U 2 es 
Ee Sy 


which are symmetrical (A,) and antisymmetrical (B,) respectively with respect 
to reflections in v. This process of constructing functions which fully participate 
in the molecular symmetry from functions which do not do this is called group- 
theoretical ‘reduction’; we shall use the notation a, a—>a*, a to denote it. 
Since a, a’ do not interact with each other, the energies associated with at, a~ 
are the same as those of a,.a’; thatis, E,=H, =H; =E,. We see that 1 remams 
invariant to reflection in v and is therefore of the + species (A,). It thus interacts 
only with functions which are also of the + species. ‘This interaction of r+ and at 
will give rise to a state g+ whose energy Ej is less than E =E,, and also to another 


state ft, for which E; > FE. The energy diagram representing this analysis is 
given in Figure 2 (a). 


wr 


Aigo V2 es {V6 yh 


eee Ps(eta") 


mone | a oe OSC ae 


Batt apn Seka en elie 
gt(atrtb*c*) 
(a) 1st approximation. (6) 2nd approximation. 
Figure 2. Correct splitting of levels of the formamidinium ion. 


We see that, if other structures are neglected, the first excited state e must 
be identified with a~, and has the same energy as that of the extreme structures ; 
e~ therefore contains no contribution from r, and the formal charge on the central 
carbon atom Q, vanishes. The ground state, on the other hand, is only more 
stable than e~ by virtue of its component from r; Qe for the ground state is 
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therefore positive. It is easily shown that the intensity of the transition g+ <>e- 
is greater than that of g+<~>f*, though the latter is not forbidden. 

The contribution of the other possible structures may also be analysed. 
‘These are enumerated below: 


~ = + ++ = - 
P N==CH__NR. bh} REN=-CH_—-NR:..cs R,N—CH—NR,, -c’; 
See 

b is essentially of the + species, andc,c’+c+,c-. Onenergetic grounds bt, c+ will 
be less important than r+ in the ground state, but will serve to diminish Q, there. 
Similarly the contribution of c~ to e~ will be small, but will make Q,<0 and 
bring E- below Ey. ‘The more complete energy diagram now follows Figure 2 (6). 

In this simple case, therefore, we can arrive at a complete description of the 
qualitative factors governing absorption. ‘The ground state has + symmetry A,; 
the first excited state has — symmetry B,. The dipole strength of transitions 
between these levels is relatively high. It is only because of structures such as 
r+ and bt, which belong essentially to the + species and have no analogues in 
the — species, that the ground state lies below the first excited state. Structures 
giving rise to both + and — components, such as a, a’ and c, c’, could not alone 
be responsible for this separation. Excited structures are therefore more 
important in the ground state than in the first excited state. Both g* and e~ 
states are associated with a positive resonance energy—that is, E> and Ef are 


5 


both less than E, or E,. Qo is positive in g+ and negative in e~ with 
lOc(g*)|>1Ocle)]. 

Next, let us see how the replacement of the central —CH= group by a 
nitrogen atom —N= will affect the absorption. To a good approximation 
we need only consider the structures analogous to a, a’ andr. ‘The problem is 
exactly the same as was the case formerly. If, therefore, we assume that the 
analogue of the matrix element 


H>= | WY HY, dv 


is not very different from that encountered earlier, it is easily shown that 
(E; —E{) is smaller here than for the formamidinium ion. For, ceteris paribus, 
the greater the energy separation (EZ; —E;)), the smaller will be the contribution 
ofr+tog+. Thus the replacement of -CH= by —N= results in a bathochromic 
shift of the first absorption maximum (to longer wavelengths). Further, the 
formal positive charge Qx(g*) on the central nitrogen atom in the ground state 
will be less than Q,(gt). These predictions are confirmed by a molecular 
orbital treatment (see below). 


-§3. GENERALIZATION 
A similar picture may be obtained of the long wave absorption of higher 
vinologues of the simple system we have been considering. ‘These vinologues, 
namely 


RONE(CH ee NR cd, 2,3), 


where the R groups on any N atom may participate, together with N, in the same 
cyclic system, contain the symmetrical cyanine-type dyes. If we consider only 
nearest-neighbour interactions, the conjugation problem may again be associated 
with a plane of symmetry v, which contains the central carbon atom and is 
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perpendicular to the molecular plane o, even when » is not a plane of symmetry 
for the actual nuclear skeleton. (‘To this approximation, the actual shape of the 
conjugating system will affect only the intensity, and not the frequency, of the 
absorption.) Accordingly we may once more speak of states as belonging to the 
+ or — species. From the observed regularities which are associated with the 
ascent of such a vinologous series (e.g. Maccoll 1947) we may infer that, whatever 
the value of m may be, the ground and first excited states are of the + and— 
species respectively, just as for the formamidinium ion. 

Let us confine ourselves to some fixed value m, and retain the symbols a, a’ 
for the ‘extreme’ structures, which have a positively charged nitrogen atom 
and contain no formal bonds. Thus a’ is obtained by reflecting ain v. Each of 
the possible m=m(n) structures may be assigned to one of the three exclusive 
groups A, A’ and B: to each structure 7 of A there corresponds a structure 2’ 
of A’, obtained by reflectingzin v. Further, if7’ is any member of A’, the matrix 
element H,, =[V,4V’,, dv vanishes; also, a belongs to A. Structures belonging 
neither to A nor to A’ belong to B and are called ‘bridge’ structures. We may 
obtain functions of the + and — species from the structures of A, A’ by means 
of the various processes 7, 7’>i+,i-. In general, structures of B will also yield 
functions both of the + and — species; only when 2=0 or 1 do we obtain no 
— components from B; the number of linearly independent + functions which 
may be constructed from the members of B exceeds the number of functions 
of the — species for all values of n. 

Now, suppose we neglect the contribution of all the (excited) bridge structures 
B. Then it is easily shown that the energies of the lowest ‘states’ of both + and 
— types are equal to E,(n), say. If E,(m) denotes the energy of an extreme 
structure, it may further be shown that, as m increases, so does {E,(n) —E,(n)}. 
(For, of the &,, ., structures of A(n +1), we may select a set of k,, structures whose 
xk,,(k, +1) different energy-matrix elements are identical with those of the set 
of k,, structures which comprise A(m) ; and this selection does not exhaust A(n + 1).} 
To this degree of approximation, in which the contribution of the bridge structures 
is neglected, the ground and first excited states have the same energy and, as we 
ascend a vinologous series, so the resonance energies associated with these 
‘states’ increases. 

Clearly, therefore, no adequate description of the long wave absorption of 
these molecules can be obtained without the consideration of the bridge 
structures. And it may be shown that these will stabilize the ground state gt 
to a greater extent than the first excited state e~, although both Ex and E; will 
lie below E,. Once more we find that it is the greater importance of the excited 
structures to g* than to e- which is responsible for the separation of the two 
levels. 

Let A(n), B(n) contain k,,, 1, structures respectively. Then it is easily shown, 
by simple enumeration, that the ratio 1,/k, decreases fairly rapidly with 
increasing values of n. (For example, if we consider only those structures with 
no negatively charged carbon atoms and at most one formal bond, then 
L,/k,=1/2(n+1). With a more representative selection of structures, the ratio 
L,)R, decreases more rapidly than this.) ‘This suggests, straight away, that as 
” increases so the extra stabilization of g*(m) with respect to e-(n) decreases. 
That is, the sequence {Ei (n) —E-(n)}4 diverges and the vinylene shift of the 
absorption maximum is ‘non-convergent’ (Maccoll 1947). Conclusions derived 
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by means of such purely qualitative arguments should, whenever possible, be 
supplemented by more rigorous proofs—such as that put forward by Herzfeld 
and Sklar (1942). But they illustrate the way in which our use of resonance 
symbolism may be applied. 

An essential feature of Brooker’s theory (1945) was that the resonance energy 
{E,(n) —E{(n)} should decrease with increasing values of n. In general, however, 
the calculations of Herzfeld and Sklar (1942) belied this conclusion, and the 
experimental evidence adduced in its favour—that the acidity required to 
decolourize a series of symmetrical cyanines decreases with n—is capable of a 
different interpretation. At best it only shows that the difference in the resonance 
energies of two systems, both of which are conjugated, decreases with increasing 
values of n. (This behaviour is, of course, not surprising because we should 
expect the differences between such odd- and even-numbered chains to decrease 
with n.) ‘The more detailed analysis, which we have adumbrated above, does 
not require a decreasing resonance energy. On the contrary, since {£,(”) —E),(n)} 
increases with m, we may infer that, at least ultimately, {£,(m) —Ej(n)} also 
increases with , For, in the limiting case of a long chain, the non-convergence 
of eae gives Ey (n)—> E.(n)> E, (7). | 

We may draw another corollary from this form of treatment. Suppose 
we have two molecules of the cyanine type, of the same chain length and differing 
only in, say, a particular substituent. If the importance of this substituent may 
be localized, as affecting only the energy of a particular bridge structure, then 
its effect on colour may be predicted. For if the substituent stabilizes the bridge 
structure, it is easily shown that such substitution is accompanied by a hypso- 
chromic shift of the absorption maximum, and conversely. This theorem 
is a particular case of a more general colour rule used by Brooker and recently 
extended in application by Knott (1950). 


§4. THE MOLECULAR ORBITAL METHOD 

In order to confirm our qualitative conclusions for the formamidinium ion, 
we shall apply the molecular orbital method to this system. We must then 
assign the four 7-electrons to molecular orbitals s which have modes ino. ‘The 
ys are expressed as linear combinations of atomic orbitals and are determined using 
the ‘self-consistent’ method recently introduced by the author (Moffitt 1949). 

The electronic structure of the ground state of the formamidinium ion may 
be described by means of the configuration 

Pye], 

where the superscripts + and — refer to the transformation properties of the 
individual orbitals (one-electron functions) with respect to v. ¢* is strongly 
bonding and % is non-bonding; [i]? therefore represents two non-bonding 
electrons, one on each nitrogen atom. ;/* has the form 


1 
pr=yedge + V2 yx(oy +n’), 
with ve =SinlU, Ve = COS 0, 
and @ so chosen as to minimize the zeroth-order energy of the configuration and 


therefore satisfying 
(5 — EX) +2/2Bon cot 20 =0; 
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ét and €¢ describe the nature of the Hartree field near the C and N atoms (see 
Moffitt 1949). In terms of valence state ionization potentials J and electron 
affinities F, this becomes 


(31, — 21, —2E, + Ex) — cos 20{2(Ig — Ec) +(Iy — Ex)} +8V/2Bexy cot 20 =0. 
The z bond order and the resultant charge distribution are 


Pix = 73 sin 26, Ox = 4(1 —cos 28), Og =cos 20. 

Numerical values of the J, E are given by Mulliken (1934): J,=11-2, 
E,=0-7; Igy =12-:2, Ey~1-5 ev. Bey, the carbon-nitrogen resonance integral, 
has-a negative value in the neighbourhood of —3 ev., but its precise value is 
unfortunately not known. We have therefore solved for 6, and hence determined 
Pix, Ox and Qgy for Boy = —2, --3 and —4 ev. respectively. Our results for 
the ground state are to be found in the first three columns of Table 1. 


Table 1 
Bond Orders and Charge Distribution in the Formamidinium Ion 
Ground state First excited state 
Bon PEN On Qc PON On Qc 
—2 ev. 0-683 0-370 0-260 0-073 0-506 —0-011 
—3 ev. 0-690 0-392 0-216 0-108 0-512 —0-024 
—4 ev. 0-695 0-408 0-184 0-141 0-521 —0-041 
Table 2 
The Effect of Replacing —CH= by —N= in the Formamidinium Ion 
Ground state First excited state 
Bun PNN On OX PNN On OX 
—2 ev. 0-700 0:428 0-144 0-067 0-505 —0-010 
—3ev. 0-702 0-439 0-121 0-100 0-510 —0-020 
—4eyv. 0-703 0-448 0-104 0-132 0-518 —0-036 


The first excited (singlet) state e~ arises from the configuration 
Pol: 
s+ is again bonding, # is non-bonding and y+ is antibonding. Setting yZ=sin 80, 
as before, we find that 
2(Iy — Ec) + cos 26{2( 1 — Eg) +(Ly — Ex)} +44/2Boy cot 20 =0. 
The z bond order and charge distribution, 


Nae 
Dex = TET 20, n=t(3 —cos 20), Qo = — 4 (1 —cos 26), 


are given in the last three columns of Table 1, for the same three values of Bon: 
The qualitative conclusions as to charge distribution, derived above by means 
of resonance theory, are completely substantiated by the molecular orbital 
method. By the method of antisymmetrized molecular orbitals, using such 


LCAO forms, it is also possible to calculate the energy separation (E, —E*), 
but this is not our immediate intention. 3 
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The effect of replacing the central —CH= group by —N= is easily seen. 
We simply replace (I,, Ey, Bey) by U%, EX, Byx) in the above formulae, where 
IX, EX and Jy, Ey refer to different valence states of the nitrogen atom. Mulliken 
gives the values /* =13-8, EX =1-0 ev. If we now take the same set of values, 
namely —2, —3 and -—4 ev. for Byy as for Boy, we obtain the results of 
Table 2; QO refers to the charge on the central nitrogen atom. ‘The differences 


between Tables 1 and 2 are in complete agreement with the results of our 
previous qualitative considerations. 
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ABSTRACT. Experimental arrangements are described by which clearly resolved 
spectra have been obtained of the luminescence, at the temperature of liquid air, of seven 
different phosphors of simple crystal structure (cubic) and of mixed crystals. ‘The main 
peaks in each case may be arranged in a band series. The vibration frequency in the lower 
state is nearly the same for the same host lattice with different activators. The energy of 
the electronic transition is determined both by the activator and the host crystal and there 
is a progression of decreasing energy in the series of simple crystals MgO, CaO, SrO, CaS, 
and in the order CaO, CaO-SrO with mixed crystals. "The most probable transitions are 
from the lowest level of the upper state. There is evidence that self-excitation occurs and 
for the view that the luminescent centres are large. 


SIS ENERO DiC VON 

MONG the many inorganic solids showing luminescence under irradiation 
by light, x-rays, cathode rays or radioactive radiations, only a few give 
emission spectra which can be resolved, even at low temperatures, into 
lines or narrow bands. ‘The rare earths, both as pure solids and in solution, 
sometimes give line-like spectra which seem to be attributable to transitions in 
a shielded shell of the metallic ion. The uranyl] salts as pure solids or in solution 
give well resolved spectra which appear to arise from transitions in a shielded 
shell of the uranyl ion. ‘The narrow bands obtained with alkali halides have 
been interpreted by von Hippel (1936) as arising from the transference of an 
electron from a halogen to an alkali ion. Deutschbein (1932) explains the sharply 
resolved luminescent spectra sometimes obtained with Cr activated solids as 
arising from transitions in the Cr+++ ion modified by vibrations of the host 
lattice. For the most part, however, the luminescent spectra of impurity activated 
solids are broad and featureless even at liquid air temperature and lower, and it 
is not possible to evaluate with any accuracy the energy of the electronic transitions 

or to interpret, with any confidence, the way in which they arise. 

Randall (1938) has pointed out that the similarity of the luminescent spectra 
of some solids activated with Mn to the spectra of pure manganous salts 
suggests that transitions in a shielded shell of the Mn*+ ion are responsible in 
such cases, and surmised that the fluorescence of pure tungstates and molybdates 
originates in the WO, and MO, ions. Kroger (1948) has offered evidence support- 
ing this view. One of us (Ewles 1938) reported previously that parts of the 
luminescent spectra of CaO—‘Specpure’, of CaO activated with Bi and of SrO 
could be resolved at the temperature of liquid air into narrow bands susceptible 
of interpretation as electron transitions modified by vibrations of the host lattice. 
The present contribution reports that this work has been extended to other cases 
which may be regarded as typical of a large class of impurity activated solids 


and suggests that the results may be of some use in the interpretation of the 
juminescence of this class. 
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§2. EXPERIMENTAL ARRANGEMENTS 

The procedures for the preparation and heat treatment of the substances 
studied are described later, along with the analyses of the spectra. ‘The source 
of excitation was a Mazda 125 watt Mercra Mercury vapour arc with the outer 
glass envelope removed. A clear silica flask of about 10 cm. diameter filled 
with a solution in distilled water of 145 gm. NiSO,.7H,O and 41°5 gm. 
CoSO,.7H,0 per litre (Bowen 1935) served to focus an image of the source on 
the specimen and to remove all light between 33414. and 5771a. ‘The strong 
yellow lines at 5771a. and 5779a. may be removed by the addition of a piece of 
B.T.H. ‘black glass’ or Corning 9863. This addition, while reducing the 
intensity considerably in the shorter ultra-violet, 1s very useful for preliminary 
visual examination for visible luminescence, providing a source virtually free 
from all visible light, extending well into the far ultra-violet, and so strong that 
under its illumination very few substances are entirely free from luminescence. 

An image of the irradiated specimen is focused on the slit of the spectrograph 
(Hilger E 383—quartz prism with Littrow mounting) with a glass lens which 
serves also to reduce scattered far ultra-violet light. A simple method of mount- 
ing the specimens for obtaining the luminescent spectra at low temperatures is 


) >) 
Cu Cn 
E A Specimen. 
B Rhodium plated copper plate. 
D Spring clip. 


Q Thin quartz plate. 

Cu, Cn Copper, constantan wires. 

E Dewar flask of transparent silica. 
WwW Heating coil. 


Figure 1. 


shown in Figure 1. A very thin layer of the specimen about 1 cm? in area is 
deposited at A on a heavily rhodium plated copper plate B 4in. long, 14 in. broad 
and 1/16 in. thick. The specimen is covered by a thin quartz plate Q held firmly 
against the copper plate by a spring clip D. This arrangement is suspended by 
thin copper and constantan wires Cu and Cn, which also serve as a thermocouple, 
in a Dewar flask of transparent silica E containing the cooling liquid, so that 
the copper plate just touches the bottom of the Dewar and the specimen is above 
the level of the surface of the coolant. It is found with, for instance, liquid air 
as the coolant that the temperature as indicated by the copper—constantan couple 
remains at that of liquid air as long as any of this remains in the Dewar. ‘To 
prevent condensation on the quartz plate covering the specimen it is necessary 
and effective first to fill the Dewar to the required depth with the coolant and to 
delay the immersion of the copper plate until the whole of the inside of the Dewar 
is well cooled—a matter of a few minutes only. Condensation on the outside 
of the Dewar in the neighbourhood of the specimen is prevented by an electrically 
heated wire W wound round the flask just below the level of the specimen. 
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§3. PREPARATION OF SPECIMENS AND ANALYSIS OF LOW 
TEMPERATURE LUMINESCENT SPECTRA 


The preparation of the phosphors described below is indicated in each case. 
A Cambridge recording microphotometer was used to analyse the spectrograms. 
Wave-numbers were obtained from interpolation from superimposed Hg lines 
using the Hartmann dispersion formula where the resolution was sufficient to 
make it worth while, and in other cases a wave-number-—distance calibration was 
obtained from a Hg spectrum. The probable error in the peaks is about 20 cm™?. 
Reduced photometer curves are shown together for comparison in Figure 2. 


Hg 3391 Hg 3650 ae geass 
| | 


Mg0-Bi 


Ca0-Bi 
(dotted line 20°c.) 


Sr0-Bi 


Ca0-Pb 


vp Sr0-Pb (?) 


CaS-Pb CaS-Bi 


Figure 2. Photometer curves of luminescent spectra (all at temperature of liquid air 
except dotted curve CaO-Bi). 


CaO-Pb and CaO-Bi 


The resolved ultra-violet and violet spectra of CaO-Pb and CaO-Bi have 
already been reported (Ewles 1938) and the results are given again here for 
convenience of reference and discussion. In the earlier report the CaO—Pb 
spectrum, since it was obtained with a Hilger Specpure specimen, was attributed 
to the pure CaO but, as recently reported by us elsewhere (Ewles and Curry 
1949), we have now reason to believe that Pb is the activator. 

The CaO used in the present work was made from the oxalate prepared by 
fractional precipitation by recrystallized analar ammonium oxalate solution from 
recrystallized calcium nitrate solution. After adding activator and flux, mixing 
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and drying, the mixture was heated for about one hour at 750°c. The analysis 
of the photometer curves shows that the peaks follow electron band formulae: 


CaO-Pb 
v=27526+ (v’+-$)264—(v” + 4)508 cm-} 
Obsd. Calcd. v, vu" Obsd. Calcd. GDh ey ol 
27404 27404 Oy 26385 26388 Oe 
27166 27160 hse tl 26142 26144 ee3 
26891 26896 Om 25886 25880 O38 
26647 26652 iS 822 
CaO-Bi 
v= 25815 + (v’ +4)298 —(v” + 4)493 cm=! 
Obsd. Caled. Cees Obsd. Caled. v, v0" 
26322 26314 Dr () 25017 25030 ih, ® 
26023 26016 i 24732 24732 O; 2 
25848 25821 Dey sl 24521 24537 EES 
25719 25718 On 24245 24239 O, 3 
25525 DSSS (Ls ail 23758 23746 0, 4 
25230 25225 Ona 


SrO-(Pb?) and SrO-Bi 


Strontium carbonate was fractionally precipitated by recrystallized analar 
ammonium carbonate solution from recrystallized strontium chloride solution. 
After washing and drying, one sample of the carbonate alone, and another after 
addition of Bi(NO;), solution, were heated in pure alumina crucibles for one 
hour at 1,100°c. to reduce to oxide. The sample containing no additional 
activator showed a weak resolved band in the near ultra-violet. The activator, 
if any, responsible for this is unknown. Although additions of Pb do not seem 
to intensify this band the recent findings with CaO—Pb referred to above suggest 
that here again some impurity such as Pb, which it is impossible to remove, may 
be responsible. 


SrO (?) 

y=26538-+ (v’ +4)230—(v’’ +4)386 cm 
Obsd. Galeds v, v0" Obsd. Caled. wee 
26460 26460 Onn0 25690 25688 @. 2 
26070 26074 Oneal 25510 25532 il 8 
25930 25918 ‘pe? 25300 25302 Oy. 2 

SrO-Bi 

v=24180+(v’ +4)190—(v” + $)390 em-* 
Obsd. Calcd. U, v0 Obsd. Calcd. v’, vu" 
4070 24080 @,. 0 23100 23100 (G8 
pee 23880 peat 22910 22910 O, 3 
23700 23690 OQ; 22730 22710 i, 4 
23480 23490 i, 2 22450 22520 0, 4 
23300 23300 Oe 
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MgO-Bi 
Analar MgO, after mixture with the activator and 10°, NaCl as flux, was 
heated for 5 hours at 700-800°c. This product gives a weak band in the ultra- 


violet with two well articulated peaks and a number of less defined peaks following 
the arrangement: 


y=28460-+(v’ +4)156—(v’” + $)346 em™ 


Obsd. Calcd. Vv, Uv’ Obsd. Calcd. v, Uv 
28350 28366 ORO 27670 27674 Ox 

28180 28176 theta 27480 27484 1 ars 

28020 28020 Os Al VI RVAY) 27328 Oo; 3 

27820 27830 ie 2 


CaS-—Bi and CaS—Pb 


Calcium sulphide was prepared by intimately mixing analar calcium sulphate 
with an amount of pure sucrose corresponding to the relation CaSO, +3 C. 
After addition of the activator and drying, the mixture was heated in an alumina 
crucible with a close fitting lid at 900° c. for about an hour. 


CaS-Bi 


This gives a fairly well articulated band in the violet whose peaks follow the 
electron band formula: 


y=23397 + (v’ +4)130—(v’”” ++4)250 cm 


Obsd. Caled. vv" Obsd. Calcd. U,V" 
DE Bi, 235011 On 0 22588 22587 ONS 
23220 RING Ibeeril 22372) 22351 0, 4 
23093 23087 Onn 22086 22087 (OSS 
22980 22967 jt? 21950 21967 ie XG} 
22820 22837 On 2 

CaS—Pb 


This gives a weak band in the ultra-violet showing two very broad ‘lines’ 
on a rising background and three fairly well articulated peaks. The analysis of 
the photometer curve gives a frequency difference for the lower state similar to 
that obtained with CaS-—Bi but there is not enough detail to make a reliable 
detailed analysis. ‘The peaks and frequency differences are: 


28140 \ 
AN 
280 


27860” > 2x250 em- 
oe 
27641 
N 
27483 > 244 
: Ye 


ie 
27397 
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Mixed Crystals: CaO—-SrO-Bi 


Mixed crystals of the oxides were attempted by co-precipitation of the car- 
bonates and, after addition of the activator and in some cases of up to 20°/, NaCl 
as flux, drying and heating for several hours at 1,100°c. The analysis of the 
photometer curves of the low temperature spectra show that, as the SrO content 
is increased, 

(a) there results a progressive loss of resolution results up to a Sr/Ca ratio 
of 1/4, where resolution is completely lost and not recovered until the SrO-Bi 
spectrum appears at about 95° SrO; 

(6) there is a slight progressive shift of the peaks to longer wavelengths. 


The measurements given below and the photometer curves in Figure 3 for 


Ca0-Sr0-Bi 
Ca0-Bi 
| | | = eee 
25500 25000 24500 24000 


[> (cm-') 


Figure 3. Comparison of resolved violet bands of CaO-(B1) and the mixed crystals CaO-SrO-(Bi). 


a Sr/Ca ratio of 1/5 show that there is also a slight change of vibrational frequency 
in this mixed crystal as compared with CaO. 


CaO-Bi CaO-SrO-(Bi) 
25719 25660 

488 480 
25231 25180 

499 480 
24732 24700 

487 480 
24245 24220 


§4. FACTORS INFLUENCING INTENSITY OF EMISSION AND 
RESOLUTION 

Samples prepared with the addition of 10°%, NaCl as flux gave more intense 
emission in all cases than those without flux, generally by a factor of some 3 or 4 
times. Resolution was also improved by a flux and, in the case of SrO and the 
mixed oxides, only obtained with the use of a flux. The spectra are, however, 
completely unaffected in electronic energy and vibrational frequency by flux— 
even up to 20% of this. The requirements regarding purity varied; in some 
cases, as with CaO, CaS and MgO-Bi analar materials were sufficiently pure to 
give well resolved spectra, while with SrO and the mixed oxides resolution was 
only obtained after careful further purification from analar materials. In general 
too, prolonged heating at a lower temperature than that normally required to 
initiate luminescent power improves both intensity and resolution. Grinding 
and rolling in an agate mortar decreases intensity. 
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§5. SURVEY AND DISCUSSION OF THE OBSERVATIONS 


The resolution of the luminescent spectra of a number of different solids of 
simple crystal structure activated by different impurities makes it possible to 
draw conclusions which may be of general application to the interpretation of 
the luminescence of solids. 

(a) The spacing of the vibrational levels in the lower state involved in emission 
is nearly the same for the same host lattice activated with different impurities. 


- 


Host Activator 

Bi Pb 
CaO 493 cm! 508 cm7} 
SrO 390 386 
Cas 250 244 
MgO 346 — 


This seems to indicate that the vibrational frequency of the luminescent 
centre is to be identified with a vibration of the host lattice. 

(b) The energy of the electronic transition is determined both by the activator 
and by the host lattice, e.g. 


Electronic energy in Electronic ener in 
peek wave-numbers (cm~?) SURES wave-numbers Cnet 
CaO-Pb DYESS CaS-Bi 23397, 
CaO-Bi 25816 CaS—Pb about 28000 
SrO-Pb ? 26538 MgO-Bi 28460 
SrO-Bi 24180 ; 


(c) There is a progression of decreasing electronic energy of transition in the 
order MgO, CaO, SrO, CaS with these simple crystals and in the order CaO, 
CaO-SrO with mixed crystals. 

(d) If the proposed analyses of the resolved bands be accepted, the strongest 
components all represent transitions from the lowest level of the upper state. 

(e) The relative weakening of some of the early members of the series (seen 
by comparison of the numerical analyses with the photometer curves) is due, it 
is suggested, to self-excitation by the overlapping of the long-wave tail of the 
excitation region with the short-wave tail of the emission region. Excitation 
curves obtained by Ewles and Curry (1949) indicate that such overlapping does 
indeed occur in these cases. 

(f) Both the levels involved in a transition giving rise to emission of a member 
of a band series must be very narrow. In CaO-Bi for example the half-width of 
a member is about 30 cm corresponding to a width of level of about 0-008 ev. 

(g) It is suggested that the observations are in accord with the view that the 
luminescent centre consists of an impurity (a term which may include a lattice 
defect) associated with a large but not too large number of lattice ions. This 
view, originally put forward in some sort by Lenard in 1904 and developed 
quantitatively by one of us in 1930 (Ewles 1930) now appears to be receiving new 
support from various sources (Ewles and Curry 1949, Gibson 1948, Garlick 1948). 
It has been shown (Ewles and Curry 1949) that in doubly activated calcium oxide 
Pb centres interfere with Bi centres when the concentration of the lead is iecwecn 
10-4 and 10-5. If the ‘size of the centre’ is m lattice spacings each way from the 
impurity then interference would occur when the chance of one lead atom occurring 
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in a volume occupied by (27)? calcium atoms is appreciable. If the chance for a 
detectable quenching is put at 10% then Pb/Ca=5~x 10-5=1/10(2n)* giving 
n=6 lattice spacings each way for the effective size of the centre in this case. 
It does not appear that there is any evidence from the observations reported here 
as to whether the association of the Bi or Pb with the crystal is by replacement of 
a lattice ion or by interstitial entry into the lattice. 

(A) Seitz (1938) and Williams (1948) have estimated the peak absorption and 
emission energies in the luminescent spectra of some impurity activated solids 
in terms of transitions in the activator atom in association with the host lattice. 
Attempts are being made to apply such methods to the spectra reported here. 
The precise experimental evaluation of electronic and vibrational energies, made 
possible by the resolution and analyses given in this paper, should afford a more 
definite test of such methods than the broad bands of ZnS(Zn), ZnS(Cu), 
ZnO(Zn) and KCl(T1) to which they have been applied so far. 
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ABSTRACT. Kirkwood’s formula for dielectric constant, €, gives it in terms of a corre- 
lation coefficient which depends on the probability that any two dipolar moleeulés have a 
certain configuration. - 

The main object of this paper is to calculate the configurational probability of the dipole 
orientations; this is dealt with in the first part of the paper by considering two dipoles and 
replacing the rest by a homogeneous dielectric. The dipolar molecules are regarded as 
point dipoles and the effect of other forces on the orientation of the molecules 1s neglected. 

In the second part of the paper the calculation of the correlation coefficient appearing in 
Kirkwood’s dielectric constant formula is carried out. This correlation coefficient is 
evaluated for water, and the value of e€ then obtained from Kirkwood’s formula is in fair 
agreement with the observed value. 


INTRODUCTION 


N early relations of the dielectric constant, the interaction between dipolar 

if molecules has not been taken into account. Kirkwood (1939) and then 
Frohlich (1948) constructed formulae which depend on a knowledge of the 

configurational probability of dipolar molecules. 

Kirkwood, in the application of his formula to water, considered that for any 
two first neighbours (owing to hydrogen bonds between them) the two dipole 
directions make a constant angle with the line joining them, and can rotate freely 
about that line. The value of the dielectric constant of water obtained on that 
assumption 1s fairly good. 

These formulae are a useful means of checking the accuracy of the configura- 
tional probability calculated by other methods. 

A knowledge of the configurational probability of the dipoles can also be used 
in calculating an additional average potential energy and an additional specific 
heat due to dipole-dipole coupling, the effect of dipole rotation on the x-ray 
scattering, the hydration of salts and many other things. 

In the first part of the paper a method for the calculation of the configurational 
probability of dipole orientations is carried out on the assumption that the 
molecules are point dipoies. The two dipoles concerned are considered as in 
a cavity (or two cavities), the mutual potential energy is formed and the 
probability for such a configuration calculated by use of Boltzmann’s law. 

In the second part of the paper the calculation of the correlation coefficient 
is given in terms of double integrals which can only be evaluated by numerical 
methods. ‘The numerical evaluation of this correlation coefficient for water 
leads to a value of the dielectric constant in fair agreement with the observed value. 
It may be noted, in the application to water, that the effect of hydrogen bonds 


has not been taken into account, since we have used a model in which all molecules 
were substituted by point dipoles. 
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THE CONFIGURATIONAL PROBABILITY 


The rigorous method of evaluating the probability p,/(pjp,) that two 
dipoles 7, 7 have simultaneously a certain configuration is to write down the 
probability for the whole set of dipoles in the form 


pZ =exp(— W/kT) 
where 


Z= % exp(—W/kT), see eee (1:1) 


Pj---Py 


and then sum p over all orientations of the dipoles except those for the two 
dipoles 7,7 in question. Since it is impracticable to evaluate such a sum, we have 
to use some approximate method to find the required probability. 

We replace the rest of the molecules macroscopically by a homogeneous 
dielectric. The form of the internal boundary will be assumed to be simply 
two separate spheres with the two molecules concerned at their centres; the 
volume of each sphere is 1/N =47a?/3 where N is the number of molecules per 
unit volume. If the two molecules are first neighbours, the two spheres will 
overlap and thus form a single cavity of volume less than 2/N; we might therefore 
consider the radius of each of the overlapping spheres to be xc where c is the 
distance between first neighbours. In this case the volume of the cavity would 
be wc3(2x + 1)?(4x—1)/12. Equating this to 87a3/3 we obtain 


x 1432@/8x 143203)’ 


which is a cubic equation to determine 1/x and hence the radius xc of each of 
the overlapping spheres. 

The evaluation of the field inside such a cavity is, however, found to be very 
difficult and therefore a cavity in the form of an ellipsoid of revolution will 
actually be considered. If the distance between the foci is 2b and the major axis 
is 2b cosht, then equating the volume of the ellipsoid to 87a*/3 one obtains 


3 
(cosh t) — (cosh t) —2 ei ZO ee (re es (1.3) 

This equation, although it has fixed the volume of the cavity, has not fixed 
the form of the ellipsoid of revolution. This will be chosen to coincide as far 
as possible with the cavity of the two overlapping spheres whose dimensions 
are given by (1.2). 

In order to calculate the probability po, we shall calculate the mutual potential 
energy between the two dipoles in a given configuration, taking into account 
the image field, or the average effect of the rest of the dipoles. 

If we write the total field acting on one of the dipoles z in the form 


=a. F = (Ay— 1p +(By— Ag)(H-K)k + Au” + (B-Ab), 
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then when the two molecules are first neighbours, the constants Ap, Bo, A,B 
will be shown in the first Appendix to be 


: ) 
eee Maley) ON Com) meas 
Ay—1= b? os (m+) {Ps (=) ae t) 1—-X,” 


i 2h . a(COSIUP eet 
Pseaee Ts (7 +3) )4P, a) } P, (cosht) 1—X 


n=0 
coms eee a \n c Q,(cosht) 1 | 
Dee aie coe (5) Pi(cosht) 1—X®” 
Wena at 1) c O,(cosht) 1 
Bee fekcigie 6 ee (n-43)4 Pas) } Se P,(cosht) 1—X,,” | 
Laren (als) 
where 


x, — Qn(cosh t)P,(cosht t) —Q/(cosh t) P,,(cosh f) 
Cay P,,(cosht)Q/(cosht) 

sinh?¢ Q/(cosht)P%(cosh t) — Q/,(cosh t) P;(cosh ¢) 
Pi(cosht)’ — Q/(cosht) cosht + Q;(cosh ¢) sinh? ¢ 


(e io IDA ay 


2b is the distance between the foci, 2b cosh t is the major axis and c is the distance 
between first neighbours. 

It may be noted that for large «, X,,, X¢?—>0 and so all Ay, By, A, B will be 
independent of «. 

The field acting on one of the two molecules given by equation (1.4) in vector 
form can be re-written in the tensor form 


—aFO= x (A— 1)o a) ee Giz) 


gush 


If the axes considered are orthogonal then A will be given by 
Nig = AF =A Soe + (By — Ao )Reky =y™, say, 
Ass! = ree! = Ady +(B—A)kRy =Yee, SAY, 


where k is a unit vector along the line joining the two molecules. Omitting 


the suffixes s, s’ then 
(0) 
ae eS 
ie @ i) ck eee (1.8) 


where 
Y= Ay + (Bo — Ao)(RR) 
Siete (1.9) 
y=A +(B—A)(ER). 
Since 
EO? ep hig EO oe oie eer (1.10) 
then, eliminating F(” between (1.7), (1.10) we obtain 
T=) AG ee 4 Sk Gah eae CESi1) 


js 
where A“! is the reciprocal of the tensor A. 
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The mutual potential energy W is given by 
ow a leah came 
eG sD wars ict) Damen ime (1.12) 


By use of (1.7), (1.11) and (1.12) we obtain 


W= —ZD(A)2pOpM +const. —... se . (£.13) 


20 ij ss’ 


It is proved in the second Appendix that if A is given by (1.8) and f(A) is any 
function of (A) then 


OSA Sayer ya yey) 
eo Aaa ny ea fy). 


In particular taking f(A) = —2“4, and omitting suffixes s, s’ 


== yl) = Fy 
yor — 2 AY Y)> Say, 
= QP = stg = FY), say. eee (1.14) 


The expression (1.13) can then be written, using (1.14), in the form 
1 * 
W = 5S [Fay roel 909? + DEP + = ELF Mae PPP oe (1.15) 


where the tensors y®, y are given by (1.9). 
It is proved in the third Appendix that if F(y) is a tensor function of the 
two tensors y and y given by (1.9), then 


F(y) = F(A, A) +[F(BoB) — F(AcA)|(RR). wv es (1.16) 


A similar equation holds for Fy(yy). Substituting in (1.15) the values of 
F(yy) and F(yy) given by (1.16) and remembering that 


> pM PPR ky = (p® .k)(p®. ks) and = pip,” 8 = (pp), 


we obtain the final formula for W 


We = [Fo(BoB) — Fo(AoA)]I(p® - k)? + (pe. k)"] 


3; : {F(A,A)(p® . p®) +[F(BoB) — F(Ao4)I(P® - k)(p® - k)}. 


The functions F, and F are given by (1.14) on replacing y, » by the 
constants Ay, A or Bo, B. 

The probability for a given configuration can then be easily found by 
Boltzmann’s law, using for W the expression given by (1.17). 
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APPEN Dial 


The internal boundary of the dielectric is an ellipsoid of revolution of major 
axis (along the axis of revolution) 2bcosht and the distance between the foci 
is 2b. A dipole w lies at distance c/2<b from the centre of the ellipsoid and on 
the line joining the two foci. It is required to find the field at the position of 
the dipole and at a point on the line joining the two foci at distance c from the 
dipole y. (on the other side of the centre). 

Using orthogonal coordinates (t, 0, f) given by 


x, =bcoshtcos6 %,=bsinhtsin 6 cos ¢ x%3=bsinhtsin@sind 
where (x,, X2, x3) are orthogonal axes, the first along the axis of revolution, the 


second in the plane of the dipole w and the first axis, the general potential V 
satisfying V?V =0 can be written in terms of (t, 0, ¢) 


V(t0d)= % & [A”Pi(cosh t) + B”O"(cosh t)] T”(cos 6) cos md 


n=0 m=0 
nee (1.18) 
where TI) = (1 — a2 inl 2/ car ym —— : 
mM m aie, x m a n\& 
P®(x)=(x¢2.—1) Ts OM (x) = (02 — 1)m2 fen NE 


The coordinates of the dipole are (0, 4g, 0) and those of the other point at which 
we want to find the field are (0, 7 —@, 0) where cos 6,=c/2b. If R is the distance 
between (t, 6, 4), (t’, 0’, 4’) then, 


1 fo) 
R75 x (2n +1) [ 2, (cosh t)P,,(cosh t’) P,,(cos 8) P,,(cos 6’) 
n=0 


: m (n—m)! m ™ , , 
42 (=) ea On'(cosh t) P™(cosh t’) T"(cos 0) T(cos 6’) 


xcosm(p—¢') |. Saeenae (La9) 


Taking a charge —e at t=0, 6=6, and a charge +e at t=t', 0=0', then if 
the dipole makes an angle « with the direction of the axis ioe 


pe ; : E ; 
he cos «=cosh ¢’ cos 6’ — cos 05: - sina =sinh tf’ sin@’, 
e 


The potential of the dipole p at (0, %, 0) is then, by use of equations (2:19) 
and (1.20), given by 


V= = > (2n +1)P/(cos 5) {cos «Q,,(cosh t) P,,(cos @) 


n=0 


sin « 
+ Lael Q'(cosh t) T\(cos 8) cos s} Sik, Se ae (1.21) 
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To find the image field inside the cavity, let the potentials outside and inside 
the cavity be, following (1.18), 


Vo= X& & CROM(cosht)Ti"(cos 0) cos mp 


n=0 m=0 


Vi= 2 [ A” P"(cosh t) + B”O'"(cosh t)| T7"(cos 6) cos m¢. 


n=0 m=9 


The second part in the expression for V, is the potential of the dipole p 
and is given by (1.21); the first part is the additional or image potential inside 
the cavity. Applying the boundary conditions V)=V,, «(@V)/dt) =0V,/dt at the 
boundary t=¢, one obtains 


Gn) == O7n'(cosh t,) 1 (m) bd 
Ae == Puchi) Take (1.23) 
(m) ae — pm te 
nt t i if 
where (e—-1) x= 3 Gi BR ee 2 (1.24) 
[Pe ee 


Using the expression for BY” from (1.21) we obtain 


Q,(cosht,) 1 cosa 


Oye D 
hat aap (cosh s,) Lo XO Prt h Ecos) abe 
ay _ Qrlcosh t)) 1 sina (2n+1) ,, ee 
ee cach) ISK. Be aad) 
where 
(-xo= (2g Ses) 
sae bars RS Wrath, ited totes ee See (1.26) 
ane eae oem) 
” ~ Pi(cosht,) \Q;,cosht+Qysinh?z/;-, 3 
The field components are given by 
—1] OV —l OV 
b,/(cosh® t—cos?0) dt’ by/(cosh® t—cos* 6) ° 00” 
and at 6=6 , t=0 the image field inside is given by 
fhe a8 x AMP*(cos 69) —>— Me wel) cos ¢ 
EG Te MM es lies calle) SR Ee (1.27) 
FO : x A P/(cos ), 
b n=0 
and at 0=7—6 , ¢=0 
oo j 
R= y (—)"ADP/(cos 9) nk "5 CES ae | 
De et he igi he See amc (1.28) 


Baily ace 
Ege a (—)"A P,(cos 69). | 
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Putting the fields at x, =c/2, x, = —c/2 as 


OES eviibieyfn (eye ah f= Ce “Diack os ae (1.29) 
respectively, then t 
Fi) = Cy sina cos d, F,= Cysin« cos 4, | dae (1.30) 


FY = —(Cy+Dp)pecosa, Fy=—(C+D)ucosa. | 
Using equations (1.25), (1.27), (1.28) and (1.30) we finally obtain 


—1e2OF(cosh?}) 1 


=> ee LUNE ES i : 9 A 
Ce 6° 9 P,(cosht,)\ (l= CS) (n +4){P/(cos y)} 


ee O (cosh, —)" 
re Scare a aD (2 ta) Feat COS By) te 


A ee (1.31) 
2 < Q,(cosh ty) 1 ; 
Cr +Do= Pave Pleat) P, (cosh i.) at (x) (n + 4){P,(cos 8o)} 
rs Q,(cosht,) (—)” ae ina | 
C+D 5D P,(cosht,) (1 =X Oye cones 


POP PE NMD lex 2 


We wish to prove that if A is a tensor given by 


(0) 
EY iP OBE 
eS ( : yop ee (1.32) 


and if f(A) is any tensor function of A, then 
[FAM =[AA)? = Ff +y) +4 fl -y) 
[FAY]? = [FA = 30 +7) -4fl(Y—y). 
n= eae y(n) ) reget sas (1.33) 


y(n) Yn) 
Multiplying the tensor A in (1.32) by the tensor A” given by (1.33) we obtain 


Bg On| 


Suppose that 


err: (1.34) 
y(n t+1)=Yy(n) +yy(n). | 
Eliminating y() from equations (1.34) we obtain 
Y(n +2) — 2y O(n +1) + (yO? —y2)yO(n) =0. ...... (1:35) 


By eliminating y(n) we obtain a similar equation for y(n). Putting y(n) a x 
we obtain from (1.35) 


= 2/0n 17/08 2a) = oe eae (1.36) 
giving 
x= y +y. 
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The general expression for y(n) is therefore 
POS Ayo ay) Ry oy); 
similarly 
Oia d ye itey )tocie ss) 
where Ay, Bo, A and B are arbitrary constants, which can easily be found to be 
A, =B, =A=-B=}. 
Therefore 


Y%(n) =H +7)" + 40/9 — 9)" = (AN) = (NY 
y(n) =4Y +7)" Hy — y= (AN)! = (NY 


Hence the required formula follows. 


AEP TON Dal Xe 3 
The tensors y, y are given by 
J) = Ades +(By—Ap)ReRe3  Vs5° = Adee +(B—A)R Ry «2.0 (1.37) 


where A,, By, A and B are constants, k is a unit vector and k,(s=1, 2, 3) are 
its components referred to an orthogonal triad. It is required to prove that if 
F(y, y) is any tensor function of Y, y then 


[Fy ee = F(ApA) 859 +[F( BoB) — F(ApA) Ry «. +. (1.38) 
We first consider the tensor (y”),.; omitting the suffixes s, s’ we have 


y =A +(B—A)(kk) 


2 ie il a a Oe ee Peer 1739 
of = Dd Ae at A)" 6) (RR) % | ( ) 
T=0 
We also have 
[(RR)?] ss i LhekekyRy =kky = (RR) sey ) 
by successive pre-multiplication by (kk) we obtain 
(kk) = (kk) (aed ig a RET AS es: (1.40) 
Substituting from (1.40) into (1.39) we obtain 
yt =A” +(kk) & A™(B—AY (") 
r=1 
pty UES IONS Dt ale Rea YE (1.41) 
similarly 
yom = An +(BR—AG)(RR) we ee ee (1.42) 
From (1.40), (1.41), and (1.42) 
Orman — A” A” + (BeBe —AGA")(RR). wwe ees (1.43) 


Hence the required formula (1.38) follows. 
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THE CORRELATION COEFFICIENT 


The equation for the dielectric constant « given by Kirkwood for isotropic 
media may be written (Kirkwood 1939) 


e— 1 meee 3e Pr’ O 1 
3 33 Deep ieasea eeenl otemonens (2. ) 
OF4[ 142 cosy) ae eee (2.2) 


where B =«/a3, f=2(e —1)/(2e +1); « is the polarizability, p the permanent moment 
of a molecule, z the coordination number, i.e. the number of first neighbours, 
and y is the angle between the two dipole directions of two first neighbours. 

We shall calculate the average of cosy by multiplying it by the probability, 
given in Part I, that two dipoles have a given configuration, and then integrating 
over the orientations of the two dipoles. Since this probability depends on « 
(because in the calculation, the action of distant molecules on a given one is replaced 
by that of a continuous distribution) Q will also depend on «, but for high values 
of « the probability depends very little on « and tends to a definite limiting value 
as « tends to oo. In the application when the formula is of interest « is very 
large indeed, therefore in the calculation of Q we can use an approximate large 
value of «, say the experimental one. By substituting this value of Q in 
equation (2.1) it becomes quadratic in e« and its positive root represents the 
theoretical value. 

The degree of coincidence of this value to the experimental value of the 
dielectric constant is a very useful means of checking the configurational 
probability and the simple model which have been used. 

Let (8, 0) and (0, 4) be the polar coordinates of the two dipoles, assumed 
to be first neighbours, then 


cosy=cos@cos@’+sin@sin@’cos¢d. —...... (2:3) 


The expression for W given by (1.17) can then be written in the form 
W 
ai a G(cos? 6 + cos? 6’) + Lcos@cos 6’ + Msin@sin6’cosP ...... (2.4) 
where 


aye sy) 
Cate 


F(BoB) — Fo(ApA)] 
p* | 

—L= SRT LPB) eee te (25) 
p? 

—~M= SAT F(A,A). 


The average of cos y can be written 


cos y=cos 6 cos &’ +sin @sin 0’ cos¢ = = aa = 9 SAV nueeee ee (2.6) 
0 0 
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where 
R= a | exp{—W/AT}sin@sin 6’ d0d0'd6 ns, (2.7) 
ee & | cos@cos@’ exp{—W/kT}sinOsin6' d0d0'db we (2.8) 
Kee = | sin 6 sin’ cosdexp {—W/RT} sin sin’ d0.d0'dp. ...... (2.9) 


Performing the integration over ¢ by use of the formula 
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aes nd exp {x cos 6} db =2n( —1)"J,(ix) =271,,(x), 
( 


where J,,(x) is a Bessel function of order m. Afterwards for the integrals over 
0, 0 we use the formula 
"mn en/2 -x/2 
| : | : F(60, 0’) do do’ = | | [F(@, 8’) + F(a — 0, 6) + F(0, 7 — 0’) 
J Q/ “0 0 
+ F(a —0, 7 —0')] dé dt" 


and obtain the following expressions for Ky, K, and K, 


rr 


1 Ses | pa : “cosh (L cos @ cos a’) exp {G(cos? 0 +cos? 0’)} 


<a (Msia@cin’ Vein? sin0-d0 dt >. ye (2.10) 


rnfl2 rn/2 
1G | | sinh (L cos 6 cos 6’) exp {G(cos? 6 + cos? 6’)} 
0 J0 
x I,(Msin 6 sin 6’) sin@ sin 6’ cos@cos6’ dbd#’ ti... se . (2.11) 


nl2 rx/2 
K,= | | cosh (L cos 6 cos 6’) exp {G(cos? 6 + cos? 6’)} 
0 Jo 
x L(M sin 6sin 6’) sin? 0 sin? 6’ dO dd’.  - ns ee (2.12) 


Application to Water 
In this section we shall calculate Q for water, on the assumption that only 
dipolar forces are essential. This seems to be rather doubtful (Pauling 1935) 
as the hydrogen bonds may influence the orientational distribution; but it will 
be shown that at least the exceptional size of « can be explained by taking into 
account only the dipolar interactions. 


The Mutual Potential Energy 
The field acting on one of the two first neighbours, included inside a cavity of 
the form of an ellipsoid of revolution is given by equations (1.4), (1.5) and (1.6), 
with the condition that the ellipsoid should coincide as nearly as possible with 
the cavity of two overlapping spheres whose dimensions are given by (1.2). 
The average distance between first neighbours c will be taken to be 
COL TN ae RR STL tar oe (2:13) 
thus.atel.—2/5_ kK. 
f 
mits 20425 a eRe ees rests (2.14) 


where 47a3/3=1/N, and N is the number of molecules per unit volume. 
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Substituting this value of a/c? in equation (1.2) we find 
M05 9235 Gi yo me at So eee (2:15) 

which determines the radius xc of each of the overlapping spheres. Putting 
cosht=...1-2, 1-3, 1-4, ... we find that for cosht = 1-3 the ellipsoid will coincide 
as nearly as possible with the cavity of the two overlapping spheres. For 
cosht=1:3 

€/2b=0-6499, Cb = 2 9625.) ae ee eee (2.16) 
Using these values of c/2b and cosht, and taking for ¢ its experimental value 
88-25 (Dorsey 1940) one obtains 

A,=0-7317; A=0:00946; B,=0-7802; —B=0-08161. 


The mutual potential energy is expressed by (2.4) and (2.5). Using the 
values provided by (2.17) one obtains 
= G=1-998 5" Eh=7 62855 = 09945) gee (2.18) 
Atel = 27 Oe 


The Correlation Coefficient Q 


The double integrals Ky, K,, Kz given by (2.10), (2.11) and (2.12) were 
evaluated numerically using the values of G, L, M given by (2.18). These were 
found to be 

Ko =(z/2)? 6:8381/81 


Ky Gri2)2 323 185 S. a ee ee (27219) 
— K, =(a/2)? 0-8578/81, 
then from (2.6) COS, = 030 0 Re ee (2.20) 
Taking zs =4 we obtain from (2.2) 
O'= 0:8433 250 “Re woke or eee (2.21) 
2 
Putting 6 =0-2248, i 1 9:6408 2 ae (2.22) 
and substituting for O from (2.21), equation (2.1) becomes 
1:20187e? —1-12373e— 2:10134=92:59288e, >) 2... (2:23) 
giving CHT So Sa ear ee eee (2.24) 
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ABSTRACT. Agreement between the results of theoretical and experimental investigations 
of the absorption of hydrogen by palladium indicates that, at least at higher temperatures, 
the hydrogen is contained in solid solution. Apparatus for observing the changes in electrical 
resistance of palladium during absorption and evolution of hydrogen is described, with 
which resistance—pressure isotherms from 75° to 150°c. were obtained for pressures 
increasing and decreasing between zero and atmospheric. Factors affecting the resistance 
of the Pd-H system are discussed and an explanation of the form of the isotherms is given; 
a decrease in resistance which is observed, when hydrogen is added to the alloy in the 
B phase, is attributed to completion of the d shells of palladium. 


I. SoLUBILITY OF HYDROGEN IN PALLADIUM 
§1. INTRODUCTION 


EW special problems have been so extensively investigated as that of the 
f occlusion of hydrogen by palladium. McKeehan (1923) gives a biblio- 

graphy of 66 titles of papers mostly published between 1900 and 1922, 
and many more are quoted by Ubbelohde in reviews of the subject in 1932 and 
1937. ‘The inconclusive nature of much of the earlier work appears to be due to 
a failure to recognize the solubility of gases in metals as simply a particular aspect 
of the general problem of pressure-temperature equilibrium of alloys. In his 
general review (1937) of the properties of metallic hydrogen and its alloys, 
Ubbelohde advanced explanations of various phenomena associated with the 
palladium-—hydrogen system in terms of the then recently developed theories of 
the metallic state. ‘The exceptional properties of transition metals generally, 
are due to the partly-filled d band. Palladium, for example, has approximately 
0-6 unoccupied states per atom in the d band, and there is now adequate evidence 
that when hydrogen dissolves in palladium, it is almost completely ionized and 
its electrons enter the d band. In a statistical interpretation of the solubility, 
Lacher (1937) calculates, from measurements of equilibrium isotherms at various 
temperatures, that the number of potential ‘holes’ per Pd atom which are 
available for hydrogen is nearly constant and equal to 0-6 approximately. From 
this result, it is to be inferred that the process of hydrogen absorption reaches 
completion when the d band is filled, that is, when the ratio r of H to Pd atoms 
is approximately 0-6. When palladium is charged with hydrogen electrolytically, 
however, values of 7 greater than 0-6 are observed at temperatures below about 
80°c. At 0°c., for example, the maximum concentration approaches unity. 
If the hydrogen is entirely located interstitially in solid solution, the additional 
absorption of hydrogen may occur with entry of electrons into the s band of 
paliadium. Such a possibility is assumed by Wagner (1944) in calculating the 
solubility from thermodynamic principles. It is also possible that the formation 
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of simple solid solutions is not the only absorptive process operating. Ina 
recent note (1948), D. P. Smith presents a formidable body of experimental 
evidence in support of his view that a large proportion of the occluded hydrogen 
is concentrated in rifts or widened spacings in the palladium lattice. Assuming 
at 0°c. a total concentration r1, of which the concentration in solid solution 
is 0-6, we find that the metal occludes about 500 times its own volume of hydrogen 
(measured at normal pressure) by some additional absorptive process. If 
Smith’s view is correct, we may suppose that when palladium is charged with 
hydrogen electrolytically, the gas (probably in a highly ionized state (Moore 1939)) 
is concentrated under enormous pressure in the intergranular fissures in the metal. 


§2. THEORETICAL FORMULAE FOR THE SOE OA MEANY 
Lacher’s formula is developed by assuming that a definite number of holes 
in the metal are available for hydrogen and that the heat of absorption increases 
with the number of holes filled. It gives a critical temperature above which 
only one solid phase can exist and below which two solid phases can be in 
equilibrium with one another. The equation for the concentration—pressure 
isotherm at any temperature T may be expressed in the form : 


log p¥(atm.) = log 5SoTy ~2-3009 — (445-6 + 1672r)/T. 


Over a certain range of pressure on any isotherm there are three roots to this 
equation, one of which, r=0-59/2, corresponds to an unstable state at the critical 
pressure at which the phase transformation occurs. Ifr, and rg are, respectively, 
the concentrations in the « and f phase at this critical pressure, for equilibrium 
conditions these roots must satisfy 7,+7,;=0:59. The 100°c. isotherm, 
calculated from the above formula, is shown in Figure 1. The continuous line 
is the equilibrium isotherm, and by putting r=0-59/2, the value 28-1 cm. Hg is 
obtained for the critical pressure at which the phase change takes place. 

Over the greater part of the range of temperature investigated experimentally, 
Lacher’s formula gives a satisfactory representation of the ‘equilibrium’ 
pressure—composition isotherms for the Pd—H system. It must be observed, 
however, that the term ‘equilibrium isotherms’ as used above applies to those 
isotherms which give the same value to the pressure for a given value of 7 no 
matter how that value ofvis approached. At higher temperatures (200° to 313°c.), 
Gillespie and Galstoun (1936) obtained such reversible isotherms without using 
any special preliminary heat treatment. Wagner contends that because of the 
inadequate mobilities of the Pd atoms at the comparatively low temperatures 
at which some isotherms have been determined, it is not legitimate to assume 
that complete thermodynamic equilibrium is established. He points out that in 
the experiments of Gillespie and Hall (1926), for example, a preliminary heat 
treatment was used in order to obtain definite single values of the pressure when 
the « and f phases co-exist, and that the values thus determined approximate 
to those for complete equilibrium probably because they lie within a ‘ hysteresis’ 
loop. In most experimental investigations it has been found that, over a certain 
range of pressure, absorption and evolution isotherms do not coincide. This 
hysteresis phenomenon indicates that thermodynamic equilibrium has not been 
completely established, and is ascribed by Wagner to restricted exchange of 
Pd atoms between the two phases at the comparatively low temperatures of 
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observation. In expressing the conditions of equilibrium by using Gibbs’ 
equation, he observes, therefore, that whereas equality of chemical potentials 
is true for transitions of H atoms between the two phases, the analogous condition 
for transitions of Pd atoms is not applicable. From this new viewpoint, Wagner 
offers a qualitative explanation of the experimental results and shows that the 
above assumptions lead also quantitatively to good agreement as to the width of 
the hysteresis loop. The solubility isotherms are expressed by the formula 


: Yr 
VE >. Kexp (Ay/RT) 


in which K is the solubility constant for ideally dilute solutions, ie. K =4/p/r 
for r<1; and Au is a function of r calculable from experimental isotherms for 
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Figure 1. Concentration—pressure isotherms at 100° c. 
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the single-phase region above the critical mixing temperature. For calculating 
the isotherms at lower temperatures, Aj must be assumed independent of tem- 
perature and K must be known. In calculating the 100°c. isotherm (Figure 1) 
from the above expression, K was chosen by trial to give the best representation 
of the experimental results of Sieverts and Danz (1937). Bearing in mind the 
idealized assumptions necessarily involved in deducing the solubility (compli- 
cations such as the effect on the equilibrium pressure of hydrogen of strains in 
the regions bordering the two phases are ignored), the correspondence with 
experiment is remarkably good. 

PROC. PHYS. SOC. LXIII, 7—A 40 
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§3. EXPERIMENTAL STUDIES OF THE HYSTERESIS BEBE Car 


An examination of experimental data reveals that reversible isotherms have 
been obtained only by using palladium in a finely divided state. With palladium 
in the form of filings, sponge, wire, etc., it has been found that the isotherms 
for increasing and decreasing hydrogen pressure do not coincide. Sieverts and 
his co-workers (1933, 1935, 1937, 1938) have also used palladium ‘black’ with 
similar results. This hysteresis effect has been observed in X-ray investigations 
of the expansion of the hydrogen-charged palladium lattice. Detailed studies 
by Owen and Jones (1937) and by Owen and Williams (1944) establish the 
presence of two co-existent expanded lattices occurring over different ranges of 
pressure, and at different mean pressures, on ascending and descending isothermals. 
In a later publication, Owen (1944) discusses the results of the x-ray work in the 
light of Lacher’s hypothesis that there exist metastable states of supersaturated « 
and undersaturated 8 phase in the palladium—hydrogen system. The experi- 
mental evidence clearly indicates that the mixed-phase or (« +f) regions within 
the hysteresis loop are characteristic of a system in which thermodynamic 
equilibrium is not completely established. The X-ray studies also indicate 
that the new phase, « or B, exists finely dispersed throughout the old phase, 
the establishment of a new phase occurring, presumably, by the process of 
nucleation and growth. 

The above conclusions support the assumptions on the basis of which Lacher 
explains the hysteresis in the pressure—concentration isotherms. Smith questions 
the validity of these assumptions and seeks to account for the behaviour of the 
Pd-H system in terms of occlusion in rifts as the principal absorptive process, 
rather than solid solution. Experimental evidence mentioned by Smith appears 
to show that the « and £ phases do not exist interdispersed but are separated by 
a single-phase boundary which migrates inwards from the surface of the metal 
as the concentration of hydrogen either increases or diminishes. The present 
work (II, §2) indicates that, at low temperatures, the (« +) phase exists over an 
extremely small range of pressure and, further, that the transition from one pure 
phase to the other takes place very rapidly when the pressure change is sufficient 
to carry the system over the two-phase region. All the experiments referred to 
by Smith are conducted at room temperature, the palladium being charged by 
making it the cathode of an electrolytic cell. Under these conditions it is probable 
that the interdispersion of phases would escape observation. 

Pressure—concentration isotherms obtained by Sieverts and his co-workers 
show that the width (pressure range) of the hysteresis loops increases with rise 
of temperature up to about 220°c., and then decreases, becoming zero at about 
300°c. Lacher’s formula predicts this feature of the isotherms, though it gives 
a greater width to the hysteresis loops than is observed experimentally. From 
Lacher’s hypothesis of the persistence of metastable states, it follows that the 
range of metastability increases with temperature up to 220°c. ‘This conclusion 
has been criticized, it being assumed on general grounds that the system should 
attain a state of equilibrium more readily at higher temperatures. It must be 
pointed out that the pressure at which the phase transformation occurs increases 
with rise of temperature, and that the persistence of the metastable state should 
be considered also in relation to this pressure variable. If P, and P, are, 
respectively, the upper and lower pressure limits of the hysteresis effect at any 
temperature 7’, we might expect the ratio (P,—P,)/$(P,+P,), rather than the 
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value of (P;—P,), to diminish as T increases. Both Lacher’s formula and the 
experimental results of Sieverts and others show that the ratio does, in fact, 
diminish rapidly as the temperature rises. 

Lacher’s calculation of the solubility, whilst using the value r =0-59 as a limiting 
hydrogen concentration set by the completion of the palladium d band, does not 
exclude the possibility of higher concentrations due to the occupation of other 
regions of energy. Such a possibility is assumed in Wagner’s thermodynamical 
theory which, as we have seen, gives a remarkably good estimate of the widths 
of the hysteresis loops. The substantial measure of agreement between theory 
and experiment may be taken as an indication that, at least at higher temperatures, 
solid solution is capable of accounting for the occlusion of hydrogen by palladium. 

The primary object of the present investigation was to discover if the hysteresis 
eflects in the electrical resistance of Pd—H could also be explained by assuming 
the hydrogen to be held in solid solution. In the experiments described below, 
the change in resistance of palladium charged from the gas phase was studied 
under carefully controlled conditions, and at temperatures lower than those 
covered by Briining and Sieverts, particular attention being given to phenomena 
at 100°c., for which temperature x-ray and solubility data were available. 


II. ELectTricaAL RESISTANCE OF PALLADIUM—HYDROGEN 


$i. EXPERIMENTAL 


Resistance was measured by the potentiometer method, all the usual precautions 
being taken to avoid errors due to thermo-electric effects. A Cambridge vernier 
potentiometer and D’Arsonval galvanometer were used and the unknown resist- 
ance measured in terms of standards from a multi-range current shunt. A 
large capacity 2-volt accumulator maintained an extremely steady current of 
about 1 ampere through the two resistances and a series-connected ballast 
resistance of rolled manganin wire. Under the least favourable conditions, the 
resistances of the palladium wire specimens could be determined to within 0-1°%. 
A 15 cm, length of thick palladium wire was clamped at both ends to small blocks 
of copper mounted on a slab of steatite. The current leads were thin copper 
rods screwed into the clamping blocks, and the potential contacts to the wire 
consisted of sharp steel knife-edges about 10 cm. apart. ‘The latter were riveted 
to copper plates which were screwed down so that the wire was gripped between 
the knife-edges and the steatite base. 

Temperature was measured with a platinum—rhodium thermocouple, the 
junction of which was flattened and clamped to the palladium wire at its mid-point. 
Alumel-chromel couples, also in contact with the wire, were located between 
each knife-edge and clamping block; these were used to measure the temperature 
gradient along the wire. ‘The whole system, enclosed in a quartz tube, occupied 
a central position in a tubular furnace which had been specially wound to minimize 
temperature gradients. Gradients and fluctuations were further reduced by 
two copper cylinders separated by a layer of magnesium oxide. These filled 
the remaining annular space between the quartz tube and the wall of the furnace. 
By means of a small electric heater situated inside the quartz tube and fed from 
the furnace supply voltage through a transformer, the gradient along the palladium 
wire was finally reduced to less than 0:05°c.percm. ‘The temperature of the wire, 
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as measured by the platinum-—rhodium couple, was maintained constant to within 
0-2°c. for periods of over a month by a continuous-control thermo-regulator 
(Wright 1947) operated by a platinum resistance thermometer. The platinum 
wire was wound between thin sheets of mica on the outer of the two copper 
cylinders mentioned above. 

Before charging with hydrogen, the palladium wires were annealed 7m vacuo 
at 400°c. for about twelve hours, and their resistances measured at various 
temperatures between 400°c. and room temperature. Wires of pure palladium, 
of diameters 0-0805 and 0-210 cm., supplied by Johnson Matthey and Co. Ltd., 
were found to have an average temperature coefficient of 0-0036 per °c., and a 
resistivity of 10-0 ohm.cm. at 0°c. In the case of wires which had been charged 
with hydrogen at 100° and 75°c., it was found that a process of annealing and 
degassing at 300°c. for several hours was sufficient to restore the wires to their 
original state as regards electrical resistance. 

Pure dry hydrogen was prepared in the manner described by Owen and Jones, 
and similar measures were taken to avoid contamination of the gas in the apparatus. 
Any pressure of hydrogen up to one atmosphere could be rapidly established 
in the observation tube by placing it in communication with a large reservoir 
containing hydrogen at the desired pressure. Normally, the reservoir was left 
in communication with the tube to reduce pressure changes due to absorption 
and evolution of gas by the palladium. After measurement of the resistance of 
the annealed specimen 7m vacuo at a known temperature, the standard procedure 
at the beginning of an experiment consisted in raising the temperature of the 
palladium to within about 5°c. of that finally required, and then admitting 
hydrogen to a pressure of about 3 cm. Hg. A slight re-adjustment of the 
thermo-regulator was then sufficient to hold the palladium at the required 
temperature. Subsequent changes of gas pressure could be made rapidly, 
since they produced only a small momentary fall in temperature. Pressures 
were measured to the nearest millimetre with a Bourdon type gauge. 


§2, RESULTS 


Resistance—pressure isotherms for pressures of hydrogen up to 73 cm. Hg are 
shown in Figure 2, and the observations at 75°c. and 100°c. are recorded in the 
‘Table in the order in which they were taken. 

With the exception of those marked with an asterisk, the figures represent 
constant final values of resistance corresponding to each value of the hydrogen 
pressure. In the case of the 100°c. isotherm, for example, the final resistance 
was found by plotting a resistance-time graph after each alteration of pressure. 
For increments of pressure in the range of ascending pressure 0 to 31 cm. Hg the - 
resistance attained its maximum value in less than 5 hours. This part of the 
isotherm is reversible, that is, it can be re-traversed by reducing the pressure 
of hydrogen. ‘The second part of the 100°c. isotherm, extending from 31 to 
about 55 cm. Hg pressure, is irreversible. The observations in this range of 
ascending pressure are plotted against time in Figure 3, and the values quoted 
in the Table are the approximate maxima obtained from these resistance-time 
curves. After each pressure change, the resistance reaches an apparent maximum 
in a few days. ‘This is followed, however, by a much slower increase. At some 
pressures, readings were taken daily for over a month; though a gradual 
diminution in the rate of increase of resistance was observed, the resistance was 
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R/Ro=ratio of resistance of Pd—H to that of pure Pd, at T° c. 
p=pressure of hydrogen in cm. of mercury. 


* Resistance changing ~0:1% per day. 


Pressure of Hydrogen (cm.Hg) 


Figure 2. Resistance-pressure isotherms for the palladium—-hydrogen system. 
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still increasing at the end of the month. ‘The constant-pressure sections of the 
curve a in Figure 3 terminate at the stage when the resistance 1s increasing at 
a rate of less than 0-1°% per day. 

The third part sone eae ascending pressures 55 to 73 cm. Hg and 
descending pressures 73 to 17 cm. Hg, is reversible, the maximum resistance 
values being attained in less than 2 hours after altering the pressure. The 
fourth part, descending pressures 17 to about 8 cm. Hg, is irreversible. Here 
the experimental procedure corresponded to that adopted in obtaining the 
second part of the isotherm, except that the pressure changes were smaller. 
The resistance-time curve, c, is shown in Figure 3, each section being plotted 
to the stage when the rate of decrease is less than 0-1% per day. ‘The fourth 
part of the 100°c. isotherm re-joins the first part at about 8 cm. Hg pressure, 
follows it to zero pressure, and (without annealing 7m vacuo) re-traverses it with 
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Figure 3. Curves showing change of resistance with time. Initial pressure changes (in cm. Hg) are 
indicated. 
a. (a+) to B, 100°c. b. (a+) to B, 75°C. c. (a+ 8) to a, 100° c. 


ascending pressures. The resistance-pressure isotherm at 75°c. exhibits the 
same general characteristics as that at 100°c. In all parts of the cycle, however, 
the system was slower in attaining maxima (or minima) of resistance following 
a pressure change. At 125°, 150° and 200°c., the maxima were reached more 
rapidly than at 100°c. 

The resistance-time curves at 100°c. shown in Figure 4 were obtained in 
traversing the hysteresis loop by making four pressure changes, each of which 
was just sufficient to carry the palladium-hydrogen system through one of the 
four well-defined stages previously considered. In the reversible ranges of 
pressure, curves showing growth and decay of resistance are given. ‘The curves 
covering the irreversible pressure ranges were determined at pressures necessary 
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to take the system to the next reversible part of the isotherm, hence in each case 
the resistance acquires a constant final value. It will be seen that the irreversible 
parts of the isotherm are thus crossed in a few hours, a longer time being taken 
on the descending part than on the ascending part. 

The (R/Ro, p) isotherms at 75° and 100°c. are the same in general form as 
those determined at 180°c. and higher temperatures by Briining and Sieverts. 
The lowest temperature isotherm, that at 160°c., determined by these investi- 
gators, rises less steeply on the first irreversible part. From the experiments 
described above, it is now clear that the slope of this part of the isotherm is 
appreciably affected by the time of exposure to hydrogen. In all of these 
experiments, the same palladium specimen, a wire 0:0805 cm. in diameter, was 
used. A brief resumé will now be given of some additional results obtained 
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Figure 4. Curves showing change of resistance with time, at 100°c. Initial pressure changes 
(in cm. Hg) are indicated. 


1:5<= 29-0, a=(a+f); 29:0 + 60:9, (a+ 8) >; 
60-9 = 20-0, B=(a+B); 20:0> 6°8, (a+B) > 4. 


by using a thin palladium rod, 0-210 cm. in diameter and 15 cm. long. From a 
knowledge of volumes of apparatus and palladium, and observation of pressures, 
a rough determination of the concentration—pressure isotherm at 100°c. was 
first made. The values of r were found to be higher than those obtained by 
Sieverts and Danz (1937) who used palladium black. ‘The high values were 
probably due to absorption of hydrogen by various materials other than palladium 
inside the apparatus. In other respects the two isotherms were very similar, 
and:are in agreement as to the pressure limits of the hysteresis loop. A careful 
re-determination of the resistance-pressure isotherm at 100°c., using the 
palladium rod, was next carried out. This isotherm was found to be identical 
with that obtained with the wire, though a decidedly longer period of time was 
required to obtain the readings in the case of the rod. Comparing the (7, p) 


736 P. Wright: 

and (R/R,,p) isotherms at 100°c., it is found that the pressure limits of the 
hysteresis are the same for concentration and resistance, in agreement with the 
observations of Briining and Sieverts at higher temperatures. Both the 
concentration and the resistance measurements just considered indicate that 
the pressure limits of the hysteresis are definite and independent of both the 
form and the bulk of the palladium. 

No detailed analysis of the form of the resistance-time curves in Figures 3 
and 4 has yet been made. Over the reversible stages of the cycle at least, these 
curves are of an exponential character. Some accurate observations at 100° Cc. 
with the palladium rod were made at constant hydrogen pressures of 8-7 and 
23-3 cm. Hg, the initial pressure changes being, respectively, 23-3 to 8-7 cm. Hg 
and 8-7 to 23-3 cm. Hg. If &R, is the difference between equilibrium values of 
resistance of palladium—hydrogen at the two pressures, and SR the change of 
resistance in ¢ minutes after making the initial alteration of pressure, it is found 
that 5R=5R,(1—e™), approximately, where « =0-004 for decreasing resistance, 
and 0-006 for increasing resistance. Using Sievert’s values for the concentration 
of hydrogen in palladium at 100°c., we find in the single-phase region under 
consideration, the relation (R/R,—1)=ar, which expresses a proportionality 
between increase of resistance and concentration in the « phase also observed 
previously by Briining and Sieverts at higher temperatures. We may, therefore, 
re-write the resistance-time equation as 


Ing=A—Bt 


where g now represents rate of change of either concentration or resistance, 
and A and B are constants. For a gas initially uniformly distributed in a homo- 
geneous wire of circular section, a similar expression has been deduced, from 
Fick’s linear diffusion law, for the rate of desorption (Smithells 1937). 


S35 DISCUSSION OR RES UES 
(1) Hysteresis Phenomena 


The experimental values for the pressure limits of the hysteresis at 100° c. 
are about 20 and 35 cm. Hg according to Sieverts and Danz, as compared with 
about 17 and 33 cm. Hg indicated by the present measurements of the electrical 
resistance of the system. Both sets of pressure values are read from the most 
steeply rising and falling parts of the isotherms. It is evident from preceding 
changes of curvature, however, that the new phase is in the process of formation 
before the pressure limits are reached. ‘This is borne out by x-ray analysis 
(Owen and Williams 1944), which reveals co-existing « and £ phases in the range 
of ascending pressures 27 to 41 cm. Hg, and in the range of descending pressures 
19 to 13 cm. Hg. We must now form a picture of the hysteresis cycle which is 
consistent with all the experimental data. As the pressure is increased from 
zero, the « phase at first exists alone. The 8 phase begins to form probably 
before a pressure of 27 cm. is reached, and up to about 31 cm. Hg it disappears 
again if the pressure is reduced (the system is reversible in the range 0-31 cm. Hg). 
If the pressure is increased beyond 31 cm. Hg, the system transforms irreversibly 
into the 8 phase, and the rate of transformation increases very rapidly with the 
pressure. Analogously, with pressures descending from atmospheric, the 
8 phase at first exists alone, and the « phase appears before the pressure falls to 
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a9 cm. Hg. ‘The system is reversible between atmospheric and about 17 cm. Hg 

and if the pressure falls below 17 cm. Hg the « phase is established inreversibly, 
Assuming that nuclei of a new phase are formed dispersed throughout the old 
phase, it Is possible, as Lacher has shown, to explain the hysteresis by considering 
the contribution to the free energy of the system made by the surface separating 
the two phases. The vapour pressure of the dispersed £-phase particles would 
be greater than the pressure they would exert in bulk, hence they could not 
form without some degree of supersaturation in the « phase. Similarly, during 
evolution of hydrogen, an undersaturated ~ phase would be necessary for the 
formation of «-phase nuclei. Over a certain range of pressure up to the point 
where the isotherms rise (or fall) sharply, it must be assumed that the two phases 
co-exist in a state of metastable equilibrium. Beyond this stage, the degree of 
supersaturation in the « phase (or undersaturation in the B phase) is such that the 
system becomes unstable and the new phase is rapidly formed in fairly large 
blocks out of the finely dispersed particles. This version of the behaviour of 
the Pd-H system is supported by the x-ray evidence and appears to be not 
inconsistent with the results of solubility measurements. We shall, therefore, 
proceed to use it as a working hypothesis in an attempt to account for certain 
characteristics of the resistance isotherms. 


(ii) Origin of Resistance Changes 

When hydrogen is alloyed with palladium by the process involved in deter- 
mining the resistance—pressure isotherms, the palladium lattice retains its face- 
centred cubic structure throughout. In transforming from the « phase to the 
hydrogen-rich or f phase, and again during the transformation back from B phase 
to « phase, the new phase appears as a conglomeration of small crystals. Crystal 
size is not, however, an important factor as regards the electrical resistance since 
the present results (II, §2) show that the resistance-pressure isotherms can be 
repeated without any initial annealing of the palladium. Hydrogen in solid 
solution may influence the resistance of palladium by disturbing the periodicity 
of the lattice or by contributing electrons to fill the incomplete d band of 
palladium. From the close correspondence between concentration and 
resistance isotherms we may deduce that the scattering centres provided by the 
hydrogen atoms are mainly responsible for increasing the resistance. The 
observed increase of resistance brought about by extensive cold-working of pure 
metals is small compared with the effect of alloying, hence we should expect 
lattice distortion at phase boundaries to be a relatively unimportant factor. 
The increase in resistivity introduced by it during the formation of a new phase 
is likely to be a maximum when the particles are very small, since the scattering 
should be a maximum when the wavelength of the conduction electrons is 
comparable with the dimensions of the particles (Mott 1940). The high 
resistance of the transition metals is attributed to scattering processes in which 
the electrons make transitions into vacant states in the d band. ‘The transition 
probability is high because of the high density of states in the partially filled 
d band of the transition metals. It has, in fact, been verified experimentally 
that the electrical resistance decreases when, as a result of alloying, the incomplete 
dshells are filled. Thus, Rosenhall (1935) has shown that the addition of hydrogen 
to palladium-silver alloys reduces the total resistance. Since the density of 
states in the d band is proportional to the cube root of the number of positive 
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holes, we should expect the reduction in resistance due to alloying with hydrogen 
to be apparent in the 6 phase of palladium—hydrogen when the hydrogen 
concentration r approaches the value 0-6. 

Consider the resistance-pressure isotherm at 100° c. As the gas pressure 
increases from zero, the « phase at first exists alone and the total resistance, 
which is here proportional to 7, consists principally of the scattering resistance 
due to hydrogen atoms. When the pressure approaches 27 cm. Hg, the B phase 
begins to form and the total resistance is slightly augmented as a result of strains 
at the phase boundaries. ‘The magnitude of this contribution depends upon 
the number and size of the f-phase nuclei and would probably increase with 
pressure up to 31 cm. Hg. At about 31 cm. Hg the system becomes unstable 
and transforms completely into the 8 phase, the rate of evolution increasing rapidly 
with the pressure. As mentioned previously, the slope of the isotherm between 
31 and about 55 cm. Hg depends upon the time of exposure of the palladium to 
hydrogen. If sufficient time were allowed, this part of the isotherm would be 
parallel to the resistance axis and would join the evolution curve at a pressure 
just exceeding 31 cm. Hg. A comparison of solubility and resistance isotherms 
in the pressure range between atmospheric and about 17 cm. Hg, reveals that a 
marked increase in resistance occurs as the hydrogen concentration diminishes. 
With descending pressures, the 8 phase does not form until the pressure approaches 
19 cm. Hg, hence the increased resistance cannot be due to strains but must be 
attributed to the increase in the number of positive holes in the d band of 
palladium. If the isotherm is re-traversed with pressures ascending from 
19 cm. Hg, to atmospheric, 7 increases from about 0-53 to 0-60, and R/R, 
diminishes from 1-57 to 1:53 as the positive holes in the d band become filled. 
If this explanation is correct, since the concentration of hydrogen in the f phase 
is less at higher temperatures, we should expect the effect to diminish with rise 
of temperature. The isotherms of Briining and Sieverts show that the expected 
diminution does actually take place. Solubility data at 75°c. are not available; 
at 80° c. the saturation value r is greater than 0-6. ‘The still greater value of r at 
75°c. may account for the fact that the res’stance isotherm at 75°c. shows a 
slight rise as the concentration in the 8 phase approaches its maximum value. 

The upward trend of the resistance isotherm at 100°c. terminates abruptly 
as the pressure descends to about 17 cm. Hg. Here the system becomes unstable 
and transforms completely into the « phase. With the evolution of hydrogen, 
the total resistance drops when the rise due to removal of electrons from the 
d band is outweighed by the fall due to reduced hydrogen scattering. ‘The 
experimental points on the descending curve between 17 cm. and 6 cm. Hg 
do not represent constant values of the resistance. As in the case of the trans- 
formation from « to 8 phase, the change should be represented by a line parallel 
to the resistance axis. 

At the temperatures considered, we conclude that the general characteristics 
of the resistance isotherms in the system Pd-H can be satisfactorily explained 
on the basis of Lacher’s suggestions, and by assuming that the occluded hydrogen 
is held entirely in solid solution. ‘The incomplete d shells of the Pd atoms 
become filled when the ratio, 7, of H to Pd atoms in solid solution is approximately 
0:6. Thus, the form of the isotherms in the 8 phase may be accepted as evidence 


in support of the theory that the high resistance of a transition metal is due to 
its partly filled d band. 
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The Barkhausen Effect 
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19th December 1949 


ABSTRACT. An account is given of further experiments on the Barkhausen effect in 
cylindrical ferromagnetic specimens. The form of the decay of induction following a 
Barkhausen discontinuity has been investigated and there is reasonable agreement between 
the theory given and the experimental results. As a result of this investigation it has been 
possible to satisfy the different conditions which are necessary for the detection of the 
discontinuities so as to obtain (a) accurate reproduction of the form of the discontinuities 
and (b) maximum sensitivity in the counting experiments. The number and size of the 
discontinuities are measured, and an estimate is made of the contribution of the discontinuous 
processes to the total change in magnetization over the relevant part of the hysteresis curve. 
Discontinuous changes in magnetic moment greater than 0-3 x 10° E.M.U. (corresponding 
to a minimum volume of 0-9 x 10- cm? for iron) account for from 47% to 86% of the 
total change in magnetization. This is considered in relation to associated work on revers- 
ible changes, and the nature of a Barkhausen discontinuity 1s briefly discussed with reference 
to other investigations on domain structure. 


§1. INTRODUCTION 
N a paper by Bush and Tebble (1948) (which will be referred to as I), an 
account is given of measurements on the contribution of the discontinuous 
processes, known as the Barkhausen effect, to the total change in magnetization 
as a ferromagnetic is taken through a hysteresis cycle. It was possible to detect 
discontinuities corresponding to reversals of magnetization in volumes of the 
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order of 10-° cm? and it was found that over the region of the magnetization 
curve considered, 10 to 30°, of the total change in magnetization could be 
accounted for by the Barkhausen discontinuities detected. Since reversible 
changes, not detected by the apparatus, contribute only 10 to 15%, of the change 
in magnetization, there was a considerable discrepancy. ‘This could have 
arisen from the presence of discontinuous changes too small to be detected, and 
an improvement in sensitivity was, therefore, considered desirable. 

The method used, both in the work described here and in the previous paper, 
is as follows. ‘The magnetic field, applied in a direction parallel to the axis of a 
cylindrical specimen, is slowly changed and the resulting discontinuous changes 
in magnetization are detected as voltage pulses, induced in a search coil wound 
round the specimen. These voltage pulses are amplified and can be displayed 
on an oscilloscope or passed to a suitable recording device as indicated in Figure 1. 
From the amplitude and shape of a pulse, an estimate can be made of the 
corresponding change in magnetization in the ferromagnetic material. By 
means of an amplitude discriminator and counting unit connected to the output 
of the amplifier, the number and size of these discontinuities can be measured 
as the specimen is taken through part of a hysteresis loop. Hence an estimate 
can be made of the contribution of the Barkhausen effect to the total change 
in magnetization over that part of the hysteresis curve examined. 


Pulse “ 
Generator 


ee Coil 
Search Coil. Specimen| 


Field Changing 
Circul 


Counting 
Unit 


Figure 1. Block diagram of the apparatus. 


The main problem in this method of approach is that the amplitude of the 
induced voltage pulses is small, of the order of a few microvolts. It is thus 
essential, if the smaller discontinuities are to be detected, that the search coil 
and amplifier should be designed to give the maximum possible signal-to-noise 
ratio; and it has been found necessary, in view of the lack of information on the 
subject, to carry out a thorough examination of the factors affecting the relation 
between the primary changes and the measured effects ($2). This involves a 
theoretical investigation into the rate of decay of induction in cylindrical specimens 
following a discontinuous change in magnetization (§2(ii)), and into the effect 
of the search coil on the form of the induced voltage pulse (§2 (ii)). The effect 
of thermal agitation noise and stray capacities on the sensitivity of the apparatus 
is also considered (S§2 (iv) and 2(v)). The conclusions which can be drawn 
from this examination have been applied to the problems arising in two distinct 
sets of experimental investigations. In the first (§3), which is subsidiary to the 
main purpose of this work, measurements have been made on the time constant 
of the discontinuities and the results are compared with those given by the 
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treatment in §2. ‘The improved apparatus has been used to carry out further 
experiments on the contribution of the Barkhausen discontinuities to the total 
change in magnetization, and the results are discussed with reference to present 
ideas on domain structure (§4). Although the contents of §§ 2 and 3 are necessary 
to justify the methods used, a detailed study of these sections is not necessary 
to an understanding of § 4, in which the main investigation is described. 


v2, LH EO RY, 
§2 (1). Introductory 


It is convenient to consider briefly the possible mechanism by which a 
Barkhausen discontinuity is produced. Ina ferromagnetic in low fields the primary 
effect of any increase in the applied field is to cause a movement of the boundaries 
of the domains, such that those magnetized with components parallel to the 
field expand, absorbing those less favourably magnetized. As the field is changed 
a given boundary is displaced, at first reversibly, until at a critical field it moves 
spontaneously (1.e. without further increase of field) over a considerable fraction, 
or even the whole of the length, of a neighbouring domain. 

It was for some time thought that such an irreversible boundary movement, 
from one energy minimum to another, corresponded to a Barkhausen 
discontinuity. The extent of these boundary movements must be limited to 
some degree by the presence of the usual irregularities in the structure of the 
material; the results of the work of Williams and Shockley (1949) on powder 
patterns on single crystals suggest that the Barkhausen effect may be related, 
at least in part, to the secondary domain structure round these imperfections. 
However, in polycrystalline specimens, particularly if the grain size is small, 
the distinction between primary and secondary domains may be less sharp, 
and the Barkhausen effect in such specimens might be a guide to the structure 
of the material. 

The case will now be considered of a cylindrical ferromagnetic specimen 
in a field applied in a direction parallel to the cylinder axis. A Barkhausen 
discontinuity in such a specimen, resulting from a small change in the applied 
field, can be considered as a change in magnetic moment AM, resolved in a 
direction parallel to the applied field, over a small volume v. This results in 
a change in the flux of induction over the cross section of the specimen given by 


4AM 
ee Che 2 pare fin) 


where d is the length of the region of the discontinuity. 
The volume v is given by 
AM 
v= =—— 
21, cos 6* 
where J, is the saturation intensity of magnetization, @ is the angle between the 


direction of J, and the applied field, and the volume is a minimum, v,,, when 


cos) = 13 
~The total flux through the search coil, and the time-varying voltage E(t) 


induced in the coil, are related by 


ee) eae ee na (2) 
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with appropriate limits of integration. ‘The amplitude and form of the voltage 
pulse will, therefore, depend, inter alia, on the rate of change of induction in 
the specimen. This, together with the original boundary movement, is 
controlled by the eddy currents set up in the specimen, and it is necessary, 
therefore, to consider first the factors controlling the decay of induction in a 
cylinder. 
§2 (ii). Decay of Induction accompanying a Barkhausen 
Discontinuity in a Cylindrical Specimen 

Bozorth and Dillinger (1930), in their work on the Barkhausen effect, assumed 
for the change in flux of induction over the cross section of the specimen an 
expression given by Wwedensky 

$ =4bou ee A Mee fhe AES eS, (3) 

1 
where A, is the nth root of J,,(A) =0, and € =t/47pa°o, o being the conductivity 
of the specimen, pz the reversible permeability, and a the radius of the specimen. 
This gives a ‘time-constant’, tw, in which the flux falls to 1/e of its total change 
at € =0-103, i.e. when 
bg = 0 aT (4) 

For a specimen of radius 0-2 mm., » =100, resistivity 1/o=10-° ohm cm., this 
gives Tw =5-4 microseconds. 

This expression, however, was originally derived by Wwedensky (1921) for 
the decay of induction following a sudden change in applied field, i.e. the initial 
change in induction takes place over the whole cross section of the specimen. 
This is not quite the same problem as in the Barkhausen effect, where there is an 
initial localized discontinuity that can take place at any position in the specimen. 

In the treatment given below, the following assumptions are made: 

(a) The area of cross section over which the initial boundary movement 
takes place is small compared with that of the specimen. ‘This holds, in the 
conditions of the present experiments whatever the precise process involved. 

(b) The movement of the boundary takes place in a time which is short 
compared with the time taken for the decay of the resulting change in induction 
throughout the specimen. ‘This is justified later. 

(c) The permeability is constant. (Since a single isolated boundary movement 
is being considered, the resulting changes in induction in the specimen are 
reversible i.e. it is reversible permeability which is involved. If this were not 
so, it would imply the occurrence of further discontinuities during the incident 
under consideration. ‘This reversible permeability is taken as constant during 
the whole process associated with a given single. discontinuity, but it will, of 
course, depend on the magnetic state of the specimen, i.e. on its position on the 
magnetization curve.) 

(d) 'The applied field is changing so slowly that it remains sensibly constant 
during what may be called the single discontinuity process. 

Consider a discontinuous change in magnetic moment taking place over a 
small area AS about. the-point-P, @=6, 0 =Uj22 —2)) seen Figures 2.enetiie 
cross section of the specimen z=2,. This corresponds to a discontinuous 
change in induction AB, over an area AS, and the change in flux over the whole 
of the specimen is given by 


bo =AB,AS 
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which will be retarded by the eddy currents set up about P. 
The equation for the propagation of this change in induction B is in 
cylindrical coordinates, 
PBS LORS 1 PB 0B 6 
BP Poe Pe apt oy, Sh. Bo 2 
where y =470, « and o are as defined above, and p is assumed constant. 
The solution of this equation is 


Da A Jaeconmd, exp, (RA) |e acess (7) 
k m=0 
where k, A,,,, are constants and m=0, 1, 2,3.... 
The boundary conditions are: (i) for all values of 6 and for ¢>0, B=0 
for r=a; (ii) att =0, B=0 except over AS at (6,0) where B=AB,. With these 
conditions, the solution (7) takes the form 


2b SS Im(Anb/a) 
Be pat) yy Ss st Ts In i\ 
TA” n=1 m=0 bF? 1 5(A,,) ( 


where € =t/4mpa"o and 6 =1 for m>0, 6 =2 for m=0. 


Tia)cos me -exp (=A,7E), io. 'see (8) 


rh 


Figure 2. 


The flux of induction over the cross section of the cylinder 1s 


Qa ra Oy) (A,,6/a) exp( —A,€) 
ne eae oe a 


0 
Thus the flux changes in a manner dependent on the distance of the initial 
discontinuity from the axis of the cylinder 6. ‘This is shown in Figure 3, where 
d/o is plotted against € for different values of b/a, together with the curve obtained 
from Wwedensky’s equation (3). In the case of a discontinuity on the cylinder 
axis, ¢ falls to d,/e when  =0-253, i.e. the ‘time constant’ is 


CAO Lari aea eae ht ee eee, ahora. (10) 


giving, for the specimen quoted above, 7)=12-7 microseconds. The time 

constant, 7,, decreases from 7, to 0 as b/a increases from 0 to 1; To, given by 

equation (10), will be referred to as the ‘time constant’ of the discontinuity. 
The rate of change of flux is given by 


db by 5 Soluble) 5.5 ¢_y,24), ; 
ais dt >» 27pca” - J (A,) exp ( n €) St) eeeheLexerere ( ) 
with Pig On mee Me On ys aan a. ats: - (11a) 


where 7 =1:047pa2c. (This can be compared with an exponential decay curve 


which would give rt =7.) 
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The form of equation (11) is shown in Figure 4 for different values of b/a, with 
E=2rpace/b). The curve of dd/dt, obtained from Wwedensky’s expression 
is also shown. These curves give the shape of the voltage pulses that would 
be induced in a search coil, wound closely round the specimen and introducing 
no distortion itself. It can be seen that the shape and amplitude of the pulse 


1-0 


Figure 3. Decay of induction following a discontinuity. 
Full line ; equation (9) for various values of b/a. 
Broken line : equation (3) (Wwedensky). 


Figure 4. E-& curves. 


Full line : equation (11) for various values of b/a. 
Broken line : equation (3) (Wwedensky). 


varies considerably with b/a, i.e. with the distance-of the discontinuity from the 
axis of the specimen. Since yw varies with the magnetic state of the material 
the form of the voltage pulse will also change a8 a specimen is taken round : 
hysteresis cycle. 

The assumption (a) of § 2 (ii) concerning the rapidity of the original boundary 
movement can now be justified. Becker (1938) shows that for a small cubical 
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domain, at the centre of a specimen, the discontinuous boundary movement 
takes place in a time which can be written 


0-321 ha 
ey aye eee Bi a ee | 6 ea ee 


where / is the length of the cube edge, and J, is the saturation intensity of 
magnetization. 

In the case under consideration, the ‘domain’ is the region of the local 
discontinuous boundary movement. Comparing the above expression with 
equation (10), it can be seen that with J,~100 pu, provided h<a/10, the condition (b) 
will be satisfied. Although equation (12) does not hold for boundary movements 
taking place at a distance from the axis of the specimen, tp would then be reduced 
in a somewhat similar manner to 7,, and the condition 7,7, would still hold. 


§2 (ii). Form of the Voltage Pulses induced in a Search Coil 
The change in flux through a search coil following a change in magnetic 
moment at the centre of the coil is given by 


27 AMN 
Ne | B@ae= amare wen el cde (13) 


where AM is the change in magnetic moment resolved in a direction parallel to 
the axis of the coil, N, is the total number of turns in the coil, 7 is the radius of 
the coil, and 2/, the length of the coil. 

If this is compared with equation (1) it can be seen that since the induced 
voltage is a function of ¢,, the sensitivity of the system is determined in the first 
place by the dimensions of the search coil, which are controlled by the 
experimental requirements. Ideally the shape and amplitude of the pulses 
would be given by equation (11), but in practice modifications are introduced 
by the search coil circuit, and it is necessary to consider this effect in some detail. 

In examining the effect of the search coil circuit on the form of the voltage 
pulses, it is convenient to consider the decay of the flux through the specimen as 
exponential, with a time constant given by equation (10) with the numerical 
factor appropriate to the position of the discontinuity in the specimen. The 
specimen can then be replaced by a single-turn coil of similar dimensions, 
inductance L,, cross section area A,, and circuit resistance R,, giving 7, =L,/R,. 
A change in current z) through the coil is equivalent to a change in magnetic 
moment given by 

INAS nn ea (14) 

The specimen and search coil can now be replaced by two coupled circuits 
of mutual inductance M, with the ‘equivalent coil’ as the primary and the search 
coil circuit as the secondary (72=L,/R,). From the theory of coupled circuits 
(Jeans 1946), the current induced in the secondary circuit by a change in current tg 
in the primary is 


ig Ge 
a he FE a pects) 
h wo eee B= {(t1 = 72)? + 4h ry Te} 
pues oe i 27,79(1 —k?) é 27,7.(1 —&?) : 


and k=M/,/(L,L,), is the coefficient of coupling. 
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It is sufficient to examine the form of the induced voltage pulse for two 
special cases which are of practical significance: 

(a) When the time constant of the search coil circuit 1s much smaller than that 
of the discontinuity, i.e. 7, >72, the solution becomes 
. My 


i) 


_ fexp(—t/7,)—exp(—t/7e')}, sever (16) 


where 7, =7,(1 —?). 

The effect of the search coil circuit on the shape of the pulse is then confined 
to the leading edge as indicated in Figure 5. With k=1, i.e. the coefficient of 
coupling a maximum, 75’ =0 and the search coil has no effect on the shape of the 
pulse, provided 7, >75. 


a T/T, =10 
b eye 
B 7/7,’ = 40 


| 


=e 
0 05 1-0 iS 
YT, 


Figure 5. Effect of search coil circuit on form of induced voltage pulse. The numbers on the 
curves give the values of 7,/72’, where 7, is the time constant of the discontinuity, and 7,’ the 
modified time constant of the search coil circuit (see equation (16)). 


If the search coil is solenoidal, then the amplitude of the pulse is 


a _ Mi, _ 4aN,d, AM _27N,AM 
Cmax 42, maxtts = T or 7,(2l,) Ag ae are de ASG Yo (17) 


which can be compared with equation (13) with r=0. Obviously if the coil is 
not solenoidal equation (17) becomes 


27N,AM 
emax = e+ el ee eee. bed Bek Sa (17 a) 


(b) When the time constant of the search coil circuit is much longer than that 
of the discontinuity, i.e. r;<z,, the shape of the pulse is largely determined by 
the constants of the search coil circuit, and is practically independent of the rate 
of decay of the induction in the specimen. The case of practical interest is that 
for a solenoidal coil when (cf. equation (17)) 


ore 27N,AM 18 
max Tol NC ee ete ( ) 
and for a non-solenoidal coil 
, _ _2aN,AM * 
max —_ To(Ip2 ak 72) : bongs ees ( a) 


When k =1 the specimen has no effect on the shape of the pulse. 
The above two cases are of particular interest in experimental work. The 
first, (a), gives the condition which must be fulfilled if the voltage pulse induced 
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in the search coils is to be that given by equation (11). The shape of the pulse 
is then a guide to the manner in which the induction changes in the specimen 
take place. ‘The second case, (b), is that in which the conversion factor from 
pulse amplitude e,,,, to flux change in the specimen is independent of the 
magnetic state of the material and of the position of the discontinuity relative 
to the axis of the specimen. The fact that the decay of induction is not a simple 
exponential but is in the form of a sum of exponentials will obviously have little 
effect on the conclusions drawn, provided the conditions given are fulfilled 
(cf. equation (10) and (11 a)). 


§2 (iv). Signal-to-Noise Ratio 


The voltage generated in the search coil circuit by the thermal agitation of the 
electrons in the conductors sets a lower limit to the size of discontinuities which 
can be detected. Ifa valve amplifier is used in connection with the search coil, 
there is also the effect of shot noise in the first valve. If the impedance is mainly 
resistive the total noise voltage produced in the first stage can be written 


2 =4RT(R+R,) do =4RT(R+R)(v,— 1), vee ee (19) 


where Rk is Boltzmann’s constant, T is the absolute temperature, v, and v, the 
upper and lower frequency limits of the apparatus, R the effective resistance 
of the search coil circuit between grid and cathode, and R, the ‘equivalent noise 
resistance’ of the valve. If the inductance is appreciable the effect is to increase 
the noise as the frequency increases. The term R,, which may be written 
R, =2-5/g,, for a triode, where g,, is the mutual conductance of the valve, allows 
for shot effect (Terman 1943). It is thus desirable that a valve with a high value 
of g,, should be used, at least in the first stage of amplification. ‘The noise 
generated in the later stages can usually be neglected. 

It is generally assumed that the lower the effective resistance of the conductors 
in a circuit, the lower the thermal noise, but it can be seen from the above that 
little useful purpose is served by reducing the resistance of the search coil circuit 
much below the ‘shot resistance’ of the first valve. A further important point 
is that if the impedance of the circuit is mainly resistive, and if, as is usually the 
case, v, > vs, it is the upper frequency limit of the apparatus which is of importance 
in controlling the noise, which can then be written 


Ea TARE Gea Me Cee Pie a rae (19 a) 


If a search coil of length 2/, and inner radius 7, is wound with wire, m turns per 
cm. (i.e. diameter 1/n), to a winding depth of ¢, the total number of turns is 
N =2ltn?, and, if t is not too large, Ro lin*(2r +t). The signal-to-noise ratio 
is then, from equations (12) and (19), (for a given time constant and frequency 
range) : 


€m ax 


pts loaner | ae (20) 


where 7 is a mean radius. It can be seen therefore that the thickness of the 
winding should be as great as possible up to the point where the effect of an 
increase in ¢ is nullified by corresponding increases in radius and inductance. 
The dimensions of the coil, r and /, should be as small as possible. 


50-2 
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§2(v). Effect of Distributed Capacity in the Search Coil Circuit 
In practice it is impossible to wind a search coil without introducing stray 
capacity into the circuit, represented by C, in Figure 6. The input capacity 
of the following amplifier is represented by C,. Because of the presence of these 
capacities, unless the resistance of the search coil R, is sufficiently high, a voltage 
pulse induced in the coil, inductance L, will produce damped oscillations with 
a frequency 


Lees ep Cian ee 
=(a-ap| / a, or, if C is small, Poa (hale 
where C=C,+C,. The condition for critical damping is 


4L 


i 
Re C.+C,. 


If the capacities are small enough the frequency may lie outside the range 
of the apparatus and the effect may then be small. Should this not be so, and if 


Figure 6. Search coil circuit. 


R, is not large enough, further damping must be provided, ‘This can be done 
either by means of a resistance R, in parallel with the stray capacities or by means 
of aseries resistance suchas R,. ‘The required values of R, and R, for a capacity C 
are 


R,>2V/(L/C) 


Li€ 


AS Roa)’ 


or, if Cas small, Re (Cee (22) 

However since the search coil capacity C, cannot be physically separated from the 
coil, only a parallel resistance can be used in that case. Whether a series or parallel 
resistance is used with C, depends on its effect on the sensitivity of the system. 


§3. MEASUREMENTS ON THE TIME CONSTANTS OF THE 
DISCONTINUITIES 


$3 (i). Introductory 


The investigation of the accuracy of equations (9)—(11) is based on an examina- 
tion of the ‘ width’ of the voltage pulses induced ina search coil by the discontinu- 
ities, which, if no distortionis produced by the apparatus, will be of the form shown 
in Figure 4. The ‘reduced width’ €, can be measured on any convenient 
reference line and that shown as B (corresponding to ¢~0:2¢,) in Figure 4 is 
suitable. The ‘wider’ voltage pulses are produced on the average by those 
discontinuities lying nearer the axis of the specimen. Now the number of 
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discontinuities, N, with a given mean cross-sectional area, lying within a cylinder 
radius 6 (the radius of the specimen being a), is given by 


N Bw 
eee hse) age: (23) 


where JN) is the total number of discontinuities in the specimen. ‘These dis- 
continuities, N, will produce voltage pulses whose reduced width €, is greater 
than that given by the (£, €,) curve, for the given value of b/a. In Figure 7, the 
variation of N/N, with the corresponding value of &, is shown. 

Unfortunately the detection and reproduction of the form of these voltage 
pulses results in practice in a low sensitivity (low signal-to-noise ratio) which 
severely limits the accuracy of the investigation. The low sensitivity is a result 
of the requirements (1) that the search coil time constant shall be short (see 
§ 2 (iii), assumption (a), and (2) that the pass band of the amplifier shall be wide 
enough to permit the recording of pulses over a frequency range at least as wide 


Figure 7. (N/No, €p) curve. 
8 16s.w.c. hard drawn iron x 2:6. 
e 16s.w.G. annealed iron x 3:0. 
@ 26s.w.G. large grained iron x 3:3. 


as that indicated by Figure 7. Under such conditions the self capacity of the 
coil must be kept very small thus limiting the number of turns with which the 
coil can be wound. Since the self capacity of the coil cannot be estimated with 
any accuracy, the most satisfactory method is to examine the effects of various 
search coils on induced voltage pulses of known form; the required damping 
resistances and resultant resolving times can then be found by systematic trial, 
and the optimum number of turns estimated. 


§3 (ii). Method 


A block diagram of the apparatus is shown in Figure 1. The search coil, 
surrounding the wire specimen, was wound with 5,000 turns of 44 s.w.G. copper 
wire, and was 0-6 cm. mean diameter and 2 cm. long. ‘The magnetic field was 
applied to the specimen by a magnetizing coil, surrounding both search coil 
and specimen, and could be slowly changed by means of a circuit which slowly 
altered the current through the coil. The circuit is an improved version of that 
referred to in I, p. 371. 
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The resistance-capacity coupled amplifier used in this and in the counting 
experiments, was constructed from a design of A.E.R.E. (Gillespie 1947). 
The main features are its high stability, low noise and wide frequency band. 
The first is obtained by the use of negative feedback in all three circuits, each 
of the three stages of which the amplifier is made up, and by the use of high 
stability components. The first valve is a miniature triode CV138 (g,, = 8-5 ma/V.) 
giving a low shot resistance (~300 ohms), and the remainder except for the output 
stage, are CV139 miniature pentodes (g,,=7-6 ma/v.). The use of these valves 
gives a low input capacity at each stage of amplification and this, together with 
the careful choice of components has given a maximum frequency range of 
1 kc/s. to 5 Mc/s. with a maximum gain of 200,000. In recording the voltage 
pulses the band width was set at 1 kc/s. to 1 Mc/s. 

The output pulses were displayed on a cathode-ray oscilloscope, and records 
taken, with the time base switched off, on a Cossor continuous film camera, 
in which the motion of the film was at right angles to the direction of the voltage 
deflection. Film speeds of 800 cm/second were obtained with this camera. A 
timing trace was provided by the output of a square pulse generator. 

The resolving time of the apparatus (6 microseconds) was found by producing 
regular rapid changes in flux in the search coil and examining the output pulses 
on the oscilloscope. The output of the square pulse generator passed through 
a suitable coil placed within the search coil was found to be satisfactory for this 
purpose. 

§3 (11). Results 

Specimens were chosen with values of 7, lying within the range of the apparatus, 
i.e. Tt) of the order of 100 useconds, and details of these are given in Table 1. 
The records were taken near H =H, and the values of the reversible permeability 
p. are given for that value of H. 


Table 1 
Properties of the Specimens 
Diameter ws 
fonts Drak (10-* sec.) Composition (%) 
16 s.w.c. hard drawn iron 0-162 83 177 C 0:03, Mn 0:05, Fe 99-9 
16 s.w.c. annealed iron 0-162 120 250 C 0:03, Mn 0:05, Fe 99-9 
20 s.w.G. large grainediron 0-046 760 123 Decarburized specimen of 


above (see § 4 (iii)) 


Unfortunately there is no means of locating the reference line mentioned 
in §3(i) on the experimental records of the voltage pulses. The widths of the 
pulses have therefore been measured at the points where they cut the oscilloscope 
base line, which, because of the effect of stray capacities, occur at finite values of ¢. 

The measurement of the pulse width on the oscilloscope base line is not 
entirely satisfactory as it would be expected that the size of a discontinuity (as 
distinct from its position in the specimen) would then affect the measured pulse 
width; further, because of the relatively high noise level, only the larger voltage 
pulses can be detected. An attempt has been made to reduce the effect of 
these errors by carrying out measurements on as many recorded pulses as 
possible, 221 in all. 

The results of these measurements on the specimens are shown in Figure 7, 
with the measured widths multiplied by an arbitrary factor. As can be seen 
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there is reasonable agreement between the estimated and measured pulse width 
distribution. The discrepancy for high values of &, is possibly caused by the 
presence of two or more discontinuities so close together as to be indistinguishable ; 
the discrepancy in the lower values is probably caused by the finite resolving time 
of the apparatus. 

The factor by which the measured results must be multiplied may appear 
to be rather large, even allowing for the errors mentioned above: this is to be 
expected, however, since the reversible permeability used in these calculations is 
that for zero frequency and it is probable that the effective permeability for a 
rapid discontinuity is much lower. 


$4. MEASUREMENTS ON THE CONTRIBUTION OF THE 
DISCON TENUOUSSRROCES SES OME LOLA, CHANGE SUN 
MAGNETIZATION 


$4 (i). Design of the Apparatus 


In this work the problems are to detect the maximum possible number of the 
discontinuities that occur over a given part of the magnetization curve, to measure 
their number and size, and to estimate their contribution to the total change in 
magnetization. The first requirement of the apparatus is therefore a high 
signal-to-noise ratio: the shape of the amplified voltage pulse is of much less 
importance since all that is required is that the amplitude of the output pulse 
shall bear some known linear relation to the change in magnetic moment in the 
specimen such as is expressed, for example, by equation (17a) and (18a). From 
these equations it would seem that a low value of either 7, or 72, respectively, 
would give a high signal voltage, but this is offset by two other effects. Firstly 
the noise is then high, as the amplifier frequency bandwidth must have a finite 
value, which in practice tends to increase at higher frequencies. Secondly the 
effect of coil capacity becomes greater at higher frequencies. If the frequency 
of the amplifier and search coil circuit can be kept low, the frequency of the 
oscillatory component of the induced pulses may lie above the range of the 
apparatus and the resistance damping required might be appreciably reduced. 

The use of a suitably low frequency has a further advantage: this is that 
the condition 7,<r, (i.e. that the time constant of the search coil circuit and 
amplifier is long compared to that of a discontinuity) can be more easily fulfilled 
for many specimens. ‘The amplitude of the output pulse for a given size of 
discontinuity is then independent of the eddy current effects in the specimen. 
If this condition is satisfied for a discontinuity at the axis of a specimen and at 
that part of the hysteresis cycle where yu is a maximum (near 1 =H,), then the 
condition tr, <r, will hold for all discontinuities in that specimen. 

In designing the search coil its dimensions are controlled by those of the 
specimen which should have a reasonably small demagnetizing factor. It is 
desirable that the sensitivity of the coil shall be the same for all discontinuities 
taking place in the specimen and this can be attained by using a solenoidal coil 
slightly longer than the specimen. Equations (13) and (18a), which hold 
strictly only for discontinuities at the centre of a coil, will then be reasonably 
accurate for all discontinuities in the specimen. The optimum number of 
turns is best found by trial as described in § 3 (ii). 
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4 (1). Method 

The basis of the experimental procedure, method of calibration and necessary 
precautions are similar to those described in paper I. 

A block diagram of the apparatus is shown in Figure 1. The method of 
magnetizing the specimen, which is in the form of a wire 7:5 cm. long, is the same 
as that referred to in $3 (ii). The search coil which was finally found to be most 
satisfactory was 10 cm. long with a mean diameter of 0-9 cm., wound with 52,000 
turns of 47 s.w.G. copper wire. ‘The amplifier was the same as that described 
in §3(i) with the frequency response set to limits of 1 kc/s. to 3 ke/s. ‘The 
amplified voltage pulses were fed through a discriminator to a scale of 1024 (21°), 
and a recorder of the post office type. The discriminator was used with its 
upper limit of discrimination set at 25 volts, which could be reduced in 
subdivisions of 0-1 volt. 

The search coil required a parallel damping resistance of 200,000 ohms, 
giving a resolving time for the coil and amplifier of 300 wsec., while the resolving 
time of the counting system was 100 psec. The noise level at the input of the 
amplifier with this assembly was about 1:5 wv. (i.e. 0-3v r.m.s. output). 

In the calibration of the apparatus, the specimen was replaced by a small 
coil of 10 turns, 0-165 cm. mean diameter, through which current pulses of 
square waveform were passed. The rapid changes in current Az in the 
calibrating coil simulate the rapid changes in flux in the specimen during a 
discontinuity, the relation being given by 


gSAt 
10 by eeeoes 


with g the number of turns in the calibrating coil and S the area of cross section 
of the calibrating coil. The discriminator was set, for each value of Ai, so that 
the amplified voltage pulses were just transmitted: a calibration curve was thus 
obtained showing the relation between discriminator setting and change in 
magnetic moment. ‘The assumption that the time constant of the search coil 
was greater than that of the change of flux in the calibrating coil was confirmed 
by examining the respective waveforms of the input current and output voltage 
on the oscilloscope. 

It should perhaps be emphasized here that although the equations (13) and 
(18) are not used directly in the calibration of the apparatus, they are needed 
to show that there is a constant ratio between change in magnetic moment and 
the amplitude of the induced voltage pulse (discriminator setting) over the 
frequency range used. 

‘The experimental procedure was as follows. The specimen was magnetized 
to near saturation in one direction, the field reversed several times to bring the 
material into a cyclic condition, and the specimen then taken through one half 
cycle of a hysteresis loop. With the discriminator set at a chosen level, readings 
of the number of discontinuities were taken at suitable values of the field. The 
rate of change of field was adjusted for each sample, insome cases so that a complete 
half cycle took up to 24 hours; the readings were repeated at each level at least 
three times and a mean taken. The lowest discriminator setting corresponding, 
in terms of AM, to AM=0-3 x 10-® E.m.u. was above the peak noise level of 
the apparatus so that there was no background count. This can be compared 


with the sensitivity of the apparatus in paper I, where with AM =3-0 x 10-° E.m.v. 
there was a noticeable background count. 


AVE 
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Results. 
Experiments were carried out on the specimens specified 1n ‘Tables 2 and 3. 


Table 2. Composition ot Specimens 


Diamet sie 
Cok Composition (%) 
Hard drawn iron (Fel) 0-384 C70- 03 Minw OZ 05,eke 99-9 
Large grained iron (Fe2) 0-378 Decarburized specimen of above 
Annealed nickel (Ni) 0-453 CxO 04551-0055 ben0-15. 


Mn 0:03, Al 0:01, Mg 0:13, 
S 0:003, Ni balance 


Table 3. Magnetic Properties of Specimens 


T 1s calculated from equation (10), using the maximum values of reversible 
susceptibility « (permeability w) given. 


: Maximum k To Average grain 
Matenal Is (#4) ie (10-8 sec.) size (cm?) 
Hard drawn iron 1720 6-9 5:7 10-1 Sp LO 
(87-7) 
Large grained iron 1720 32 0-39 46-0 Ome 
(403) 
Annealed nickel 485 8-0 163 16°6 case 107° 
(102) core 10~° 
1200 Se rane Gear | ae 
Fe 2 
Fel 
800 ;— 
FERS Ni a 
sty oe 
400+ | 
800 [- = 
is a 
(GSS Le ie rene Sar ese (ees oN ce Sys ea ee eee | ees | 
25 20 15 10 5 0 5 10 15 20 25 


Figure 8. Hysteresis curves for hard drawn iron (Fe 1), large grained iron (Fe 2), 
and annealed nickel (Ni). 

The large-grained iron was prepared by the method of Edwards and Pfeil 
(1924): a specimen of hard drawn iron was decarburized in an atmosphere of 
moist hydrogen for 48 hours at 950° c., stretched by 3:5°%, then annealed in dry 
hydrogen at 830°c. for 72 hours. The nickel was annealed in hydrogen 
at 930°c. for 1 hour and then allowed to cool over a period of 6 hours; 
this treatment is similar to that used by Bates and his collaborators. The average 
grain sizes shown in Table 3 were obtained from measurements on etched specimens. 
The hysteresis loops for the three specimens are shown in Figure 8, for which 
I is plotted against the applied field. For the present purpose a knowledge of 
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Figure 9. Number-level curves for iron and nickel, showing the number of discontinuities 
per cm? with a magnetic moment greater than AM, 
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Figure 10. Distribution curves for iron and nickel. (The minimum volume associated with the 


change in resolved magnetic moment AM can be obtained by multiplying by 2-91 4 
for iton and by 1-03 x 10- for nickel.) Mens eae ie AEE 
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the true field (i.e. corrected for demagnetization effects) is not required directly. 
The corrected curves do not differ greatly from those shown except for the large 
grained iron specimen. 

From the readings taken of the number of discontinuities at each level of the 
discriminator, a curve was drawn for each specimen showing the number of 
discontinuities per cm’, above a given size AM, plotted against the corresponding 
discriminator setting (in terms of AM). Number-—level curves for the three 
specimens are shown in Figure 9. By differentiating these curves, a (size, 
number) distribution curve was obtained for each specimen, i.e. (AM, dN/dAM). 
‘These are shown in Figure 10. 

Figure 11 shows the contribution per unit range of size, AM.dN/dAM, to the 
total change in magnetization, A>, plotted against AM: for the purpose of 
comparison the ordinates have been divided by the appropriate value of AI». 
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Figure 11. Contribution curves for iron and nickel showing the relative contributions to the total 
change in magnetization per unit range of AM. 


‘The contribution of the measured discontinuities, over the measured range of 
AM, is given by the area under the contribution curve 


F dN ; 
| Amy daM) = | AM.dN, 

with appropriate limits, which is thus the corresponding area between the 

number-—level curve and the N axis. This is the measured change in intensity 

Alp, produced by the Barkhausen effect, over the field limits considered. ‘The 

results of this integration carried out directly on the number-level curves are 

shown for each specimen in Table 4. 

The relative contribution of discontinuities, above a given size, to the total 
change in magnetization AJ,/AJ; is shown in Figure 15, which has been 
obtained by integrating the area under the contribution curve above chosen 
values of AM. 

An estimate has also been made of the contribution of reversible processes 
to the change in magnetization, by the measurement of the reversible 
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susceptibility, «, at different points on the hysteresis curve. (It is hoped to 
publish details of this work later.) A (x, H) curve can then be drawn and the area 
under this curve is given by 


He 
| «d= len dH=Al, 


with appropriate limits, which can be referred to as the reversible contribution. 
This is also shown in Table 4, and the sum of the two contributions AJ, +ATJ, is 
compared with the change in intensity taken from the hysteresis loop AJ;. The 
value of AJ, +AZ, obtained by extrapolating AJ, to zero AM is shown for the 
iron specimens, but this cannot be done satisfactorily for nickel as the 
contribution curve shows no maximum. 


Table 4. Contribution of Discontinuous and Reversible Processes (AJ, and AJR) 
to the Total Change in Magnetization AI, 


Limits 


Limits 


Contribution (and percentage 


Material aE ofAM Ale (49- a ie ofa) 
field (10-*:.m.Uv.) : Al, AI, AIk+AR AR+AR 
extrapolated 

Hard drawn +18°6 0-3-35 2150 1-20 1844 180 2024 2109 

iron (86) (8:4) (94) (98) 
Large + 3-0 0-3-20 2180 0-2-4 12505 110 1400 1390 

grained iron (57) (4-9) (62) (64) 
Annealed +14-9 0-3-12 700 1-20 S20meelalS 441 

nickel (AZ) 6-4) (63) 


The minimum values of wv (v,,) (from equation (16)) corresponding to 
AM x 10®= 1 are, for hard drawn and large grained iron (J, = 1720), 2:9 x 10-10 cm3, 
and for nickel, (J,=485) 1-0; x 10-® cm’. The corresponding number of atoms 
can be readily found. The minimum values which could be detected with 
certainty above the noise level in the experiments, also the maximum values. 
(as given in the results), are shown in Table 5. 


Table 5. Minimum and Maximum Changes in Resolved Moment AM, and 
the Associated Minimum Volumes v,, and Numbers of Atoms n, 


Fe 1 Fe 2 Nickel 
Min. Max. Min. Max. Min. Max. 
AM x 108 0-30 40 0:30 20 0-3 12 
Brey SC MOO 0:87 116 0:87 58 3-1 124 
Lie NOR 0-74 99 0:74 49 2:6 106 


An ‘average value’ of AM (or v,,) can be calculated for each specimen by 
dividing the total measured contribution by the total number of measured 
discontinuities but, as was pointed out in I, p. 379, this value is meaningless 
because of the probability of a large number of much smaller discontinuities 
which contribute very little to the total change in magnetization. For purposes 
of comparison, however, an ‘average value’ of AM (or v,,) can be defined as the 
value above which the discontinuities account for half of the total discontinuous 
magnetization (i.e. }(AJ,—AJ,)). This average value, which has been calculated 
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from Figure 15, together with the corresponding values of v, and nm, are shown 


in Table 6. 
Table 6., Average Values of AM, v,, and n, 


AM Om Nyy 

(10-§ E.m.v.) (10-1° cm?) (101 atoms) 
Hard drawn iron 1:6 4-6 (55) 
Large grained iron 0:7 2:0 2°8 
Annealed nickel 0-4 4-4 4:8 


In order to demonstrate the variation in contribution with field, curves are 
shown giving dJ,/dH plotted against H, compared with the corresponding 
curves for dl,/dH obtained from the hysteresis loop, and with dJ,/dH, the 
reversible susceptibility (Figures 12-14). 
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100 


Figure 12. dJ/dH curves for hard[drawn iron. 
—— dIy/dH. —-—dI3/dH. ---- dIp/dH x10. 


Figure 13. dI/dH curves for large grained iron. 
— dIn/dH. —e— dIp/dH. ---- dIp/dH x10. 
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The measured sum AJ, +AZ, for hard drawn iron is 94% of AZ, and 
extrapolation of the contribution curve from AM=0-3x 107° E.M.U. to zero 
accounts for another 4°%, giving a total of 98%. Since the probable error in 
the measured Barkhausen contribution is about 5%, and the extrapolation can 
give only the order of magnitude of the contribution of the smaller discontinuities, 
the agreement is partly fortuitous and there is a possibility that there is still some 
further unaccounted contribution for the specimen (see also below). 
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Figure 14. dJ/dH curves for annealed nickel. 
—— dIy/dH. —-— dlIz/dH. - --- dIp/dH x 10. 
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Figure 15. Fractional contribution AJp/ AJ» of discontinuities above a given size AM. 


The results for the large grain iron show that 38% of the total magnetization 
remains unaccounted for by discontinuous processes. Extrapolation below 
AM=0-3 x 10-® E.m.u. would indicate that these smaller discontinuities can 
contribute to only a small proportion of this discrepancy. One possible source 
of error lies in the fact that the Barkhausen discontinuities do not take place at 
random, but tend to occur in groups, i.e. the production of one discontinuity 
tends to trigger a number of others. This is the main reason for the extremely 
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low rates of change of field (10-4 oersted/sec.) required in order to keep the rate 
of counting within the range of the apparatus and to reduce the number of 
coincidences between the discontinuities. In the experiments on the hard 
drawn iron and the annealed nickel the effect was small, but with the large grained 
iron the presence of the ‘ grouping’ was still noticeable and this may be the reason 
for the sharp fall in contribution for low values of AM. The effect may be associated 
with the steepness of the (J, H) curve over a wide range of J, which means that 
at any particular (J, ) value in this range, a very small increase in H is sufficient 
to induce a relatively large change in J. The local change in field due to a single 
discontinuity may therefore be sufficient to initiate others. It is possible that 
the effect would be reduced by using a specimen with a higher demagnetizing 
factor, and experiments to test this are projected. Should this explanation be 
correct, the contribution of the discontinuities below the level AM=1 x 10-6 cm. 
is greater than is indicated in the results. 

The lower contribution for annealed nickel (47°) is quite definite but shows 
a marked increase over that given in paper I, 11° for a hard drawn specimen. 
The percentage contribution over the steep part of the hysteresis curve appears 
to be somewhat greater, but it is difficult to estimate this accurately because of 
the rapid change in intensity over this region. Because of the absence of a 
maximum in the contribution curve for this specimen, it was not possible to 
make an accurate extrapolation over the region below AM=0-3 x 10-® E.m.u., 
but an examination of the curve indicates that the smaller discontinuities could 
account for a considerable part of the discrepancy. 

The results for the ‘average size’ of the discontinuities show that, for the 
hard drawn iron, half of the estimated total discontinuous change is accounted 
for by discontinuities for which AM>1-6 x 10-* E.m.u. This value is more than 
twice that for the large grained specimen in which presumably flaws and 
inclusions are smaller. This supports the view that the larger Barkhausen 
discontinuities may be associated with the secondary domain structure round the 
imperfections in the material, as suggested by Williams and Shockley. Where 
the difference in the values of J, is taken into account in estimating the average 
volume, there is no significant difference in the values for annealed nickel as 
compared with the other specimens. ‘The order of magnitude of these average 
minimum volumes is about one-tenth those given in paper I, and by Bozorth 
and Dillinger, and this is a measure of the improvement in the sensitivity of 
the method. . 

The dI/dH curves are of interest in showing a maximum which, for a given 
specimen, occurs at nearly the same value of H for the Barkhausen, reversible, 
and total contributions. For the hard drawn iron it seems that the Barkhausen 
contribution may be an overestimate, where d/,/dH is greater than d/,/dH, but 
the difference is considerably greater than the experimental error, particularly 
as in this region the rate of counting is not very great. It is more probable that 
this is caused by magnetic viscosity, i.e. by the difference in the rates at which 
the specimen is magnetized in obtaining the hysteresis curve and in the counting 
experiments. The dJ,/dH and dl,/dH curves for the large grained iron both 
exhibit a well defined region of constant d//dH, corresponding to the steep 
part of the loop: the noticeable curvature in the dI,/dH curve may perhaps 
indicate that the method of measuring reversible susceptibilitv is more sensitive 
in detecting small changes in slope on the (J, H) curve. 
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There is one possible source of error in any work of this nature which there 
seems to be no means of avoiding, except by increased sensitivity. This is that 
if there are a large number of small discontinuities below the noise level, these 
will raise the level of the noise as soon as the field is changed and discontinuities 
take place. The result would of course lead to an increase in the measured 
number of discontinuities. ‘There was, however, no obvious increase in the noise 
as examined on the oscilloscope, although it would be difficult to differentiate 
between genuine discontinuities and an increase in the background noise, nor 
is there any indication, except for the hard drawn iron, Figure 12, that the 
measured AJ, is too large. It is unlikely that the steep rise in the number-level 
curve at low values of AM is caused by noise; as can be seen by comparison 
with the inset curves, the rise is just as pronounced at higher values of AM if 
the scale is suitably chosen. 

It has been considered desirable in this work to keep the demagnetizing 
factor for a specimen low in order to keep the internal field uniform and so that the 
fields required to saturate the material would be reasonably low. Specimens with 
the above dimensions (7:5 cm. long and 0-376 mm. diameter) have therefore been 
used. This has resulted in the large number of counts which have to be taken, with 
the rather long time (up to 3 months) necessary to investigate a single specimen. 
It is proposed in any future work to use a shorter specimen about 1 cm. long, 
and this should reduce considerably the time taken. Because of the higher 
demagnetizing factors and the fact that the effective slope of the (J, H) curve will 
be reduced, it will only be possible to carry out experiments over a more restricted 
part of the magnetization curve at any one time; it will, however, have the 
advantages that lower effective rates of change of field will be possible and the 
‘grouping’ of the discontinuities should be reduced. Further it will then be 
possible to use a shorter search coil with consequent reduction in the thermal 
agitation noise. 


§4. CONCLUSION 


It may be convenient to summarize some of the main results and tentative 
conclusions. 

(1) A theory of pulse shape has been developed, with which the experimental 
results are in reasonable agreement. 

(2) The sum of the measured Barkhausen and reversible contributions 
agrees quantitatively with the measured total change for hard drawn iron and the 
discrepancy for the large grained iron and annealed nickel can be explained. 

(3) 50°% of the discontinuous change in magnetization is due to individual 
discontinuities with AM>1-6 x 10-* g.m.v., 0-7 x 10-8 E.m.u., and 0-4 x 10- E.M.vu., 
for Fe 1, Fe2, and Ni, respectively. The Barkhausen volumes (~4 x 107-!° cm$) 
are in general much smaller than primary domains, at least those in single crystals, 
and may be attributed to localized irregularities in the displacement of primary 
domain boundaries. 

(4) The number of larger discontinuities is greater for hard drawn iron than 
for a decarburized and annealed specimen but the distribution curves are not 
widely different. Measurement on a wider range of materials is necessary in 
order to see whether there are differences which can be related to the other 
properties of the materials. 
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Domain Wall Movement in a Single Crystal 


By K. H. STEWART 
Cavendish Laboratory, Cambridge 


Communicated by D. Shoenberg; MS. received 12th Fanuary 1950 


ABSTRACT. A single crystal of silicon—iron was shaped so as to encourage the formation 
of large ferromagnetic domains. Its magnetization curve was found to contain a large 
Barkhausen jump, indicating reversal of magnetization in a volume of the order of 1/10 cm’. 
The rate of change of magnetization was found to be proportional to the excess of the applied 
magnetic field over a certain critical field, and could be made very slow. It is suggested that 
the rate of change is controlled by the braking effect of eddy currents on the movement of 
domain walls. 


magnetization in a single crystal of 3%, silicon iron cut into the form of a 

hollow rectangle (‘ picture frame’) all of whose faces were parallel to (100) 
planes of the lattice. Williams (Williams and Shockley 1949, Williams, Bozorth 
and Shockley 1949) has demonstrated the existence of remarkably simple domain 
arrangements in such a specimen and it was hoped that the study of magnetization 
changes with such a simple arrangement would give useful information about the 
factors controlling the movement of domain walls. . 

The present specimen had the dimensions shown in Figure 1. The orienta- 
tions of its faces were tested by etch pit reflections and by the ‘fir-tree’ Bitter 
patterns (Williams, Bozorth and Shockley 1949), and the specimen was filed 
until all the faces were within 1° of (100) over most of their surface. The 
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Toms note reports some observations that have been made recently on 
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characteristic domain structure found by Williams is indicated by the dotted lines 
in Figure 1, directions of domain magnetization being shown by arrows. ‘The 
Bitter patterns formed when colloidal magnetite was applied to the present 
crystal revealed no such definite pattern; all that could be seen was a number of 
isolated ‘fir-trees’, on parts of the specimen where the orientation was slightly 
imperfect. This examination was not pressed very far as the electrolytic polishing 
necessary before the patterns could be seen would have consumed too much of 
the specimen. It may be worth recording that it was found that ‘ fir-tree’ and 
straight line patterns, i.e. patterns characteristic of the true domain structure and 
not merely of the state of strain of the surface, could be obtained on this and other 
specimens if careful mechanical polishing was followed by an anneal at 1200°c. 
ina high vacuum (~10-4mm. Hg). If the vacuum was less good, a film of tarnish 
formed on the polished surface and prevented the patterns from appearing. 
This method has advantages over electrolytic polishing in cases where the exact 
shape of surface is important. 


Figure 1. The single crystal specimen: domain boundaries are shown by broken lines, 
directions of magnetization by arrows. 


The main interest of the specimen lay in its magnetization curve, which was 
obtained by the usual ballistic galvanometer method. The field winding consisted 
of six turns of 368.w.G. copper wire on each of the long limbs and three turns on 
each of the short ones. The flux coil was normally ten turns of 42s.w.c. wire on 
one of the limbs. It was found that the results did not depend on the exact 
arrangement of these coils, in fact in some of the later experiments the eighteen- 
turn field winding was replaced by a single thicker wire threaded through the 
specimen, without affecting the results for a given magnetomotive force. 

The magnetization curve is shown in Figure 2. Its chief feature is the large 
instability in fields of about 0-05 gauss; it can be said that almost half the total 
change in magnetization takes place in a single Barkhausen jump, occurring at a 
reproducible value of the field. Large Barkhausen jumps have been observed 
before, notably in the experiments of Forrer (1926) and Sixtus and Tonks (1932), 
but always in material subject to severe stresses and not in material where the 
magnetization directions are controlled by the magneto-crystalline anisotropy, a- 
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in the present well-annealed silicon—iron crystal. The presence of such a large 
discontinuity in the magnetization curve suggests that very large domains are 
involved and that the change in magnetization occurs chiefly by the movement, 
perpendicular to its plane, of a single 180° wall such as that shown in Figure 1. 

The ‘jump’ from A to B in Figure 2 took place quite slowly. It proved 
possible to arrest it by switching off the (negative) field as soon as movement of 
the galvanometer showed that the jump had begun. If the negative field was 
applied again gradually, the jump would re-start at a lower value of the field than 
the one corresponding to point A. The true hysteresis loop is therefore not the 
dotted curve of Figure 2, but rather the re-entrant full curve. This was at first 
plotted by manual control of the field as just described, but later by an automatic 
feedback method which adjusted the field so that there was a constant change of 


0-05 0:10 
H (oersted) 


Figure 2._ (J, H) curve of the specimen. The detail of the re-entrant part of the loop is taken 
from one ‘run’; other runs give the same main features but differ in detail. 


raagnetization, at a very slow rate, from A to B. Figure 2 can be interpreted in 
terms of the movement of the 180° wall of Figure 1 if it is assumed that the movement 
is opposed by various obstacles such as those due to local internal stresses (Becker 
and Déring 1939) or inclusions of impurity (Kersten 1943). The curve suggests 
that there are many obstacles of about the same size in this specimen, that at A 
being larger than all the others as far as B. ‘The hysteresis loss of the specimen 
will, of course, depend on the way in which the field 1s changed, being considerably 
larger if the dotted line AB is followed than if the field is adjusted to follow the 
re-entrant curve. Similar effects may well occur in more normal materials, but 
in these it would be very difficult to adjust the field to follow the small and rapid 


discontinuities. 
sie 
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The speed of the change from A to B was found to depend on the field strength. 
It was investigated by bringing the field to a point just short of A and then increasing 
it suddenly to some larger negative value; the subsequent changes in magnetization 
were measured by a photographically recording galvanometer. The results for 
the initial rate of change of magnetization are shown in Figure 3. If the change is 
interpreted as due to the movement of a single 180° wall, we obtain a linear relation 
between the field H, and the velocity of the wall V, given by V = A(H — Hy) where 
A and H, are constants equal to 6-3 cm/sec/oersted and 0-049 oersted respectively. 
The experiment was repeated with the specimen cooled in liquid oxygen and the 
lower line of Figure 3 was obtained, with A=5:7 cm/sec/oersted, Hy =0-056 
oersted. 

These values of wall velocity are much lower than those that must occur in 
ordinary polycrystalline silicon—iron, since such low values would give rise to 
viscosity effects much larger than those commonly observed. Various effects 
have been suggested as controlling the speed of wall movement (Snoek 1938, 1947) 
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Figure 3. Rate of change of flux density as a function of field strength at room temperature 
(circles) and at the temperature of liquid oxygen (crosses). 


but the only one which provides a ready explanation for the very slow movement 
in the present specimen is the control by eddy currents. Since eddy current 
effects depend markedly on dimensions, the present very large wall would be 
expected to move much more slowly than the usual small ones. A rough calcula- 
tion, assuming that the wall moves as a rigid plane and induces the same eddy 
currents as would an infinite plane wall, gives a value of A of 1:0 cm/sec/oersted 
forthe specimenofFigurel. ‘Thisis only onesixth of the value found from Figure 3, 
but the very rough assumptions of the calculation would certainly lead to an 
underestimate of A, so that the discrepancy is not serious. ‘The idea that A is 
determined by eddy current effects is supported by its variation with temperature. 
If the speed of movement were controlled by effects such as the relief of elastic 
strains, it should have a very large temperature dependence; if by eddy currents 
then the speed should be proportional simply to the resistivity of the Seca 
In fact, Figure 3 shows that A decreases by 10° when the specimen is cooled from 
room temperature to liquid oxygen temperature, while the measured decrease in 
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resistivity for a strip sample of the same material cooled in the same way was 18°, 
In view of various experimental inaccuracies this order-of-magnitude agreement 
seems to support the eddy current hypothesis. 

The linear relation between V and H shown in Figure 3 is similar to that found 
by Sixtus and Tonks (1932) in their experiments on stretched nickel-iron wires, 
but in their case the values of the constant A were much greater, of the order of 
10,000 cm/sec/oersted. ‘Their original interpretation of the results suggested 
that A was determined by eddy currents, but Snoek (1947) has shown that this is 
probably incorrect and that some additional ‘frictional’ mechanism must be 
present. In the single crystal specimen it is apparently unnecessary to postulate 
such a mechanism to account for the observed effects; if it exists it is masked by 
the much larger ‘friction’ provided by eddy currents. 

Further experiments on various phenomena in the crystal would have been 
desirable, but unfortunately it was destroyed while being annealed in an attempt 
to improve its properties. 
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Sommerfeld’s Polynomial Method Simplified 


By A. RUBINOWICZ 


Institute of Theoretical Mechanics, University of Warsaw 
MS. received 30th September 1949 


ABSTRACT. From two earlier papers by the author it follows that eigenfunctions 
obtainable by Sommerfeld’s polynomial method contain either a Riemann P-function or a 
degenerate form of it when the fundamental interval is either finite or infinite respectively. 
Comparing the differential equation for the P-function or its degenerate form with the 
differential equation which we get for P from the original differential equation of the 
eigenvalue problem, we obtain standard formulae for the constants determining the eigen- 
values and eigenfunctions. ‘The calculations necessary to give the complete solution of an 
eigenvalue problem, which can be treated by Sommerfeld’s polynomial method, are thus 
reduced to a minimum of differentiations and simple algebraic operations. It is also 
possible using appropriately chosen P-functions to obtain the eigenfunctions directly in 
the required form. The method yields also eigenfunctions for eigenvalues in the continuous 
spectrum. 


§1. INTRODUCTION 


N two earlier papers (Rubinowicz 1949a,b) it was shown that the 
| solutions of the eigenvalue problem belonging to the self-adjoint differential 


equation 
d [{_ df. 
z(Px) = (¢—Ap)f{=0 1 a a eee (1) 
and obtainable by Sommerfeld’s (1939) polynomial method are given by 
fare? 2 (2) 


The function P is here either of the form 


(1) P=E(1 (5) Nr b, Cc, é) 


or of the form 
(II), Fer exp (5/2) i (@5cGe), 
where E=Kxh, (3) 


eeeue vale, 


F and F, are the ordinary and the confluent hypergeometric functions and 
his a positive integer. We use function (I) where the fundamental interval of the 
eigenvalue problem (bounded by two singular points of the equation (1)) is finite 
and its boundary points are given by €=0 and €=1, and function (II) when this 
interval is infinite and bounded by €=0 and €= +00. But ifhis an even integer 
and «= =0 is not a singular point of the differential equation (1) the boundaries 
of the fundamental intervals will be given in the two cases (1) and (II) only by 
€=1 and = + ~ respectively. 

In this paper the above results will be used to show that we can obtain the 
eigenfunctions and eigenvalues quite simply in all cases where Sommerfeld’s 
polynomial method can be applied. ‘The same procedure can also of course be 
used independently of its connection with Sommerfeld’s polynomial method. 
= The starting point of our considerations is the remark that the function P 
is in case (I) a Riemann P-function and in case (II) obviously a confluent form of it. 
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§2. GENERAL RELATIONS 


Inserting the function f as defined by (2) into the differential equation (1) we 


get for P, as a function of the variable é defined in equation (3), the differential 
equation 


PP h-11dP 1487/1 dp\?_ 1 dp g~d 
d2 h &dé ' RE |\2p dx 2p dx? =] 


In case (I) P is a Riemann P-function and therefore a solution of the differ- 
ential equation 


ht — l-—y- Ye (= + 


apt (eee + (te my 
wets (5) 
where Oto +B +S yty al eee (5 a) 
The constants a, b,c in (I) are then given by 
a=at+B+y; b=a+ f'+y; Calera =e a, | ances (5d) 


Comparing the coefficients of dP/d€ in (4) and (5) we get (h—-1)/h=1—a-«’, 
0=1-—y-—~y’ and hence from (5 a) 


a +a’ =1/h, B+f' =—-1/h, Vey La” ae (6) 


Comparing coefficients of P in (4) and (5), multiplying by A?é? and resolving 
into partial fractions in €, we finally obtain our standard formula for case (1) 


1 dp apes 1] 
“l(s a DE: 
-ie[ - hoe " (2 Re +8(8-+1/h) |. 
cee: (7) 


In order to apply the same procedure in case (II) we remark that the confluent 
form (II) of the P-function is obtained in the limit «0 if in (I) we put 


ae B=By+1/e, v=yot1/e, 
ee nee, 12 pa ae . 
For the coefficients a’,c’ in (II) we then get 
a =a+ fo +70, CS eon Oval oR eter hee (9) 


and writing é instead of é’ we obtain for the function P in case (II) the differential 
equation 


BP | 1-a-o! dP (Gee tis 1) po, 
de UE & 2g “t 
From (4) and (10) follows 
Cee Sale eel eS A eee ere (11) 
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and therefore also the standard formula for case (II) 

indp\e> SAedsp | f ( aol 7) ¥ 4 
21 [ — ~ =h?} —a(a—I1/h ~ ~=— -~)é—-+é" |. 
® (aH) 2p dx p HE ol HEN Oe Since cme 


In deriving (12) we have used the fact that By — By’ —yo +¥o =2(a—@’) +1—I/h 
which follows from (5a), (8), (9) and (11). 

For A =1 equation (10) is identical with the differential equation for Whittaker’s 
confluent P-function W. 

When, using Sommerfeld’s polynomial method, we obtain a solution of an 
eigenvalue differential equation of the form (1), its coefficients p, g and p must be 
such that they identically satisfy equation (7) for case (I), and equation (12) for case 
(II). If we split the left-hand side of (7) or (12) into partial fractions and equate 
the two sides of the equation in question, using (3), we obtain the value of h, and 
in case (I) also that of x. Further, since in both these equations we have three 
terms on the right hand side, we obtain three equations for the six constants 
(x, h, «, B, y, 4) involved in case (1) and for the five constants (x, h, «, a’, A) in case 
(11). ‘These equations are quadratic in « and in case (I) also in B and y, being of 
the form a(a—1/h)=A, B(B+1/h)=B and y(y—1)=C. Generally we have 
therefore two values of «, 8, y at our disposal. For discrete eigenvalues these 
constants must have fixed values or values which depend upon a quantum number. 
In both our cases we still need for this purpose one equation corresponding to 
Sommerfeld’s condition for the breaking off of the power series and containing the 
quantum number. ‘This equation we get, possibly together with inequalities for 
our constants, by investigating the behaviour of the eigenfunctions at the boundary 
points of the fundamental interval. In case (I) these points coincide with 
€=0 and €=1 and in case (II) with €=0 and €=+00. The convergence of the 
normalization and orthogonality integrals 


| fitfep dx = | PirPLh ae We srs (13) 


at the point x = =0 depends, in both cases (I) and (II), upon the factor £*(p/p)¥? 
(also if «=O is aninner point). This factor may suppress one of our two «-values. 
One at least of the boundary points is given in case (1) by €=1 and in case (II) by 
= +00. In case (I) the convergence of (13) at €=1 depends upon the factor 
(1—£)” (p/p)"? if F(a, 6, c, €) converges for€=1. This factor can cause a suppres- 
sion of one of our two y-values given by y=$}+(C+4)!? and y=1— (C+4)'2, 
But if a or b is not a negative integer the hypergeometric series F (a, b, c, €) does 
not terminate and converges for €=1 only if Re(a+b—c)<0, i.e. from (5 5) and 
(6), if Re(y)<4. In every special case in which we are concerned with discrete 
eigenvalues it must be impossible to fulfil this inequality so that we must choose 


for a or ba negative integer —n(n=0,1,2....). Otherwise we would not have a 
quantum number. For a= —n we obtain according to (55) and (6) 
p=-n-a-y, Bo=n+at+y-lh  ...... (14 a) 
and for b= —n 
B=n+a+y—1/h, Bo=—n-a-y canes. (145) 


Since the constants 8 and f’ enter symmetrically in our formulae we get the same 
final results for both the cases a= —n and b= —n. 
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We must apply a similar procedure in case (II). Here the behaviour of 

oe, c’,€) for €> + © is given by the semiconvergent expansion (Sommerfeld 
F(a‘,¢',£)= mea Olas a7 ean 2) 

Le Perr Fy ey 

aves EV C'e*G(a’' — 1, c’ —a’, — 1/8), 

a(a+ zl +1) ge 

In order that the integrals (13) may converge for real values of €, we must 
therefore assume a’ = —n which gives us a discontinuous eigenvalue spectrum. 
In the case of continuous eigenvalues the variable € and hence also x must be 
imaginary. In §4 we shall show by an example how to treat such a case. 

We can use formulae (7) and (12) also for the determination of the limits of 
applicability of Sommerfeld’s polynomial method (cf. Rubinowicz 1949 a). 
Having given p and p we can determine, from these formulae, the functions q 
for which an eigenvalue problem (1) can be solved. 

Finally we remark that, instead of the special functions (I) and (II) suggested 
by Sommerfeld’s original method, we can also take for P other solutions of the 
differential equations (5) or (10). Such a procedure of course cannot change the 
equations for our constants but it will in general require us to take other values of 
the solutions of these equations in order that the integrals (13) may be convergent. 
It will also be necessary to replace the equations (14) or the equation a’ = —n by 
other equations representing Sommerfeld’s condition for the breaking off for the 
special P-functions we use. 

It is clear that for non-degenerate eigenvalues we obtain, by this method, the 
same eigenfunction for a given eigenvalue as in the case of the P-functions (I) or 
(II), since according to our assumptions these functions are uniquely determined. 
The advantage of using P-functions other than (I) or (II) is that we can obtain 
the eigenfunctions in different forms containing, for example, a polynomial in 
descending powers of €. 


where G(a,b,£)=14 ge 


$3. ASSOCIATED LEGENDRE FUNCTIONS 
An example of an eigenvalue problem in a finite fundamental interval (case (I)) 
is that for the associated Legendre functions. Its differential equation is given by 


(0-08) - (ya 


where m is an integer. We have therefore here the differential equation (1) with 
p=H=l-x*, p=, Geil Raa ae ls Sachse (15) 
The boundary points of the fundamental interval of this eigenvalue problem 
are the singular points of our differential equation situated at w= + gb 
If we put (15) into equation (7) and split into partial fractions in the variable x, 
we get 
ie a ein es Oa d 
ae x?)2 1— x? 
y(y—1) , —a(a—1/h) + B(B +1/h)+yy—1) _ ] 
=) ae oa p(e+1/n) |. 
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Hence using (3) we have =x? so that 
[aes KeLS Se ee tte (16) 
and therefore 4y(y —1) =m? —1, 48(8 +4) =A, 
4{—a(a—}) + B(B +4) +r(y—D) =m? tA=1. 


From these equations it follows that 


a(a—4)=O0sothatxa=Oora=}, = = 9 see (17 a) 
A=Kl+ lat we put 28 =h | eee (17 5) 
vy=(n + 1)/2 fonm=0, = ee eee re (17 c) 


because y=(—m-+1)/2 with m>0 would yield a function f with singular points 
at x= +1 and hence an infinite normalization integral (13). 

Since, according to (17 a) and (17c), « and y are integers or half odd integers 
and n is an arbitrary integer, the numbers /=2f are, according to (145), arbitrary 
integers. Hence the eigenvalues A=/((/+1) are determined. 

In order to get the eigenfunctions we must calculate the constants a, 8, c. 
From (56), (16) and (17) it follows that 


a=(l+m+1)/2+a, b= —(l—m)/2 +a, c=2e +4. 
According to (2), (I) and (15) we obtain therefore the eigenfunctions 


lim+1 l—m 
2 


fie (ae) cee +a, — —s— +a, 2«+4, «*) ees (18) 


2 
where « is equal to 0 or 4 according as (J—m) is an even or odd integer. We 
obtain (175) and (18) also, if we put 28 = —/—1. 

Equation (18) represents the associated Legendre functions apart from 
a constant factor. We can see this by using the relation 


F(a, b,c,z)= Tirana (na Fe l=e4 alba) 
Trea ao 1 = hiiea ah ae (19) 


which, when a or 6 is a negative integer (when I'(a) or I'(b) = c) transforms a 
hypergeometric polynomial in ascending powers into one in descending powers. 
Applying (19) to (18), we get the well-known formula 


ven ieee ae l—-m = 
P(x) = DTNe mye (1 — x?) oF(- >" >-—_— 5 -1+4, 1/3") 


for the Legendre functions. 


We can also obtain this expression directly if we put the values (16) and (17) 
for our constants into the P-function 


P=EP-"(1—€) Fa +f’ +y,0' +R +y, 1 +f’ —B, 1/2). 


This is possible because either « +’ +y or «’ +8’ +y is in our case a negative 
integer and therefore the above P-function is regular in the fundamental interval 
given by €=x*=1. The fact that the same values of the constants which we used 


formerly are applicable also in case of our new P-function is due to its special 
choice. 


Sommerfeld’s Polynomial Method Simplified TE 


§4. RADIAL FUNCTIONS OF A NON-RELATIVISTIC ONE-ELECTRON 
ATOM 


This is an example of an infinite fundamental interval. Equation (1) is here 


given by 
= azl* Z) (4 a (E+ =)- oe |f- (eae. (20) 


Since « =0 is a singular point of (20) we must take the fundamental interval from 
€=0 to €=+0. 
A comparison of (20) with (1) gives 


p=p=x, A= B,  g=—2xja,+l+1), ...... (21) 


where a, =h?/me?Z is the radius of the first Bohr orbit. 
With regard to (12) we have therefore 
—U(L+1) +2x/a, +Ax? =h?[ — a(x —1/h) +(a—a’ +4 —1/2h)E — €2/4]. 


Hence in accordance with (3) 


n= Cay, 
and therefore 
Use a ented (seed @ enn eee anaes (22 a) 
Bid (Ga eet = ee Ea) oe, (22 b) 
Nate) Gate a) Oye nk Levee. ate (22 c) 


From (22a) it follows that «=/+1 so that, since a’=—n, we have 
a—a’=n+l+1. We get, therefore, according to (225), 


=e (tata dja aarti ele (23) 


For the eigenvalues of E we thus obtain, using (21) and (22), the well-known 
Balmer formula. The eigenfunctions of (20) are given according to (II), (9), 
a’ =—n, «=1+1 and h=1 by 

f=x'lexp (—Kx/2)F,(—n, 21 +2, xx) 


where x is given by (23). 

In order to get positive values for E and A the constant « must, by (22c¢) be 
imaginary. If we put «= +2zk then according to (22) we have «—a’= ¥1/ka, 
and hence using «=/+1 we have a'=/+1+i/ka,. We get therefore for the 
eigenfunctions f in the continuous eigenvalue spectrum the expression 


fewlexp (Fike) F,(1 +1 £i/ka,, 2142, +2ike), 


k/27 being the wave number at infinity. 
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Investigation of the Spectrum of the High Voltage Arc in 
Carbon Dioxide: the CO Flame Spectrum 
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ABSTRACT. The spectrum emitted by a high voltage arc in CO, at atmospheric pressure 
has been investigated in-the region 2000 a. to 9000 a. It was found to consist, mainly, of a 
continuum and the ‘ CO flame bands’. A study of the spectrum at different dispersions 
and on photographic plates of various kinds shows that although the Schumann-Runge 
O, bands do form a part of the spectrum they are relatively weak compared with the bands 
attributed to CO,. The excitation of the Schumann—Runge O, system in flames is discussed. 


Se UENZER OD) Cueto 


N earlier paper (Feast 1949) described a high voltage, low current arc 
and also discussed the use of this source in the investigation of the 
spectrum of the O, molecule. More recently the are source has been 

used in the investigation of other gases, viz. N,, H,, N,-H, mixtures, and NH3 
(Feast 1950). In continuation of these investigations the spectrum emitted by 
carbon dioxide in the arc source has now been studied and the results obtained 
~ are discussed in the present paper. Experiments on carbon dioxide are of 
particular interest at this time as there is some controversy over the question of 
the identity of the spectrum emitted by the CO-O, flame. This spectrum has 
usually been attributed to CO, (see Gaydon 1948 and references given there). 
Recently, however, it has been suggested that the major part, if not the whole, 
of this spectrum may simply be the emission Schumann—Runge O, band system 
(Hornbeck 1948, Hornbeck and Hopfield 1949). The fact that the Schumann— 
Runge bands were found to be so easily excited in the high voltage arc in oxygen 
at atmospheric pressure and to be so sensitive to the pressure of the gas (Feast 
1949) made a study of carbon dioxide in this source even more interesting. 


§2. EXPERIMENTAL 


The arc source used in the investigation of oxygen was again used in these 
experiments. ‘he arc was struck between wires of a platinum-—rhodium alloy 
in a litre ‘Pyrex’ flask containing the gas. Small dispersion spectrograms of 
the arc in CO, at atmospheric pressure were taken on medium contrast plates 
with a Hilger E.3 spectrograph. Apart from some weak CO Third Positive 
bands and a few lines (oxygen, carbon, platinum and rhodium) the spectrum 
consisted solely of an apparent continuum stretching from the ultra-violet region 
to the infra-red region. By using a prism-grating spectrograph it was possible 
to follow the continuum as far as 9000a. When photographic plates of higher 
contrast were employed it was found that in the region 3000a. to 5000a. there 
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was a large number of bands superimposed on the continuum. By a comparison 
with Gaydon’s published spectrograms they were identified as the CO flame 
bands. An investigation of the spectrum at various pressures showed that the 
flame bands decreased rapidly in intensity as the gas pressure was reduced and 
that they gave way to the normal low pressure discharge type CO spectrum. 
It should be noted that the appearance of the CO flame spectrum thus found is 
very different from the appearance of the Schumann—Runge O, bands 
photographed with the same spectrograph from the arc in oxygen. 

Following on the low dispersion work the region 3000 a. to 50004. was inves- 
tigated with a 10 ft. Eagle Mounting grating. Using high contrast (process) 
plates the Schumann—Runge bands were found and nothing else appeared to 
be present. These bands were easily identified by comparison with those obtained 
in the oxygen investigation. However, using low contrast plates (Ilford Zenith) 
the band lines were lost in continuum. 


§3. DISCUSSION OF OBSERVATIONS 


In order to discuss these observations it is necessary first to review briefly 
the present state of the work on the CO flame bands. ‘These bands have usually 
been studied in the spectrum of the CO-O, flame. A summary of all but the 
most recent work on these bands is given by Gaydon (1948). The conclusion 
reached by Gaydon and others was that the bands were to be attributed to a CO, 
electronic band system. However the analysis of the band system could not be 
made with sufficient accuracy for this assignment to be quite definite. Other 
possible emitters such as C,O, and CO; were considered unlikely whilst the 
spectrum appeared too complex to be assigned to any diatomic emitter. Quite 
recently, on the other hand, it has been suggested (Hornbeck 1948, Hornbeck 
and Hopfield 1949) that the major part if not all of the flame bands are to be 
attributed to the Schumann—Runge O, band system. High dispersion spectro- 
grams of CO-O, flames and explosions have revealed the emission Schumann-— 
Runge bands and also (in the infra-red region) emission bands belonging to the 
the same system as the O, atmospheric absorption bands. 

The present low dispersion spectrograms of CO, in the high voltage arc show 
little or no resemblance to the Schumann—Runge bands emitted by O, in the 
same source. Besides the other differences in the two spectra it is important to 
notice that though the rotational line structure of the Schumann—Runge bands 
is reasonably clear even on Hilger E.3 spectrograms, the flame bands are diffuse 
and without trace of fine structure. Hence if the Schumann—Runge system is 
emitted even weakly by the arc in CO, it will dominate the spectrum at high 
dispersion (especially with high contrast plates) for, as is well known, the apparent 
intensity of a continuum or quasi-continuum falls rapidly compared with the 
apparent intensity of a resolved line as the dispersion of the spectrum is increased. 
It should be noted that the differences between the high voltage arc spectra of 
O, and CO ,cannot be attributed to differences in rotational or vibrational intensity 
distributions in the two cases since the high dispersion work shows that the 
structure of the Schumann—Runge bands is the same in the two spectra. 

Gaydon and Wolfhard (1949) state that the CO flame spectrum is emitted from 
a different part of the CO—O, flame from the Schumann—Runge bands. Inthe arc 
however no such separation of the various spectra was observed. 
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§4. EMISSION MECHANISMS FOR THE CO FLAME AND THE 
SCHUMANN-RUNGE BANDS 


In attributing the CO flame bands to CO, Gaydon has suggested a mechanism 
for their emission. ‘The normal state of CO, is not formed from CO and O in 
their ground states CO(#X) and O(?P) but probably from CO(a*II) and O(P). 
Since there is no evidence for excited states of CO in flames, Gaydon concludes 
that the CO, is formed from normal CO(1X) and O(?P) giving an excited state 
of CO,. ‘The transition from this excited state to the ground state gives the flame 
bands according to this theory. Thus if M is a third body we have : 


CO+0-+M— €O;*-- Mb eee (1) 
CO;* CO, thy (banded). 7) eee (2) 


There is nothing in this theory to disagree with the present observations on 
the high voltage arc. Indeed the flame bands would be expected to occur quite 
strongly in the arc since this mechanism is very similar to the one used to explain 
the pressure dependence of the Schumann—Runge bands (Feast 1949). The 
fact that a similar pressure sensitivity has been found for the flame bands 
indicates a similarity in excitation mechanism and tends to strengthen both 
theories. 

It may be noted that the presence of the Schumann—Runge bands in flame 
spectra may be expected if the emission process for the system is considered. 
The emission mechanism suggested by the author involves a three-body collision 
resulting in the recombination of an O(?P) atom with a O(/D) atom into the 
upper state of the band system. Since it is known that there is a considerable 
_ concentration of oxygen atoms in flames (Gaydon 1948) it seems likely that some 
proportion of these will be in the O(:D) state. This state is metastable and 
hence there will be only a very small depopulation of the state due to radiative 
return to thegroundstate. Furthermore, the state is only 1-96 ev. above the ground 
state. ‘Thus the excitation of the system by this mechanism does not require 
the large amount of energy (6.09 ev.) needed if the system is directly excited and 
hence the appearance of the system in a source of low excitation such as a flame 
is more readily understood. This process is further supported by the fact that 
since Gaydon and Wolfhard (1949) find the Schumann—Runge bands in H,-O, 
flames also, it must be concluded that only oxygen is involved in the emission 
mechanism. 
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- Radioactive '!©Er 


There is no stable erbium isotope of mass number 165. Pool and Quill (1938) 
bombarded erbium with fast neutrons and found a positron activity with a half-life of 
1-1 min., for which a tentative assignment to !°Er was suggested. There was, however, 
no direct evidence for the chemical identity or mass number. 

To prepare ‘Er by the reaction '°Ho(p, n)!Er, 5 mg. of holmium oxide of over 
99-97°% purity was bombarded with 10 Mev. protons. Erbium carrier was then added and 
a chemical purification of the rare earth mixture from radioactive contaminants made. 
The erbium was then isolated by means of an ion-exchange column of the type described 
by Ketelle and Boyd (1947). 

The radioactivity of the erbium decayed entirely with a half-life of 10-0-40-1 hrs. The 
decay was measured, by means of a thin window (1:9 mg.cm~?) Geiger counter, over a 
period of ten half-lives. 

Experiments in a magnetic deflection device failed to indicate any electrons or positrons 
in the radiation, which was entirely electromagnetic. 

Measurements of the absorption of the electromagnetic radiation in aluminium showed 
two components with energies of 6:6 kev. and 52 kev. Absorption measurements of the 
harder component in copper gave a value of 47 kev. These values agree well with those 
for the L and K x-radiations of holmium. 'There was no harder gamma-radiation. The 
isotope 1®°Er therefore decays purely by orbital electron capture. 

The author wishes to thank the members of the Cyclotron Group for carrying out the 
irradiation. 


Atomic Energy Research Establishment, F. D. S. BUTEMENT. 
Harwell, Berks. 
7th May 1950. 


Kerete, B. H., and Boyp, G. E., 1947, 7. Amer. Chem. Soc., 69, 2800. 
Poor Mei, and Our, Ly Iv 1938. Phys. Rev, 53,437. 


Proton Capture Radiation of Fluorine-19 


When 1°F captures a proton of low energy, it is well known that the compound nucleus 
of 2°Ne which is formed has an excess energy of 13 to 14 Mev. which may be emitted in the 
form of gamma-radiation. The most probable process, however, for the decay of this 
excited nucleus is the emission of an alpha-particle. This process usually leaves the 
residual 16O nucleus in an excited state which subsequently emits gamma-radiation of about 
6 Mev. The yield of gamma-radiation due to the decay of ?°Ne to the ground state is likely 
to be small compared to that of 6 Mev. radiation, and Walker and McDaniel (1948), who 
searched for energetic gamma-rays with a pair spectrometer, were unable to detect any such 
hard radiation. They stated that this radiation, if it existed, had an intensity of less than 
0-3% of that of the 6 Mev. gamma-rays accompanying the alpha-decay for a thick target 
bombarded with protons of energy 1,150 kev. 

Devons and Hereward (1948) made a search for this hard gamma-radiation, using the 
photo-nuclear effect in ®*Cu, which has a threshold at about 11 Mev., and found that gamma- 
rays having energy sufficient to activate the Cu were emitted at a proton bombarding energy 
of 660 kev., a value agreeing with that for the short-range alpha-particle resonance. These 
authors, who were using a thick target, found that the intensity of the hard radiation increased 
only very slowly as the bombarding energy was increased from 700 to 960 kev., and certainly 
possessed no further detectable resonances in this region. They estimated the intensity of 
the radiation from the 660 kev. resonance, which they assumed to be of energy 13-4 Mev., 
corresponding to a transition to the ground state of °°Ne, to be 2:2 -+-0:8% of the intensity 
of the 6 Mev. radiation from that resonance. 

In view of the results of Devons and Hereward, we have used a pair spectrometer to 
search for hard gamma-rays from this reaction using a thinner target than that of Walker 
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and McDaniel, and a proton bombarding energy corresponding to the 660 kev. resonance. 
The u.T. set at Glasgow consists of an 800 kev. cascade generator, with an accelerator tube 
and a high-frequency ion source. Magnetically resolved proton beams of about 100 pa. 
were used. The targets had a stopping power of about 50 kev. and were prepared by 
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Figure 1. Gamma-rays from °F +p, 340 kev. resonance. 
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Figure 2. Gamma-rays from °F +p, 660-kev. resonance. 


evaporation of calcium fluoride on to a brass plate. Energy measurements were made with 
a pair spectrometer which is similar to that described by Walker and McDaniel, but which 
has ten counters placed in the field of a twelve-inch electromagnet. The electron pairs were 
produced in a lead sheet of thickness 0:0015 inch. Typical spectra obtained under these 
conditions are shown in Figures 1 and 2. At the 340 kev. resonance (Figure 1) it will be 
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seen that there is a strong gamma-ray component having an energy of about 6 Mev. together 
with a very weak component at about 7 Mev. Above this energy no further gamma-ray 
lines are evident. Figure 2 shows the spectrum of the radiation emitted at the 660 kev. 
resonance, and here it will be seen that the 7 Mev. component has become appreciable, and 
that there has also appeared a weak line having an energy which has been determined as 
12-0-+.0-2 Mev. by comparison with the *Li(p, y) radiation. The intensity of the 12 mev. 
line has been estimated as 2°% of that at 6 Mev., which is in good agreement with the results 
of Devons and Hereward. It should also be noted that this value is not in contradiction 
with the findings of Walker and McDaniel, since their value of less than 0-3°% of the intensity 
of the 6 Mey. radiation was determined for a bombarding energy of 1,150 kev. With this 
high bombarding energy and a thick target, the intensity of the 6 Mev. radiation would be 
about 15 times as great as that due to the 660 kev. resonance alone (Bonner and Evans 1948), 
while the intensity of the capture radiation is known to increase very little when the bombard- 
ing energy is increased from 700 to 960 kev. (Devons and Hereward 1948). For bombarding 
energies between 960 and 1,150 kev. there are two weak gamma-ray resonances which 
contribute only about 1:5°% to the total output of gamma-rays (Bonner and Evans), and so 
unless these two resonances produce a high proportion of capture radiation, it is unlikely 
that its intensity, for a bombarding energy of 1,150 kev., exceeds the figure of 0:3% given 
by Walker and McDaniel. 

The most interesting feature of the present experiment, however, is the determination of 
the energy of the capture radiation as 12:0 Mev. This suggests that the radiation is not due 
to a transition to the ground state of ?°Ne as supposed by Devons and Hereward, but is 
emitted during a transition to some intermediate level which is about 1:4 Mev. above the 
ground state as shown in Figure 3. There seem to be good grounds for this interpretation 
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Figure 3. Energy-level diagram, '°F +p. 


of the experimental results, for such a state has already been detected as a result of the study 
of other reactions (Curran and Strothers 1940, Powell 1943, Powell, May, Chadwick and 
Pickayance 1940). 
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26th April 1950. R. D. SMITH. 
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Isotope Effect in Superconductivity 


In a paper recently submitted for publication it has been shown that investigation of the 
‘ self-energy ’ of the free electrons of a metal in the vibrational field of the lattice leads to a 
quantatitive understanding of superconductivity. So far the theory has been developed for 
the absolute zero of temperature only ; there the energy difference between the superconducting 
and the normal state is 
C?mne) 


cM ’ 


|SJ=c 


where c,~0-1, m and M are the electronic and ionic mass respectively, € is the kinetic 
energy of the fastest of the 7, free electrons and C is an energy connected with the interaction 
between an electron and the displaced lattice. Using the empirical fact that apart from a 
numerical factor |S | is equal to the thermal energy of the free electrons in the normal state 
at the transition temperature 7, i.e. 


[S| ns €(RT 9) nel 


(c, is a numerical constant), it follows that 


For different isotopes of the same element C has the same value, so that RT) 0c 1/4/M. 
This is in excellent agreement with the experiments by Reynolds et al. (1950) and Maxwell 
(1950) on mercury which have just come to my notice. This agreement forms a direct 
check on the basic idea of the theory : interaction between electrons and lattice vibrations. 


Department of Theoretical Physics, H. FROHLICH. 


University of Liverpool. 
19th May 1950. 
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Excitation of Spectra in the Inner Cones of Flames 


Flames of most organic compounds premixed with air or oxygen show a relatively bright 
inner cone, corresponding to the reaction zone of the flame. The luminosity is no doubt 
partially confined to this region because active radicals emitting the light, such as C, and 
CH, are only present here, and also in some cases molecules may be formed in electronically 
excited states by the chemical processes. However, we have recently accumulated evidence 
showing that there is some active species present in the reaction zone which causes 
abnormally high excitation of inert atoms introduced into the flame. 

We have studied the effect by means of the spectrum-line reversal method, using a 
carbon arc with adjustable sector as background. Most of the work has been done on flames 
at very low pressure, because these can be made to give a thick flat reaction zone (Gaydon 
and Wolfhard 1948, 1949). When a metal vapour such as Fe (as iron carbonyl) or Pb (as 
lead ethyl) is introduced it is found that most of the lines of the spectrum are obtained in 
emission even against the background of the full light of the arc (around 3,800° k.), but that 
resonance lines may be reversed. The reversal temperatures are usually much nee than 
the expected flame temperature, and lines of short wavelength give higher temperatures than 
lines of longer wavelength. Results for various resonance lines for an acetylene—air flame 
at 24 mm, Hg pressure are tabulated, 
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: Reversal 7 
Tae ' d Reversal : 
temp. (" K.) ae as temp. (° K.) hee 
Na 5890-6 2010 pain Fe 2983°5 3050 
Fe 3859-9 2565 14-4 Fe 2966-9 apee ee 
Fe 3824-4 2580 14-45 Fe 2936-9 3100 15-7 
TI 3775-7 2350 1602 Pb 2833-0 2970 17-0 
Fe 3719-9 2650 14-5 Fe 2719-0 3310 15-9 
Fe 3679-9 2630 14-8 Fe 2522-8 3365 16-8 
Fe 3440-6 2880 14-4 Fe 2501-1 3410 16°75 
Fe 3020-6 3058 15-5 Fe 24833 3450 16-7 


The actual observation is of the effective temperature corresponding to the ratio of the 
populations of atoms in the ground and excited electronic states. This may be represented 
by the value of —E/kT, which is equal to the natural logarithm of the ratio of populations in 
the ground state and excited states; E is the energy of the excited state expressed in cm~, 
and k is the Boltzmann constant. In the Figure the values of E/RT are plotted against E. 
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For equilibrium the points should lie on a straight line passing through the origin, the slope 
giving the temperature. Clearly there is a major departure from equilibrium. The actual 
flame temperature from reversal measurements just above the reaction zone appears to be 
about 2,000° k. For energies below about 18,000 cm! the excitation may be mainly due 
to thermal causes, but at higher energy the population is much higher than could be accounted 
for thermally, and indeed the value of E/RT for the Fe lines varies surprisingly little with E. 
The extent of the abnormal excitation appears to vary from element to element, Pb and Tl 
being less highly excited than Fe. The low slope of the graph at high FE explains the 
extremely high apparent temperatures for non-resonance lines. The Tl green line at 
5350 a., which is not a true resonance line, reverses at 2,750° K.;_ by combining this value 
with that for Tl 3775, which has the same initial level, we obtain E/RT for the lower level of 
the 5,350 line at 7793 cm—. 

The amount of energy available for the abnormal excitation in these flames is surprisingly 
high. The emission lines of Fe 1 from the flames of C,H, C,H, and CH;0H have an 
intensity distribution very similar to that in an iron arc. Lines requiring an excitation 
energy up to at least 60,757 cm occur, corresponding to 174 kcal/mole. It is very difficult 
to find any chemical reaction liberating this amount of energy. No lines of Fe 11 have been 
observed. 

The effect is shown by the reaction zones of all the organic flames we have studied 
(Cy»H,, C.Hs, CH,, CH;OH) except formaldehyde. Flames of hydrogen and carbon 
monoxide do not show abnormally high excitation. Flames of organic compounds sup- 
ported by N.O show the effect, but not as strongly as flames with air or O,. The effect 
occurs at atmospheric as well as at low pressure, although we have not succeeded in making 
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quantitative measurements at 1 atmosphere; the measured reversal temperatures do not 
vary much with pressure in the range 10 to 200 mm. Hg. Diffusion flames of the type 
studied by Wolfhard and Parker (1949) at atmospheric pressure do not show abnormally 
high excitation; at very low pressure even diffusion flames show high excitation, probably 
because of the rapid mixing by diffusion under these conditions. In the pre-mixed flames 
the effect is most pronounced for slightly rich mixtures, but it decreases surprisingly slowly 
on the side of lean mixtures. 


Imperial College, A. G. GAYDON. 
London, S.W. 7. H. G. Wo.LFHARD. 
5th May 1950. 
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CORRIGEND: 


“The Angular Distribution of Synchrotron Target Radiation: A Preliminary 
Experimental Study ”’, by J. D. Lawson (Proc. Phys. Soc. A, 1950, 63, 653). 


Page 654, line 25 should begin \/2m/E instead of \/2p/E. 
Page 659, equation (1) should read 


EMCO exp (— E?6?/4402). 


EF? 
440nt 
Insert under Figure 1: 

The values of t shown on the diagram should be multiplied by 2. 


REVIEWS OF BOOKS 


lone Electrons and Ionizing Radiations, by J. A. CROWTHER. Pp. ix+322+7 
plates. 8th Edition. (London: Edward Arnold and Co., 1949.) 21s. 


The preface to the eighth edition of this well-known book states that it is intended to be 
“* a textbook from which students who have been grounded in the more elementary portions 
of physics might obtain a systematic knowledge of its later developments ”’; the level of the 
students is said to b2 that of the second year ot a pass degree course. The previous editions 
have introduced to a whole generation of students the subjects of Atomic and Nuclear Physics 
and the book has been deservedly popular with this past generation because of the clear style 
in which it is written and the interesting manner of presentation. These characteristics of 
Professor Crowther’s writings are, of course, present in the new edition, but the present 
reviewer considers that the book does not fulfil its stated purpose for the new, post-war, 
generation of students. 

One, to some extent, admirable feature of the book is the description of ‘ classical ’ 
experiments and apparatus. However, this is grossly overdone; it should be possible to 
retain the more important of these descriptions while giving accounts of modern apparatus 
and methods, which are just as easy to understand. ‘This would give a much better balance 
between a historical description of the subject and a logical one—that is, a better balance 
judged from the standpoint of our present knowledge of Atomic and Nuclear Physics. 
Another modification which is necessary in order to achieve this better balance is in the 
amount of space allotted to the various subjects. ‘Too much space is devoted to a subject 
like “ the conduction of electricity in gases’, whose importance is historical and technical, 
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and too little is given to such fundamental topics as ‘ the structure of the atom’ and ‘ the 
physics of the nucleus’. The present reviewer believes that the book would be improved 
if it were substantially revised so as to strike a new balance along the lines indicated above. 

There are also many serious errors and misleading accounts, particularly in the chapters 
on the nucleus and related subjects. Professor Crowther describes the Geiger counter for 
detecting single B-particles as a ‘ point’ counter. If the neutrino is to be discussed, it is 
surely possible in a book written in 1949 to give a brief account of the experiments of Sherwin 
and of Allen, rather than mention only the twelve-year-old cloud chamber work of Crane and 
Halpern, which Professor Crowther describes as recent. The section on the properties of 
the meson could have been written much better in 1949. It is certainly not correct to 
describe the {D(d, n)?He reaction as “ by far the most prolific source of neutrons’. 
Professor Crowther should have known that the 184 m. cyclotron at Berkeley has been 
working for two years now, and that the difficulties associated with variation of mass with 
velocity have been overcome by frequency modulation and not, as he implies, by shaping 
the pole pieces of the cyclotron magnet. Another error which is inexcusable is the state- 
ment that the half-value periods of RaC’, ThC’ and AcC’ are 10-8 second or less. It has 
been known for at least ten years that the half-value periods are 1-5 x 10-4 sec., 3 x 10~? sec., 
and 5 x 10~* sec., respectively. 

This is only a small selection of the errors and misleading accounts in which the book 
abounds. Most of them occur in the sections on more recent work and the net result is to 
give the impression that there has been no serious attempt to bring the book up to date. 
It is, therefore, with regret that the present reviewer feels that he cannot recommend the 
eighth edition of this book. J. DAINTY. 


Materiewellen und thre Interferenzen, by M. von Laue. Pp. 392. Second 
Edition. (Leipzig: Akademische Verlagschaft Geist and Portig K-G., 
1948.) 


Only a few copies of the first edition of Professor yon Laue’s book, published in 1943, 
reached this country. The second edition is therefore all the more welcome. ‘The 
development of theoretical aspects of electron diffraction is due, in no small degree, to 
Professor von Laue and his collaborators, and all physicists interested in this subject are 
greatly indebted to the author for his masterly and comprehensive treatise. 

The publication comes at a particularly appropriate time as experimental work is in full 
swing again after a lapse during the war years, and a survey of the oustanding problems will 
be a great help to the experimental physicist. It is also some ten years since the only other 
textbook on this subject was published, and an up-to-date account of all theoretical develop- 
ments will satisfy a definite need. 

The subject matter of the book can be divided into three main parts, an historical and a 
theoretical introduction, an account of the geometrical theory of the diffraction of electrons 
by crystals, and a study of the dynamical theory of diffraction. The diffraction of molecular 
rays and of neutrons is treated in a separate chapter. In the main the practical applications 
of electron diffraction are not dealt with—experimental results are considered only in so far 
as they have theoretical significance. ' 

The preliminary sections, apart from an excellent historical survey, are concerned chiefly 
with quantum mechanics (Schrédinger theory and Dirac theory of the electron), and in 
particular with the Born theory of elastic and inelastic collisions. It is rather striking that 
the theory of inelastic collisions should have so little quantitative information. to offer as to the 
probability of small energy losses, which are of particular interest, for example, in connection 
with the origin of Kikuchi lines or the intensity of the background in electron diffraction 
photographs. The position is unsatisfactory enough as regards atoms—where solids are 
concerned there seems to be no position at all. However, von Laue draws attention to some 
very interesting experiments on discrete energy losses of cathode rays after passing through 
thin films of solids, and it is from experiments of this type that one can hope to obtain the 
necessary information. 

The section on the geometrical theory of diffraction contains a complete account of 
reciprocal lattice theory, effects of structure factor etc. The most important feature is a 
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comprehensive discussion of the effect of the external form of a small crystal on the Structure 
of the intensity pattern in the reciprocal lattice. The theory of these ‘ shape factors , due to 
von Laue, has been of material assistance in the interpretation of electron diffraction patterns 
and will prove increasingly useful as the study of surface structure proceeds. The influence 
of crystal deformations on the diffraction pattern is also discussed, though somewhat briefly ; 
it seems likely that such features, hitherto rather neglected, would be well worth further 
consideration. 

Almost half the subject matter of the book is devoted to a discussion of the interaction of 
diffracted rays inside a crystal—the dynamical effects. The two-ray problem is worked out in 
detail, and a complete account is given of the many-ray problem. The use of the dynamical 
theory is well illustrated in a detailed discussion of the experiments of Kossel and Mollen- 
stedt, in which special types of interference patterns are obtained in transmission (using 
convergent beams) through thin films of single crystals such as mica. Apart from checking 
many theoretical points, these experiments lead to the determination of the coefficients of 
the Fourier expansion of the potential of the crystal. 

The last sections contain an account of the theory of Kikuchilines developed by von Laue, 
together with an interpretation of the features of the structure of Kikuchi lines and bands. It 
appears from the discussion that there are still little-understood features, including, for 
example, the relative intensity of bands at widely separated angles. 

From the pcint of view of experimental work on electron diffraction, the discussion of the 
mean electrostatic potential or ‘ inner potential ’ of a crystal is of particular interest. ‘There 
are numerous experimental measurements of this quantity. von Laue shows that, ignoring 
surface effects, the potential inside a crystal is not constant except where there are special 
symmetry properties, e.g. for cubic crystals. The potential is shown to be determined by 
the components of the quadrupole moment tensor of the cell unit. What is however not 
made clear, is how large the deviations from a constant value are likely to be. This inform- 
ation might be obtained from the diamagnetic susceptibility of a solid, in certain cases, as 
the same tensor is formally involved here. In the case of quartz, for instance, the dia- 
magnetic susceptibilities, parallel and perpendicular to the principal axis, do not differ by 
more than one per cent, and one would not expect appreciable variations of the ‘ inner 
potential ’. 

However, even in the case of cubic crystals, the surface distortions may, in von Laue’s 
opinion, produce appreciable electrostatic fields. In fact, one might draw the conclusion 
that the experimental measurements of ‘ inner potential ’ have very little meaning, and the 
discrepancies in published values of the ‘ inner potential’ for a given crystal are taken as 
evidence on this point. This is rather unfair to the experimental physicist, as the early— 
admittedly rough— values are given the same weight as later values. In particular, care 
has been taken in recent experiments using cathode rays to avoid conditions in which strong 
rays, other than the one used for the measurement, appear in the diffraction pattern—which 
covers one of the criticisms of the experimental work made by von Laue. All recent 
experimental values of the ‘ inner potential’ of rocksalt, for instance, agree extremely well. 
The points raised by von Laue, particularly those on surface effects, do however need 
careful consideration by experimental workers, and further work will be required to see how 
important they are in practice. 

Physicists concerned with electron diffraction will find Professor von Laue’s book an 
indispensable guide to the theory. It possesses, in addition, the very rare feature of 
suggesting experiments on almost every aspect of the theory. M. BLACKMAN. 
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ABSTRACTS FOR SECTION B 
The Sliding Surface (34th Guthrie Lecture), by G. I. FINcu. 


The Influence of Concentration-Dependent Diffusion on Rate of Evaporation, by 
J. CRANK. 


ABSTRACT. It has been suggested that evaporation through a membrane or from a 
polymer film has a characteristic behaviour different from that commonly associated with a 
free liquid surface. In particular it is claimed that an increase in surface humidity produces 
an increased rate of evaporation through a membrane, though the opposite effect is well 
known for a liquid surface. This behaviour is said to be observed when diffusion within the 
membrane or film is the rate controlling process in systems for which the diffusion coefficient 
increases with increasing concentration. 

In this paper the behaviour of such diffusion-controlled evaporation processes is examined 
from the point of view of the diffusion equations and it is found that there is no theoretical 
reason to expect such processes to show anomalous behaviour. 

It is seen that when the experimental data for evaporation through a membrane are 
correctly interpreted they are in agreement with theoretical prediction. 

The so-called “ skin effect’ in the dry spinning of fibres is examined and appears to be 
due to complicating features other than diffusion. 


The E.M.F.-Temperature Calibration of Platinum, 10°, Rhodium—Platinum and 
Platinum, 13°% Rhodium-Platinum Thermocouples over the Temperature 
Range 0°-1,760 °c., by C. R. Barer. 


ABSTRACT. E.M.F. temperature reference tables are given for platinum, 10° rhodium— 
platinum and platinum, 13% rhodium—platinum thermocouples over the range from 0° to 
1,750°c. ‘The tables are based on the calibration of twelve thermocouples of each kind, six 
from each of two manufacturers. Each set of six thermocouples comprised two samples 
from each of three batches of wire. The calibrations are in terms of the International 
Temperature Scale as revised at the Ninth General Conference of Weights and Measures 
held in Paris in October 1948. 
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The methods of intercomparison between the thermocouples and the various inter- 
polation instruments specified in the International Temperature Scale are described. It 
is shown that the differences between the new tables and the old are mainly due to changes in 
the International Temperature Scale, and that the new tables represent the shape of the 
E.M.F.-temperature characteristics of the thermocouples of either make sufficiently closely 
to give interpolation between calibration points with accuracy within = 11° Cup tol 14008 2e% 
+2° c. up to 1,550° c., and +3° c. above this temperature. 


A Joule-Thomson Cascade Liquefier for Helium, by J. ASHMEAD. 


ABSTRACT. A helium liquefier has been constructed with an output of 3-8 litres per 
hour and a starting time of 70 minutes. The liquid hydrogen needed to cool the helium 
below its inversion temperature is provided continuously by a hydrogen liquefier within the 
same machine. The advantages of the machine are (i) short starting time, (ii) simplicity of 
operation and ease of servicing, (iii) closed hydrogen circuit. 


Time Decrease of Magnetic Permeability in Alnico, by R. Street and J. C. WOOLLEY. 


ABSTRACT. The formal theory, based on the concept of thermal activation of domain 
processes, previously proposed to account for the phenomenon of magnetic viscosity, has 
been extended to describe time decrease of permeability. Within certain limits the pre- 
dictions of the theoretical analysis are verified by measurements made on time decrease of 
permeability using an alnico rod specimen. It is shown that the quantitative results of these 
measurements are in good. agreement with those obtained from the magnetic viscosity 
experiments. The possible physical mechanisms responsible for both magnetic viscosity 
and time decrease of permeability in alnico are discussed. 


The Magnetic Properties of Uranium and Uranium—Iron Alloys, by L. F. BATES 
and J. R. MALLARD. 


ABSTRACT. Measurements have been made of the magnetic susceptibility of fairly 
pure uranium and of 5 and 10 atomic per cent alloys of iron with uranium. Using a modifi- 
cation of the Curie method, with special pole-pieces and electromagnetic dynamometer 
control, the specimens were found to be paramagnetic with no trace of ferromagnetism in the 
temperature range 20° to 350°c. At 20°c. the susceptibility of pure uranium is 
+1-:740 x 10-* E.M.U. per gm., increasing to +1:804 x 10-* at 350° c.; the results fit the 
equation 


Xy=32 x 10 T+ 1-564 x 10+ 24-0 « 10-4 T. 


Assuming that the iron is present in form of the compound U,.Fe, it is found that the latter 
has a susceptibility of +2-036x10-* E.M.U. per gm. at 20° c., which remains practically 
constant with temperature. 


On Phase-Contrast Testing with a Sht Source, by E. H. Linroor. 


ABSTRACT. In the standard form of the phase contrast test, the light-source is a 
pinhole comparable in size with the Airy disc of the system under test. As a consequence 
it is difficult to obtain enough light through the pinhole to carry out the test with full 
efficiency. Ifa slit source and a phase-changing strip are used this difficulty is overcome, 
but the properties of the test are altered; for example, the brightness distribution seen on 
the surface of a true mirror under the test is no longer radially symmetrical. 

In the present paper the properties of the phase-strip test are examined in the mathe- 
matically simplest case where the light-source is a slit Of negligible width, passing mono- 
chromatic light, and it is concluded that the test can safely be used in figuring an optical 
surface to an accuracy of about one twentieth of a fringe. Below this limit, the systematic 
error inherent in the usual interpretation of the test begins to be appreciable in the case of 
es errors, though rapid local errors are correctly seen right up to the limit of sensitivity of 
the test. 
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